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Abstract. In this paper, we present a recent developments in the theory of transport 
phenomena based on the Fermi liquid theory. In conventional metals, various transport 
*y ■ coefficients are scaled according to the quasiparticles relaxation time, t, which implies 

that the relaxation time approximation (RTA) holds well. However, such a simple 
<^ ■ scaling does not hold in many strongly correlated electron systems. The most famous 

example would be high-T c superconductors (HTSCs), where almost all the transport 
r^ coefficients exhibit a significant deviation from the RTA results. This issue has 

I , been one of the most significant unresolved problems in HTSCs for a long time. 

Similar anomalous transport phenomena have been observed in metals near their 
antiferromagnetic (AF) quantum critical point (QCP). The main goal of this study is 
QJ , to demonstrate whether the anomalous transport phenomena in HTSC is the evidence 

of non-Fermi liquid ground state, or it is just the RTA violation in strongly correlated 

£NJ ■ Fermi liquids. Another goal is to establish a unified theory of anomalous transport 

phenomena in metals with strong magnetic fluctuations. For these purposes, we 

^_. ■ develop a method for calculating various transport coefficients beyond the RTA by 

\ S Q . employing field theoretical techniques. 

In a Fermi liquid, an excited quasiparticle induces other excited quasiparticles by 
collision, and current due to these excitations is called a current vertex correction 
(CVC). Landau had noticed the existence of CVC first, which is indispensable for 
calculating transport coefficient in accord with the conservation laws. Here, we 
develop a transport theory involving resistivity and Hall coefficient on the basis of the 
microscopic Fermi liquid theory, by considering the CVC. In nearly AF Fermi liquids, 

j^ ■ we find that the strong backward scattering due to AF fluctuations induces the CVC 

with prominent momentum dependence. This feature of the CVC can account for 
the significant enhancements in the Hall coefficient, magnetoresistance, thermoelectric 
power, and Nernst coefficient in nearly AF metals. According to the present numerical 
study, aspects of anomalous transport phenomena in HTSC are explained in a unified 
way by considering the CVC, without introducing any fitting parameters; this strongly 
supports the idea that HTSCs are Fermi liquids with strong AF fluctuations. Further, 
the present theory also explains very similar anomalous transport phenomena occurring 
in CeMIn5 (M=Co or Rh), which is a heavy- fermion system near the AF QCP, and 
in the organic superconductor k-(BEDT-TTF). 

In addition, the striking w-dependence of the AC Hall coefficient and the remarkable 
effects of impurities on the transport coefficients in HTSCs appear to fit naturally 
into the present theory. Many aspects of the present theory are in accord with 
the anomalous transport phenomena in HTSCs, organic superconductors, and heavy- 
fermion systems near their AF QCPs. We discuss some of the open questions for future 
work. 
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1. Introduction 

1.1. Relaxation time approximation (RTA) and current vertex correction (CVC) 

The investigation of transport phenomena in metals is very significant since we can 
extract a large amount of important information with regard to the electronic states of 
metals. In conventional metals with weak electron-electron correlation, various transport 
phenomena are governed by a single parameter, namely, the quasiparticle relaxation time 
t. That is, the RTA holds well in such conventional metals [1]. For example, resistivity 
p is proportional to r -1 , which is proportional to T 2 in conventional Fermi liquids. The 
Hall coefficient /?h is independent of r, and l/e|i?n| expresses the approximate carrier 
density n. The magnetoresistance Ap/po is proportional to r~ 2 . The thermoelectric 
power S is proportional to T for a wide range of temperatures. The signs of i?n and S 
represent the type of carrier (i.e., electrons or holes). The behaviors of these transport 
coefficients, which are referred to as Fermi liquid behaviors, are well explained by the 
RTA. 

In strongly correlated electron systems, however, transport coefficients frequently 
exhibit a prominent deviation from conventional Fermi liquid behaviors. Because of 
this, it is difficult to obtain information on the electronic states from the transport 
phenomena. For example, in high-T c superconductors (HTSCs), both i?H and S/T are 
very sensitive to temperature; in some compounds, sign changes occur with changes 
in the temperature [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. Evidently, a simple RTA does 
not work well for HTSCs. For example, l/e|i?n| is considerably greater than the 
electron density in optimally doped or under-doped systems. Prominent "non-Fermi- 
liquid-like behaviors" are also observed in other strongly correlated systems such as 
organic superconductors [13, 14, 15, 16, 17, 18] and heavy-fermion systems (/-electron 
systems) [19, 20, 21, 22]. Although these behaviors reveal essential information about 
the electronic states, only a little was understood about them for a long time. To obtain 
a significant amount of information from these wealth of treasures, we had to develop 
the theory of transport phenomena in strongly correlated Fermi liquids. 

The main objective of the present study is to investigate the non-Fermi-liquid- 
like transport properties in HTSCs, which have been intensively studied as one of the 
central issues of HTSC. Thus far, two different kinds of theoretical models have been 
actively investigated: in the first kind of models, non-Fermi-liquid ground states are 
assumed. For example, Anderson proposed that the elementary excitations are described 
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by spinons and holons. It is considered that holons, whose density is small in under- 
doped systems, are responsible for the anomalous transport [5, 23]. In the second kind 
of models, Fermi-liquid ground states are assumed. For example, transport coefficients 
were calculated based on the nearly antiferromagnetic (AF) Fermi liquid picture [24]. 
In the previous studies, the RTA has been frequently used [25, 26, 27]. However, both 
these models can explain only a limited number of experimental facts involved in HTSC: 
many other anomalous transport properties remain unsolved. To understand the true 
ground state of HTSC, it is highly desirable to solve the issue of anomalous transport 
phenomena in HTSCs in a unified way. 

The main aim of the present study is to explain that the rich variety of anomalous 
transport phenomena observed in HTSCs can be understood in a unified way in terms 
of the Fermi liquid picture. For this purpose, we have to take the CVC into account 
correctly, which is totally dropped in the RTA. In interacting electron systems, an excited 
electron induces other particle-hole excitations by collisions. The CVC represents the 
induced current due to these particle-hole excitations. The CVC is closely related to 
the momentum conservation law, which is mathematically described using the Ward 
identity [28, 29, 30, 31]. In fact, Landau proved the existence of the CVC, which is 
called backflow in the phenomenological Fermi liquid theory, as a natural consequence 
of the conservation law [28]. The CVC can be significant in strongly correlated Fermi 
liquids owing to strong electron-electron scattering. However, its effect on the transport 
phenomena has not been studied well until recently. 

In the present study, we discuss the role of CVC in nearly AF Fermi liquids, such as 
HTSCs, organic metals, and heavy-fermion systems near the magnetic quantum critical 
point (QCP). For example, the RTA for the highly anisotropic t^ model cannot explain 
the transport anomaly in HTSC. We find that the RTA is unreliable in the presence 
of strong AF fluctuations, because the prominently developed CVC entirely modifies the 
RTA results [32, 33, 34, 35, 36, 37, 38, 39]. In a Fermi liquid, the transport coefficients are 
described by the total current J^, which is expressed as a sum of quasiparticle velocity 
v k and CVC [40, 41, 42, 43, 44, 45, 46]. In the present study, we find that the J k in 
nearly AF metals shows anomalous ^-dependence due to the CVC. This is the origin of 
the transport anomalies in nearly AF metals. Based on the microscopic Fermi liquid 
theory, we investigate the important role of the CVC in i?n, ^p/po, S, and the Nernst 
coefficient (z/). Furthermore, we confirm this idea by performing a numerical study 
based on the FLEX+CVC theory. In this approximation, the Coulomb interaction U is 
the only fitting parameter. 

To demonstrate the physical meaning of CVC, we discuss the scattering processes 
between quasiparticles in an isotropic model [28]. In the RTA, it is assumed that the 
conductivity due to a quasiparticle at k is proportional to the mean free path v k r k . Since 
Tk oc T~ 2 in a Fermi liquid [28], the resistivity according to the RTA — pR- TA ^ T -i - 
is finite. However, this result is not true since the momentum conservation law ensures 
that p = in a spherical model in the absence of the Umklapp process [42]. The failure 
of the RTA originates from the assumption that the velocity of the excited quasiparticle 
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Figure 1. Decay process of a quasiparticle at k. After the relaxation time r, the 
quasiparticle at k collides with another quasiparticle at k' in the Fermi sea. As a 
result of this collision, the two quasiparticles at k + q and k' — q, and a quasi-hole at 
k' are created. The total momentum is conserved in this process. 



disappears after scattering. Here, we explain that the correct answer (zero resistivity) 
in the absence of the Umklapp processes is recovered by considering all the relevant 
normal scattering process as shown in Fig. 1. When a quasiparticle at k is scattered 
to k + q after the relaxation time Tk, a particle-hole pair (at k' — q and k') should be 
created according to the momentum and energy conservation laws. The CVC represents 
the current conveyed by the particles at (k + q, k' — qj and a hole at k' ; which emerge 
during the scattering process. Therefore, the momentum and energy conservation laws, 
which are violated in the RTA, are restored by considering the CVC [42]. The CVC is 
necessary to reproduce the zero-resistivity in the absence of the Umklapp process. 

In the RTA, we consider an excited quasiparticle at k (in Fig. 1) as if it annihilates 
after the relaxation time. Since the RTA allows such an unphysical process, the 
conservation law dp/dt + V • j = is violated in the RTA. Since the conservation 
law is a very important constraint on the transport properties, the RTA frequently 
yields unphysical results [42, 45]. In later sections, we find that the CVC is significant 
in nearly AF metals: we explain that the CVC is the origin of a variety of anomalous 
transport phenomena in such metals. 



1.2. N on- Fermi-liquid-like transport phenomena in high-T c cuprates 

In HTSCs, almost all the transport phenomena deviate from the conventional Fermi 
liquid behaviors, which are referred to as the non- Fermi liquid behaviors [47, 48]. 
These anomalous transport phenomena have been studied intensively as one of the most 
important issues in HTSCs, since they offer important clues to reveal the true ground 
state in HTSCs, which has been unsolved until now. To answer this question, many 
analytical and numerical studies have been performed to determine the ground state 
of two-dimensional (2D) systems. For example, Anderson et al. considered that the 
ground state of a square-lattice Hubbard model is the resonating-valence-bond (RVB) 
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Figure 2. Schematic phase diagram of HTSC. AF spin fluctuations start to increase 
below To. At the same time, Rn starts to increase. Below T* (in the pseudo-gap 
region), the AF fluctuations are suppressed since the strong SC fluctuations reduce 
the density of states at the Fermi level, which is called the pseudo-gap. In this study, 
the filling area is |1 — n\ > 0.1. 



state, where spin and charge degrees of freedom are separated [23]. In the RVB state, 
concept of quasiparticle is not valid. According to his philosophy, possibility of the RVB 
state in the square-lattice t-J model was studied, by using the mean- field theory [49] 
and the gauge theory [50]. In §1.4, we will introduce significant works performed in 
the t-J model. Unfortunately, it is difficult to perform quantitative studies of transport 
coefficients using the t-J model. 

On the other hand, the Fermi liquid theory has been developed and applied to 
analyzing the HTSCs, considering that the Fermi liquid ground state is realized [51, 52]. 
By using the self-consistent renormalization (SCR) theory, Moriya, Takahashi and Ueda 
explained that the correct superconducting (SC) order parameter d x 2- y 2 as well as T- 
linear resistivity are derived from the strong AF fluctuations [53, 54, 55]. Based on the 
phenomenological spin fluctuation model, Monthoux and Pines performed a quantitative 
analysis for optimally-doped YBCO [56]. Further, Bickers et al. studied the square- 
lattice Hubbard model according to a self-consistent random-phase-approximation 
(RPA), which is now called the fluctuation-exchange (FLEX) approximation [57]. 
These spin fluctuation theories have succeeded in reproducing various non-Fermi-liquid 
behaviors in the normal state of HTSCs: For example, temperature dependence of 
nuclear spin-lattice relaxation rate given by NMR/NQR measurements, \jT{T oc 
Xl q I m Xq(^) A^U=o ^ T _1 , is reproduced well. Moreover, famous T-linear resistivity 
in HTSC, 

p oc 1/r oc T, (1) 

is also explained. (The Fermi liquid behavior 1/r oc T 2 will recover at very low T 
if we suppress the superconductivity.) The FLEX approximation is also useful for the 
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Figure 3. Experimental Hall coefficient and the thermoelectric power for LSCO and 
NCCO. [Rcf. [2]] 



study of electron-doped systems [33, 58, 59]. These spin fluctuation theories successfully 
reproduce various non-Fermi-liquid-like behaviors in HTSCs except for the under-doped 
region |1 — n\ < 0.1, where n is the number of electrons per site. A schematic phase 
diagram of HTSC is shown in Fig. 2. The Fermi liquid description of HTSCs is still 
in progress [60, 61, 62, 63]. A rf x 2_ y 2-wave symmetry in the SC state was confirmed by 
phase-sensitive measurements [64, 65, 66], and tunneling spectroscopy [67]. rf x 2_ y 2-wave 
SC state was also derived according to the third-order-perturbation theory with respect 
to U [68]. 

Regardless of the remarkable success of the Fermi liquid theory during the 90s, the 
non-Fermi-liquid-like transport phenomena observed in HTSCs remained unresolved for 
a long time. The significant deviation from the RTA results in HTSCs was frequently 
considered as a hallmark of the breakdown of the Fermi liquid state. This issue has 
remained as one of the most important problems in HTSCs for a long time. In HTSCs, 
|i?n| increases below T ~ 600K as 

Rr oc 1/T (2) 

above the pseudo-gap temperature T* ~ 200K, and |i?n| ^ 1/ne at T* [2, 69]. The 
sign of i?H is positive in hole-doped systems such as La 2 _5Sr5Cu0 4 (LSCO) [2, 3], 
YBa 2 Cu 3 7 _ 5 (YBCO) [4, 5] and Tl 2 Ba 2 Cu0 6+5 (TBCO) [6]; however, it is negative in 
electron-doped systems such as Nd 2 _5Ce,5Cu04 (NCCO) [7] and Pr 2 _5Ce,5Cu04 (PCCO) 
[8], even though the angle resolved photoemission (ARPES) measurements resolved hole- 
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like Fermi surfaces [70] . The experimental T-dependences of i?n for LSCO and NCCO 
are shown in Fig. 3. The magnetoresistance of HTSCs also shows strong temperature 
dependence as 

Ap/po oc T" 4 (3) 

for a wide range of temperatures in LSCO [12, 71, 72], YBCO [71] and TBCO [73]. 
These results completely contradict with Kohler's rule (Ru occonst. and Ap/po oc p^" 2 ) 
which is derived using the RTA for a single-band model. Interestingly, the following 
"modified Kohler's rule" holds well for optimally-doped LSCO [12, 71], 90K and 60K 
YBCO [71] and TBCO [73]: 

Ap/po oc tan 2 H , (4) 

where Oh = o~ xy /o~ xx is the Hall angle. This experimental fact strongly suggests that the 
anomalous behaviors of the Hall effect and the magnetoresistance have the same origin. 
Below T*, i?H decreases whereas Ap/p increases further [74]. Therefore, modified 
Kohler's rule is not very well applicable for under-doped HTSCs at low temperatures 
[75, 76, 77]. 

For a long time, anomalous transport phenomena have been considered as one of 
the strongest objection against the Fermi liquid description of HTSCs. For example, 
to explain eqs. (l)-(4), Anderson introduced the Tomonaga-Luttinger model with two 
types of relaxation times. However, it is not obvious how to describe the crossover 
from the Tomonaga-Luttinger liquid state in the under-doped region to the Landau- 
Fermi liquid state with doping. On the other hand, the Fermi liquid description for 
HTSC, starting from the well-established Fermi liquid state in the over-doped region, 
appears to account for a wide doping range in HTSCs [52]. Here, we conform to the 
following principle: before abandoning the Fermi liquid picture in HTSCs, we have to 
verify whether the RTA is really applicable in HTSCs or not. We stress that the RTA 
may be unreliable for strongly correlated Fermi liquids since the CVC is not considered, 
regardless of its importance to satisfy the conservation laws. Due to this inadequacy, 
the RTA sometimes yields unphysical results in correlated metals [42] . In later sections, 
we will explain that the non-Fermi liquid-like behaviors in HTSCs can be understood 
based on the Fermi liquid theory by including the CVC. 

Here, we discuss the pseudo-gap phenomena in slightly under-doped systems below 
T* ~ 200K, which is also one of the most important issues in HTSC. \jT\T starts 
to decrease below T*, which means that the AF fluctuations are suppressed in the 
pseudo-gap region [78, 79, 80, 81, 82]. According to ARPES measurements, prominent 
deep pseudo-gap appears in the density of states (DOS) below T* [83, 84, 85, 86]. 
A simple spin fluctuation theory cannot explain the various anomalous phenomena in 
the pseudo-gap region. Recent theoretical studies using the T-matrix theory [87] and 
the FLEX+T-matrix theory [88, 89, 51, 90] have shown that the strong SC amplitude 
fluctuations, which are induced by the AF fluctuations, are a promising candidate for the 
origin of pseudo-gap [51, 52, 89]. In §5, we study transport phenomena below T* using 
the FLEX+T-matrix theory by including CVCs, and show that the various anomalous 
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transport coefficients are well reproduced in a unified way. The present study strongly 
supports the idea that the pseudo-gap phenomena in under-doped HTSCs are induced 
by the strong SC amplitude fluctuations with a d x 2. y 2 -symmetry [52, 51, 89]. 

We note that wide and shallow pseudo-gap in the DOS (weak pseudo-gap) 
is observed by ARPES even above T*, which is considered to originate from AF 
fluctuations that appear below T ~ 600 K. We will discuss the weak pseudo-gap 
phenomena in §9.2. 

1.3. Non-Fermi liquid transport phenomena in CeMIn^ (M=Co,Rh, or Ir) and 
k-(BEDT-TTF) 

During the last decade, it has been found that in strongly correlated materials, 
including heavy-fermion systems and organic metals, various transport coefficients 
exhibit striking deviations from the Fermi liquid behaviors. In particular, anomalous 
transport properties similar to those in HTSCs (eqs. (l)-(4)) have been observed 
in many systems with strong magnetic fluctuations such as k- (BEDT-TTF) 2 X [ 
X=Cu[N(CN) 2 ]Br [13], X=Cu[N(CN) 2 ]Cl [14, 15], X=Cu(NCS) 2 [16, 17] ], /t-(BEDT- 
TTF) 4 Hg 2 . 89 Br 8 [18], and CeMIn 5 (M=Co,Rh) [19, 20]. BEDT-TTF is an abbreviation 
of bis(ethylenedithio)tetrathiafulvalene. These experimental facts strongly suggest that 
the transport anomaly given by eqs. (l)-(4) is not a problem specific to HTSCs, but 
a universal property in the nearly AF Fermi liquids [20]. The study of these transport 
phenomena in such systems will serve to resolve the origin of the non-Fermi-liquid-like 
behaviors in HTSCs. 

k- (BEDT-TTF) is a layered organic compound made of BEDT-TTF molecules. 
The c?-wave superconductivity can be realized in a wide region of the pressure- 
temperature (P-T) phase diagram, adjacent to the AF insulating states [91]. For 
example, Cu[N(CN) 2 ]Cl salt at ambient pressure is an AF insulator, with its Neel 
temperature T N = 27 K. With increasing pressure, T N decreases and superconductivity 
appears via a weak first order transition; The maximum T c is 13 K at 200 bar. An 
effective theoretical model can be given by the anisotropic triangular lattice Hubbard 
model at half-filling [92]. According to this model, the phase diagram of /t-(BEDT- 
TTF) can be well reproduced by using the FLEX approximation, by assuming that 
^eff/H^band decreases by pressure [93, 94, 95]. These studies revealed that the d x 2_ y 2 -wave 
superconductivity occurs due to the strong AF fluctuations q ~ (ir, tt). In the metallic 
phase of k- (BEDT-TTF), the relationships i?n oc T _1 and tan0H oc T~ 2 are observed 
[13, 14, 15, 16, 17, 18]. Moreover, the magnitude of Rn decreases with increasing the 
pressure [15, 18, 17]. These behaviors are quantitatively reproduced by the FLEX+CVC 
theory [96]. 

CeMIns is a quasi two-dimensional (2D) heavy-fermion compound [97, 98, 99]. At 
ambient pressure, CeCoIn 5 is a superconductor with T c = 2.3 K. The electronic specific 
heat coefficient 7 has been measured to be 300 mJ/K 2 mol at T > T c . Celrlns is also a 
superconductor at ambient pressure with T c = 0.4 K and 7 = 680 mJ/K 2 mol. CeRhIn 5 
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is an AF metal with Neel temperature T^ = 3.8 K at ambient pressure. Under the 
pressure, CeRhIn 5 undergoes a SC transition at P c = 2 GPa, indicating that the AF 
quantum critical point (QCP) is located at or in the vicinity of P c . The NMR relaxation 
rate T^ 1 measurements indicate the presence of quasi 2D AF spin fluctuations in the 
normal state [100, 101]. The measurements of the angle resolved thermal conductivity 
[102, 103] and specific heat [104, 105] revealed that the symmetry of the SC state exhibits 
rf-wave symmetry. Furthermore, the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) SC state 
has been observed [106]. 

In CeMIn 5 , eqs. (l)-(4) are well satisfied for a wide range of temperatures, and 
the value of \Ru/ne\ reaches ~ 50 in CeRhlns near the AF QCP [19, 20]. Similarly, 
the Nernst coefficient v also takes a large positive value at low temperatures [107]. 
Recently, we studied both i?n and v in a quasi 2D Hubbard model, and found that 
they are prominently enhanced, as large as those in pure 2D systems [108]. In CeMIns, 
modified Kohler's rule (4) is well satisfied for < H < 3 Tesla, whereas both o xy j ' H 
and (Ap/po)/H 2 are drastically suppressed by a very weak magnetic field (H > 0.1 
Tesla) near the AF QCP [19, 20]. This surprising fact can be understood in terms of 
the field-dependence of the CVC. 

1.4- Fermi liquid or non-Fermi liquid? 

In the present article, we assume that the Fermi liquid state is realized in a wide range 
of the phase diagram in HTSCs, and the non-Fermi-liquid-like behaviors are created 
by strong spin (and SC) fluctuations. According to this idea, we study the anomalous 
transport phenomena in HTSC, by carefully analyzing the main-body effects that had 
been overlooked in previous studies. Especially, we intensively study the CVC based on 
the microscopic Fermi liquid theory. Our final aim is to explain the anomalous transport 
phenomena in various nearly AF Fermi liquids in a unified way, including HTSC, heavy 
fermions and organic metals. 

However, it is a nontrivial question whether the Fermi liquid state is realized in 
two-dimensional strongly correlated electron systems near the half-filling. In fact, the 
Fermi liquid state seems to be broken in heavily under-doped HTSCs. In the infinite 
dimension Hubbard model, the Fermi liquid state with heavy mass is realized next to 
the Mott insulating state [109, 110]. In 2D Hubbard models [111], it was rigorously 
proved that the limit value of the interaction U, below which the Fermi liquid state 
is realized, is finite. However, in strongly correlated 2D systems, an exotic non-Fermi 
liquid ground state may be realized next to the Mott insulating state due to the strong 
quantum fluctuations. That is, removal of large tracts of the Hilbert space due to strong 
correlation effect may lead to a violation of the Fermi liquid state. 

To study the strong correlation effect in HTSCs, the t-J model has been frequently 
analyzed. The t-J model is derived from the Hubbard model (or d-p model) by a 
canonical transformation, by excluding the double occupancy of holes to represent the 
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strong Coulomb interaction. It is given by 

ii.ii 
H t-J = J2 P G (t hJ clc, a + h.c)P G + J J2 Si ■ S J (5) 

{i,j),<f (i,j) 

where Pq represents the exclusion of the doubly occupied state, tij is the hopping 
integral between (i,j) sites, and J is the superexchange energy between the neighboring 
spins. J = 0.10 ~ 0.14 meV in real HTSCs. In the t-J model, the Hilbert space with 
high-energy state is eliminated by Pq, which enables us to perform numeral calculations 
easier or to invent new approximations. Based on the t-J model, various versions 
of the new fluid have been proposed, with novel kinds of excitation, many involving 
gauge theories of spin-charge separated spinons and holons [23, 49, 50]. The exact 
diagonalization technique [112, 113] has been applied to the square-lattice t-J model. 
The ground state phase diagram of the t-J model has been studied using the variational 
Monte Carlo method [114, 115], and it was found that the d^.^-wave SC state is realized 
in a wide range of the phase diagram. Unfortunately, quantum Monte Carlo (QMC) 
simulation for t-J model is difficult because of a serious negative sign problem. Instead, 
QMC simulation for the Hubbard model with moderate value of U [116, 117] have been 
performed intensively. 

On the other hand, many authors have considered non-Fermi liquid ground states 
due to novel quantum criticalities, other than a conventional (SCR-type) quantum 
criticality near the spin density wave (SDW) state [61, 118]. For example, Varma et 
al. proposed that the marginal Fermi liquid state is realized in HTSCs, where the 
k-independent self-energy is given by [118] 

£(u; + iS) = \(u\n(uj c /x) — ix), (6) 

where x = max(|u;|, irT), A is a coupling constant, and u c is a cutoff energy. This state is 
not a Fermi liquid since the quasiparticle renormalization factor z = (1 — <9ReS/(9cj) _1 = 
(1 + \ln(uj c /x))~ l vanishes logarithmically as (w,T) — > 0. The self-energy in eq. (6) can 
be derived if electrons couple to the following k-independent charge and spin fluctuations 
that are singular at T = 0: 

lmP(co + iS) cxmin(|w|/T,l). (7) 

In this model, p oc — ImE(i5) = XnT, \jT{T oc ImP(u) /lj\ cj=0 oc 1/T, and the optical 
conductivity is cr(u) oc (uo — iXx)^ 1 . They are typical non-Fermi liquid behaviors 
in HTSCs. Now, a microscopic derivation of the k-independent quantum critical 
fluctuations is an important issue. Varma proposed [119] that the circulating current 
phase exists in the under-doped regime, and the current fluctuations near the QCP (~ 
at the optimum doping) are the origin of the marginal Fermi liquid state. 

In the present article, we will argue that transport anomaly near the AF-QCP can 
happen even if the Fermi liquid state (1 > z > 0) remains intact at the QCP. We will 
show that the CVC in the Landau Fermi liquid theory causes various striking quantum 
critical behaviors. 
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2. Spin fluctuation theory and model Hamiltonian 

2.1. Phenomenological spin fluctuation model 

Here, we discuss the functional form of the dynamical spin susceptibility Xq( u ) m 
nearly AF metals. Hereafter, we use the unit c = H = k-Q = 1. This is the most 
important physical quantity in such metals since it is the origin of various non-Fermi 
liquid behaviors in HTSCs. The phenomenological form of Xq( u )i which can be obtained 
by using NMR/NQR spectroscopy and the neutron diffraction measurement, is given 
by [53, 56, 120, 25] 

4H = ?i + a F (q- X Q) 2 --M f ' (8) 

where Q = (±7r, ±7r) is the antiferromagnetic (AF) wavevector, and £af is the AF 
correlation length. This is referred to as the Millis-Monien-Pines model [120]. In HTSCs 
above the pseudo-gap temperature T*, both xq an d l/^ s f are scaled by £^f as follows 
[121]: 

el F ^a o /(T + 0), (9) 

Xq « «i • £af> 1/wrf « a 2 ■ £ AF , (10) 

where 6, a , a.\ and a 2 are constants. Since XQ^sf oc £ AF in eq. (10), the corresponding 
dynamical exponent z is 2. The coefficient «o rapidly increases in the under-doped 
region: £af reaches ~ 2 a in optimally-doped YBCO, and it exceeds 100 a in slightly 
under-doped NCCO just above T c [122]. (a is the unit-cell length; we put a = 1 
hereafter.) The relationship u 8 t > T (u s f < T) is satisfied in the over-doped (under- 
doped) YBCO. On the basis of the phenomenological x S q( u ) model, the SC transition 
temperature was successfully reproduced by solving the strong-coupling Eliashberg 
equation [56]. 

Theoretically, the relationships in eqs. (9) and (10) can be explained according to 
the SCR theory, where we consider the renormalization of the dynamical susceptibility 
due to both self-energy correction and vertex correction. This renormalization 
effect is referred to as the "mode- mode coupling effect" [123], which represents the 
destruction of the AF long-range order due to thermal and quantum fluctuations. 
Relationships (9) and (10) are also reproduced according to the FLEX approximation 
[57, 124, 125, 126, 127, 128, 129]. In this approximation, only the mode-mode coupling 
effect due to the self-energy correction is taken into account: Since the large imaginary 
part of the self-energy reduces the DOS, the spin susceptibility given by the mean- 
field approximation (which is equivalent to the RPA) is drastically suppressed in the 
FLEX approximation. As a result, the magnetic long-range order does not occur in 
pure 2D systems, which means that the Mermin- Wagner theorem is satisfied in the 
FLEX approximation [130]. In nearly AF metals, the FLEX approximation and the 
SCR theory yield similar results since the mode-mode coupling effect due to the self- 
energy is dominant [123] . On the other hand, the mode-mode coupling effect due to the 
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second-order vertex correction (Aslamazov-Larkin term) can be important in the nearly 
ferromagnetic metals. 

2.2. Model Hamiltonian and FLEX approximation 

In the main part of this study, we investigate the transport phenomena in HTSCs based 
on the Fermi liquid theory, and perform a numerical study using the spin fluctuation 
theory. In this study, we analyze the following Hubbard model: 

H =J2 £ k C L C ka + U J2 C k+qT C k'-q| C k'|C kT , (11) 

kcr kk'q 

where U is the Coulomb interaction, and e£ is the spectrum of the conduction electron. 
In a square lattice, e£ is given by 

e£ = 2to(cos k x + cos k y ) + 4ti cos k x cos k y + 2£ 2 (cos2/c;e + cos2/c ?/ ), (12) 

where c ka is the creation operator of an electron with momentum k and spin a, and 
U is the on-site Coulomb repulsion. We represent the electron filling by n, and n = 1 
corresponds to the half-filling. To fit the band structures given by the local density 
approximation (LDA) band calculations for YBCO [131], NCCO [132], and LSCO 
[133, 134] and by the angle resolved photoemission (ARPES) experiments for YBCO 
[135], NCCO [70] and LSCO [136], we select the following set of parameters [137, 33]. (I) 
YBCO (hole-doping) and NCCO (electron-doping): t = — 1) h = 1/6, and £ 2 = —1/5. 
(II) LSCO (hole-doping): t = — 1, t\ = 1/10, and £ 2 — —1/10. Here, n is smaller 
(larger) than unity in YBCO and LSCO (NCCO). The Fermi surfaces for YBCO and 
NCCO without interaction are shown in Fig. 4 (i). The deformation of the Fermi surface 
in the presence of U has been discussed in Ref. [33]. Note that the Fermi surface in 
Bi 2 Sr 2 CaCu 2 8 (BSCCO) is similar to that in YBCO [138]. Since \t Q \ ~ 4000K in 
HTSCs, T = 0.1 in the present study corresponds to 400K. 

Here, we study the self-energy, dynamical spin susceptibility and the CVC using 
the FLEX approximation [57]. The FLEX approximation is classified as a conserving 
approximation whose framework was constructed by Baym and Kadanoff [139] and by 
Baym [140]. For this reason, we can calculate the CVC without ambiguity, by following 
the Ward identity T 1 = STj/SG. Here, the Green function and the self-energy are given 

by 

G k (e n ) = (ie n + \x - e£ - S k (e n )) _1 , (13) 

S k (e„) = T J2 C k -q(en - w,) • V q (w,), (14) 

VqM) = U 2 Qxq(^) + \X\{UI) ~ 4M) + U, (15) 

4 c) (^) = xl-{l-(+)Uxl(^)y\ (16) 

Xq (^) = -T J2 G q+k (^ + e n )G k (e n ), (17) 

k.rt 



Anomalous transport phenomena in Fermi liquids with strong magnetic fluctuations 14 



BZ bound 




(7t,3l) 



YBCO 
NCCO 




1 V"? N 

o: hot spot \ 5 
o : cold spot ) |_ 



(-3T-X) 



Elk, w) 



(i) 



AF fluctuations 



(JI-7I) 



k+q 

m 



Figure 4. (i) The Fermi surfaces for YBCO (n < 1) and NCCO (n > 1). The location 
of the hot spots and the cold spots are shown. According to the FLEX approximation, 
the hot spot in YBCO shifts to point B, by reflecting the large DOS at (ir,0). (ii) 
Diagrammatic representation of the self-energy in the one- loop (FLEX) approximation. 



21ttT, respectively. The self-energy in eq. (14) 
ii). We solve the eqs. (14)-(17) self-consistently, 



T Ekr,. G k(enJ- ! ' 



i-6 



where e n = (2n + 1)ttT and ui - 
is schematically shown by Fig. 4 
choosing \i so as to satisfy n — j. / , k n 

The FLEX approximation is suitable for the analysis of nearly AF Fermi liquids. 
Many authors have applied this approximation to the square-lattice Hubbard model 
[57, 124, 125, 126]. Although it is an approximation, the imaginary-time Green function 
obtained from the FLEX approximation is in good agreement with the results obtained 
from the quantum Monte Carlo simulations for a moderate value of U [57]. It has 
also been applied to the SC ladder compound, Sri 4 _ a; Ca a ;Cu2404i [127], and the organic 
SC k-(BEDT-TTF) compounds [93, 94, 95, 141]. Using the FLEX approximation, T c 
was studied in various types of tight-binding models [142] . The FLEX approximation 
predicts that the AF fluctuations are predominant in a square-lattice Hubbard model 
near half-filling. This result has been confirmed by renormalization group analyses, 
which considers both spin and charge fluctuations on the same footing by solving the 
parquet equation [143, 144, 145]. 

The FLEX approximation can not reproduce the pseudo-gap behaviors below 
T* ~ 200K in slightly under-doped systems. However, they are well explained 
by the FLEX+T-matrix approximation [88, 51, 89], where the self-energy correction 
due to strong SC fluctuations, which are induced by AF fluctuations, are considered 
self-consistently. We will explain about the FLEX+T-matrix approximation in §5. 
In the present study, we perform numerical studies from the over-doped region to 
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Figure 5. x(q,w = 0) for (i) YBCO (n = 0.85) at T = 0.02 and (ii) for NCCO 
(n = 1.20) at T = 0.04, given by the FLEX approximation. In YBCO, £af = 2a - 3a, 
whereas £af for NCCO exceeds ~ 10a at T = 0.02. a is the lattice spacing. [Ref. [33]] 



the slightly under-doped region (i.e., n < 0.9 or n > 1.1), where the FLEX(+T- 
matrix) approximation yields reasonable results. Note that the FLEX approximation is 
inappropriate to describe the "Mott physics" in the heavily under-doped region. 

Figure 5 shows the spin susceptibility Xq( u = 0) given by the FLEX approximation, 
both for YBCO (n = 0.85; optimum doping) at T = 0.02 and for NCCO (n = 1.20, 
slightly over-doping) at T = 0.04, respectively [33]. Since the nesting of the Fermi 
surface is not good as shown in Fig. 4, £af oc Jxq is moderate even at low temperatures. 
In optimum YBCO, the Stoner factor Ux%(0) ~ 0.98 at T = 0.02, and the AF 
correlation length £af is approximately 2 ~ 3 a (a denotes the lattice spacing). On 
the other hand, £af in NCCO (n = 1.20) exceeds 10 a at T = 0.02. It should be 
noted that Xq( u = 0) given by the RPA shows an incommensurate structure, which is 
inconsistent with the experimental results. I verified that \jT{T of the Cu nuclei under 
H _L c given by the FLEX approximation for n = 0.85 is 1 ~ 2 [1/K msec] at 200K, 
which is consistent with the experimental results [146]. 

One of the advantages of the FLEX approximation is that the Mermin- Wagner 
theorem with respect to the magnetic instability is satisfied: its analytic proof is given 
in the appendix A of Ref. [130]. Hence, the critical region (Uxq(0) > 0.99) is stable 
in 2D systems since the SDW order {Uxq{0) — l) is prevented by the Mermin- Wagner 
theorem. That is, U CI = oo in 2D systems. As a result, the U dependence on Rn (and 
Xq) given by the FLEX+CVC approximation is rather moderate, as shown in Refs. 
[33, 96]. 

In the next stage, we derive the two-particle Green function L(k, k') which is 
indispensable for the study of transport phenomena. According to the microscopic 
Fermi liquid theory [30, 31], L(k,k') can be obtained by the solution of the following 
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Figure 6. Bethe-Salpeter equation for two-particle Green function L(k, k'). T 1 (k, k') 
is the irreducible four point vertex with respect to the particle-hole channel Gk"+ q Gk" ■ 



Bethe-Salpeter equation: 
L(k,k';q) = 



— Gk+ q GkO~k,k'/T — G k+q G k T(k, k ; q)Gk'+ q Gk i 

— Gk+qGkO~k,k'/T 

-TJ2 G k+q G k T\k, k"; q)G k „ +q G k „L(k", k'; q), 



k" 
where k = (k, e n ) and T(k, k'\ q) is the full four-point vertex. r 7 (/c, k'\ q) is the irreducible 
four-point vertex, which is given by the Fourier transformation of the Ward identity in 
real space; V 1 = 5Y>/5G. The Bethe-Salpeter equation is expressed by Fig. 6. In 
later sections, we show that various linear transport coefficients are described in terms 
of L. (e.g., see. eq. (51).) In the conserving approximation, transport coefficient 
obtained by L(k,k') in eq. (18) automatically satisfies conservation laws [139, 140]. 
This is the reason why we refer to it as the conserving approximation. This is a great 
advantage of the FLEX approximation for the study of transport phenomena. In the 
FLEX approximation, irreducible four-point vertex will be given in eq. (89) in §4.1. 

2.3. Hot/cold-spot structure and T -linear resistivity in nearly AF metals 

As we have explained, important advance in HTSCs has been achieved by using the 
Fermi liquid theory with strong AF fluctuations [54, 56, 57, 52, 25, 24]. One of the 
most important predictions given by these theories is the "hot/cold-spot structure" of 
the quasiparticle damping rate, 7k = ImSk(— iS) (rk = l/27k). That is, 7k becomes 
anisotropic in the presence of AF fluctuations. This fact is very important to understand 
the transport phenomena in HTSCs. The portions of the Fermi surface at which 7k 
takes the maximum and minimum values are referred to as hot spots and [cold spots), 
respectively [25, 24]. According to the spin fluctuation theory, the hot spots usually 
exist around the crossing points with the AF Brillouin zone (AFBZ)-boundary, whereas 
the cold spot is at the points where the distance from the AFBZ-boundary is the largest. 
Their positions are shown in Fig. 4 (i). The electronic states around the cold spots play 
the major role for various transport phenomena. Note that the hot spot in YBCO shifts 
to the Brillouin zone boundary (point B), by reflecting the large DOS at (vr,0). 

Figure 7 shows the k-dependence of 7k on the Fermi surface given by the FLEX 
approximation [33]. In YBCO, the hot spot shifts from point C in Fig. 4 (i) to the 
BZ-boundary [point B], by reflecting the large DOS at the van- Hove singularity point 
(7r,0). Therefore, the spectral weight at the Fermi energy is strongly reduced around 
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Figure 7. The k-depcndcnce of 7k on the Fermi surface at various temperatures. The 
cold spot in YBCO (NCCO) is point A (B). 



(ir, 0) due to large 7hot 5 which is consistent with the ARPES experiments [86]. In 
NCCO, in contrast, the hot spot and the cold spot are located at point C and point 
B, respectively. We emphasize that the location of the cold spot in NCCO was first 
predicted by the FLEX approximation in 1999 [33], and it was later confirmed by ARPES 
[70]. According to the spin fluctuation theories [25], 7hot oc y/T, and 7 co id oc T except 
that £af ^> (A/^c)" 1 ; A/c c is shown in Fig. 4 (i). 

The critical value of U for a spin density wave (SDW) transition in the RPA (i.e., 
the mean-field approximation) is £/* TA ~ 2.3 in LSCO (n = 0.9). In YBCO and NCCO, 
U^ TA takes a much larger value, U* TA ~ 3.5 for both YBCO (n = 0.9) and NCCO 
(n = 1.1), since the nesting is not good due to the large next-nearest and third-nearest 
hopping integrals (t\ and £2)- Here, we consider to put the same U for both YBCO and 
NCCO. In YBCO, the dimensionless coupling constant UN(0) is large, since the saddle 
point is close to the Fermi level. Here, N(0) is the DOS at the Fermi level. Therefore, 
ImEk(O) takes a large value, which significantly reduces the interacting DOS and x q (0) 
at q = (-7T, 7r) [33]. In fact, the obtained ./Vflex(O) /N u=0 (0) is considerably smaller than 
unity in YBCO. On the other hand, the reduction in x q (0) due to ImEk(O) is small 
in NCCO since the saddle point is far below the Fermi level. For this reason, £af for 
NCCO is much larger than that for YBCO in the FLEX approximation. This result is 
consistent with experiments. 

The real-frequency Green function Gk( e ) is given by the analytic continuation of 



Gk(e n ) in eq. (13). In a Fermi liquid, the advanced Green function G( 



G 



k e 



iS) 



in the vicinity of e *- 

G£(e) 
£ k = 

7k = 



and |k| ~ kp can be represented as 

= V(e-£k-i7 k ), 
e k + S k (0)-^, E^ = z k E k , 



(19) 

(20) 

7k = ^k7k, (21) 

where z^ is the renormalization factor given by z^ = 1/(1 — <9Re£k(e)/<9e) e= o, an d -^k ^ s 
the renormalized quasi-particle spectrum, which is the solution of ReGj^ 1 (££) = 0. The 



ImE^O), 
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quasiparticle weight is given by 

Pk (e) = hmG£{e). (22) 

7T 

The DOS is expressed as N(e) = J2 k Pk(t)- I n the case of Zk7k <C /i, 

Pk(e) = ^k5(e - El) (23) 

for e f» 0. 

In the FLEX approximation, 7k is obtained by the analytic continuation of eq. 

(14): 

^ = \llj de [ cth ^ - th ^] ImV ^ e + ^)^n-q(e) (24) 



q 



where Vq(a> + i#) is given by the analytic continuation of eq. (15). In the spin 
fluctuation model (Millis-Monien-Pines model in eq. (8)), ImVq(O) in eq. (24) is replaced 
with (3[/ 2 /2)im Xq (0) = (3U 2 /2)u X q^/^1 + oo 2 ), where cu q = uj si + cu sf £ 2 F (q - Q) 2 . 
According to Refs. [25, 33], eq. (24) is approximately transformed to 

3U 2 f dk\ (nT) 2 

7k « -— / XQ^si-, - ( ■ — Tfrny; (25) 

4vr J FS Vk' 4w k -k'(^k-k' + 7rT/2) 

According to eq. (25), 7hot (7k at the hot spot) in 2D systems is given by [25]: 

7 hot oc T£ AF for ttT/2uj s{ > 1 (26) 

7hot ocT 2 £ AF for irT/2u s{ < 1. (27) 

Since £^ F oc T" 1 , 7 hot oc \/T for any value of cj si /T: This result is recognized in the 
numerical study in Fig. 7. Also, 7 co id in 2D systems is obtained as [25]: 

7cold oc T for 7iT/2u si ~ (£afAA; c ) 2 , (28) 

7coidOcT 2 for 7rT/2^ sf < (a F AA; c ) 2 . (29) 

According to eqs. (9) and (10), nT/2u s f is constant if 6 ~ 0, and it is of the order 
of 0(1) in optimally-doped HTSCs. Therefore, 7 co id oc T when £afA/c c ~ 0(1), and 
7coid oc T 2 when £afA/c c ^> 1. [Here, Ak c represents the "deviation from the nesting 
condition at the cold spot"; see Fig. 4.] We comment that the almost all part of the 
Fermi surface becomes the cold spot in the case of £afAA; c ~^> 1. 

Now, we discuss the temperature dependence of resistivity according to the SCR 
theory [54], by dropping the CVC. In the SCR theory, the resistivity is derived from 
the Born approximation p oc (7k)FS = Sk7kPk(0). In the case of u s f ^> T, eq. (25) is 
simplified as 

^ « E ^ImV q (0) Pk+q (0), (30) 

q 
where Vq(0) = dVqiuj + i5)/du\ LLl=0 . It can also be derived directly from eq. (24) by 
using the relation fj^ de[cth(e/2T) - th(e/2T)]e = (ttT) 2 . Using eqs. (8) and (30), the 
temperature dependence of the resistivity is given by [54, 147] 

Pscr oc ( 7k )FS oc T 2 ^, (31) 
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where d is the dimension of the system. In deriving eq. (31), we utilized the fact that 
the q-dependence of ImXq(O) = (n/2) X]kPk(0)Pk+ q (0) is moderate. Since £af oc T -0 - 5 
near the AF-QCP [54] , pscr is proportional to T d l 2 [d = 2, 3). We stress that the Fermi 
liquid behavior Pscr oc T 2 is recovered when £af =constant away from the AF-QCP, 
like in over-doped systems at low temperatures. 

In §4.2, we will calculate the resistivity using the FLEX approximation (and 
FLEX+T-matrix approximation), based on the linear response theory. In Fig. 18, 
"FLEX" represents the resistivity obtained by dropping the CVC, and "FLEX+CVC" 
represents p given by the FLEX+CVC approximation; the latter gives the correct result. 
Consistently with eq. (31), the obtained p follows an approximate T-linear behavior in 
under-doped LSCO and NCCO, and it shows a T 2 -like behavior in the over-doped 
NCCO. As for the resistivity, the CVC is quantitatively important. In later sections, 
we explain that the CVC completely changes the T-dependence of R-u, Ap/p and v. 

However, the Born approximation (eq. (31)) gives overestimated values when 7^ is 
highly anisotropic, as pointed out by Refs. [24, 148]: According to the linear response 
theory, the correct resistivity is give by 

P oc l/(7k lw JL)FS ~ 7coid- (32) 

According to eqs. (28) and (29), p shows the T-linear behavior in the case of 
£afA/c c ~ 0(1), whereas p oc T 2 in the case of ^\pAk c ^ 1- 

In optimally doped or slightly under-doped YBCO and LSCO, p shows an 
approximate T-linear behavior, which means that the relation £afAA; c ~ 0(1) is satisfied 
in these compounds above T c . In fact, if £afA/c c ^> 1, the dominant part of the Fermi 
surface should be the cold spot. However, this result is inconsistent with the ARPES 
measurements [86]. In electron-doped systems, on the other hand, £afAA; c ^> 1 seems 
to be realized at low temperatures even in optimally-doped systems, since £af reaches 
100 just above T c in NCCO [122]. The relationship ^ApAk c ~^> 1 is also satisfied in the 
present numerical study. In this case, 7k becomes highly anisotropic and p ~ 7 co id ~ T 2 
is realized [25]. 

According to eq. (32), the relationship p oc T 2 holds in the close vicinity of the AF- 
QCP; £afA/c c ~^> 1. However, Rosch [148] pointed out that the relationship p = a + bT 
(a, b > 0) holds even at the AF-QCP when the quasiparticle damping rate due to 
impurities 7i mp is finite: In fact, at sufficiently low temperatures where elastic scattering 
is dominant (7i mp ~^> 7k), then eq. (32) becomes 

p OC l/((7k + 7imp) _1 )FS OC 7 imp + (7k)FS- (33) 

Therefore, according to eq. (31), p — p oc T d l 2 (d = 2,3) holds near the AF-QCP in 
the presence of impurities. 

Here, we explained the T-linear resistivity in HTSCs based on the spin fluctuation 
theory. There are other theories which accounts for the T-linear resistivity. As discussed 
in §1.4, the marginal Fermi liquid hypothesis was proposed by Varma et al. [118, 119]. 
Also, the Tomonaga-Luttinger model with two types of relaxation times was proposed 
by Anderson [23]. Unfortunately, based on these theories, it is difficult to calculate other 
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transport coefficients such as i?H and Ap/p . In later sections, we will explain various 
anomalous transport phenomena in HTSCs in a unified way based on the Fermi liquid 
theory. 

3. Anomalous transport phenomena in nearly AF Fermi liquids 

3.1. Results by RTA 

Before investigating the transport coefficients using the microscopic Fermi liquid theory, 
we briefly review the relaxation time approximation (RTA) based on the Bloch- 
Boltzmann theory [1]. The CVC provides correction for the RTA. In later sections, 
we will see that the CVC becomes crucial near the AF QCP. The Boltzmann equation 
in a nonequilibrium steady state is expressed as [1] 

()J]A -e(E + 4xH)-V k / k , (34) 



°^ / scatt. 

where we put c = K — 1, and — e (e > 0) is the charge of electron. / k is the 
distribution function in a nonequilibrium steady state, and E and H are the electric 
and magnetic fields, respectively. (<9/ k /<9£) scat t. represents the rate of change in / k due 
to scattering between quasiparticles, which is called the collision integral. Using the 
scattering amplitude /(k, k; q) for (k, k') <-> (k + q, k' — q), it is given by [1] 

(§r) = -£ 7 ( k ' k ;<o \f*f*o- - J"-*)* 1 - /*-*) 

\ / scatt. k/q 

-(l-/ k )(l-/ k ,)/ k+q /k'-q] (35) 

where the first (second) term represents the outgoing (incoming) scattering process. To 
derive the conductivity, we have to linearize the Boltzmann equations (34) and (35) 
with respect to E. 

In solving these equations, we frequently apply the RTA; (<9/ k /(9t) scatt . = — g k /T k , 
where g k = / k — /£ and /£ = (e^ k ~^/ fcBT + 1) _1 is the equilibrium distribution function. 
Then, the linearized Boltzmann equation is simplified as 

& = _ eE ■ v (-^) + e (4 x H) • V k £ k . (36) 



r k \ <9 £ k 

The solution of eq. (36) is given by 

<?k = - (l - er k (v x H) • V k ) "' eE • r k v k ("|^) • (37) 

Although it is a crude approximation, the RTA can successfully explain the various 
transport phenomena in metals with weak correlation. 

In the RTA, the conductivity is given by cr^J A = — 2e J2 v kM g k / 'E v , where the factor 
2 is attributed to the spin degeneracy. When the k-dependence of r k is moderate, the 
functional form of a^ under the magnetic field is a^ u = tF^ u (tH). Here, we assume 
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H || z. Then, (vx H) • Vk = —H z (v k x Vk) z - In the RTA, the longitudinal conductivity, 
Hall conductivity, and magnetoconductivity are given by 

*T A =2e 2 ^(-£^ k ,-r k ^, (38) 

a^ A = - 2e 3 H z ^ (-^- J ^(r k 4 X V k ) z r k v ky , (39) 

Aa*J A = 2e A H 2 z £ (-$£) ^(r k 4 x V k )^ k „ (40) 

where (rki>k x Vk) 2 = Tk(v x d y — v y d x ). When the k-dependence of Tk is moderate, 

a xx oc r, (7^ oc t 2 # z , A^J A oc t 3 H 2 z , (41) 

which is called the Kohler's rule. 

Based on the RTA, Hussey et al. studied transport coefficients in heavily over- 
doped T12201 [149] and LSCO [150]. They determined the anisotropy of r k in T12201 
by measuring the Angle-dependent Magnetoresistance Oscillation (AMRO) along c-axis 
[151], and calculated transport coefficients using eqs. (38)- (40). The derived i?§ TA agrees 
with the experimental Hall coefficient R^ p in heavily over-doped samples (T c m 10 K), 
where the Hall coefficient is small and its temperature dependence is tiny. As the doping 
decreases, R^ p increases quickly and its temperature dependence becomes prominent, 
whereas the doping dependence of i?y TA is more moderate. 

In optimally doped HTSCs, however, the RTA does not work: Stojkovic and Pines 
[25] attempted to explain the violation of Kohler's rule in HTSCs based on the RTA. 
They assumed the highly anisotropic Tk model (hot /cold-spot model), where only the 
quasiparticles near the cold spot contribute toward the transport phenomena. In this 
model, i?H can take a large value since the "effective carrier density for the transport 
phenomena" is reduced. They calculated Tk based on the Millis-Monien-Pines model; 
the anisotropy of Tk reaches 100 in the optimally doped YBCO, which is too large to be 
consistent with the ARPES measurements. The reason for this overestimation is that 
self-consistency is not imposed in their calculations. In spite of the large anisotropy of Tk, 
the obtained enhancement ratio of R^ is about two according to Ref. [27]. Therefore, 
this scenario cannot account for the large i?H in under-doped systems. It should be 
stressed that the magnetoresistance becomes 100 times greater than the experimental 
value when the anisotropy of Tk is of the order of 100 [26]. Therefore, the highly 
anisotropic r^ model is not applicable for optimally- doped HTSCs. In the next section, 
we explain the various anomalous transport phenomena in HTSCs "all together", by 
considering the CVC. 

3.2. Physical meaning of the CVC in nearly AF metals 

In the RTA where gk = — S'k/Tk is assumed, the deviation from the equilibrium 
distribution function g^ dissipates with time since g^ oc e _ */ Tk . This oversimplification 
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Figure 8. Decay process of a quasiparticle at k in HTSCs. (i) Due to strong AF 
fluctuations, this quasiparticle at k is scattered to k + q (q ~ Q), creating a particle- 
hole pair (k', k' — q). (ii) The current of the particle-hole excitation in (i) is almost 
canceled by another particle-hole excitation (k", k" — q) for k" = — k'. Therefore, the 
CVC is given by the quasiparticle at k + Q alone. The process (ii) is the Umklapp 
scattering. 



frequently leads to serious unphysical results. For example, in the absence of electric 
field, the RTA predicts that g^ always vanishes when t — > oo since r^ > at finite 
temperatures. However, g^ should remain finite when the Umklapp scattering process 
is absent, because of the momentum conservation laws [1, 42]. Therefore, the RTA 
violates the conservation laws. To satisfy the conservation laws, the incoming scattering 
from the other states to k has to be taken into account by solving eq. (35). This 
scattering process is represented by the CVC in the microscopic Fermi liquid theory. 
Here, we intuitively discuss the CVC in nearly AF metals, where the quasiparticles are 
scattered by strong AF fluctuations with q ~ Q = (jr, n). Then, the momentum transfer 
q in the scattering process is restricted to ~ Q, as shown in Fig. 8. According to the 
momentum conservation law, the quasiparticle at k' is scattered to k' — q as shown in 
Fig. 8 (i). The current due to the quasiparticle at k' — q and the hole at k' is given by 
£V-q — ?V- This current almost cancels after performing the k'-summation: In fact, the 
current of the particle- hole pair (k" — q, k") for k" = — k', which is shown in Fig. 8 (ii), 
is W-k'-q — W-k' ~ — #k'-q + ^k' since 2q « 2Q is a reciprocal lattice vector. Therefore, 
the CVC is given only by the quasiparticle at k + q; t>k+q- In this case, the conductivity 
does not diverge because of the existence of the Umklapp processes, e.g., the process in 
Fig. 8 (ii). 

In nearly AF metals, the total current Jk is approximately parallel to v^ + ^k+Q, 



Anomalous transport phenomena in Fermi liquids with strong magnetic fluctuations 23 

which implies that J^ is not perpendicular to the Fermi surface. This fact is the origin 
of the enhancement of /?h [33]. The anomalous k-dependence of J^ becomes much more 
prominent near the AF QCP due to the multiple backscattering of the quasiparticle 
between that at k and that at k + Q. [The schematic behavior of J^ is shown in 
Fig. 14.] Thus far, we discussed only the two-body scattering process, and ignored the 
higher-order processes. This simplicity is justified in good metals where 'y/Ep <C 1, 
which is satisfied in HTSCs except for heavily under-doped compounds [28]. 

3.3. Analysis of the CVC based on the Fermi liquid theory 

Here, we calculate the CVC based on the microscopic Fermi liquid theory. In principle, 
we can also calculate the CVC based on the Bloch-Boltzmann theory, by analyzing eq. 
(35) using the variational principle [1]. However, a systematic calculation of the CVC 
is very difficult, particularly in the presence of a magnetic field. Therefore, we analyze 
the CVC based on the linear response theory since we can utilize the powerful field 
theoretical techniques. In the linear response theory [152, 153, 154], the conductivity is 
given by 

<vM = ^ K» - KM , (42) 

where K^ u (u) is the retarded current-current correlation function. 

Here, we consider the conductivities in the presence of the uniform magnetic field. 
Then, the hopping integral between site % and j exhibits the Peierls phase: 

tm,j = t° m>j exp[-ie(A m + Aj) • (r m - r^/2], (43) 

where £° ■ is the original hopping integral, and A m is the vector potential at site m. 
Then, the Hamiltonian and the velocity operator under the magnetic field, Ha and j^ 
respectively, are given by 

h a = Yl Wi/K + YE"" 1 ^™!' ( 44 ) 

(m,j),a m 

tf(r m ) = 'Y^i[H A ,v m c\ no c ma \ (45) 

(7 

Here, we assume that the vector potential is given by Aj = Ae lcl ' Tj , and take the limit 
q — ► at the final stage of the calculation [44]. (The magnetic field is H = iqx A.) Then, 
eq. (43) can be expanded as t mJ = t° mj [l - i(e/2)A ■ (r m - r^e^^ + e** T *)\ + 0{A 2 ). 
After the Fourier transformation, both the Hamiltonian and the velocity operator of the 
order of 0(A) are given by [38] 

H A =H-A-j(-q)+0(A 2 ), (46) 

x,y 

tf (q ; ) = JM) ~ £ Aj^q! - q) + 0(A 2 ), (47) 

a 

where ^(q) = -eE^^^c^^c^,, and j Q/i (q) = e 2 ^ k 9 a/1 e^'c[_ q/vCk+qA(; . 
Here, <9 M = d/dk^, d^ a = d 2 /dk^dk a , and — e (e > 0) is the electron charge. Therefore, 
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Kp V (u) is given by the analytic continuation of the following function for uoi > 

[43, 44, 38]: 

^ r 1/T 

K iaf {w l )= J2 / dre-^(T T ^(m q ,0)^(0,r)) A (48) 

1/T 

dre'^ (7^(0,0)^(0,^) (49) 

o 

+ f] A a f f ' drdr'e-^ {-T ■ (T r ^(q, 0)j Qi ,(-q, r)> 

+ (T T j M (q, 0)^(0, r) j Q (-q, r')>} (50) 

+ 0(A 2 ), 

where c^z = 27rT7 is the Matsubara frequency; here we promise / represents an integer 
and n is a half-integer, respectively. 

Hereafter, we ignore the spin indices to simplify expressions. According to eq. (49), 
K xx (iui) without the magnetic field is given by [40] 

K xx (iu>i) = - 2e 2 T ^ v{J a .u£, a Z(kn, k ' m ; iuj i) 

n,m, k,k' 

= -2e 2 T"£vlgX, (51) 

n,k 

where g£ = G£G£ + '- G£ = Gk(e n ) is the Green function where n is a half-integer, and 
L{kn,k'm;iuJi) is the two-particle Green function in eq. (18). i>£ = de^/dk n is the 
velocity of the free electron, and 

A? = v£ + T J2 LOm, k'm; «*,)<£"< (52) 

k',m 

is the three-point vertex. 

In the same way, Kohno and Yamada [43] derived the if -linear term of K xy (iui) 
from eq. (50) as 

K xy {iu x ) = i ■ e 3 HT J2 E i d ^ +l ■ G l ~ G l +l ■ d ^1 

n,k fi,u 

+ [6 point VC term], (53) 

where e uu \ is an antisymmetric tensor with e xyz — 1. In deriving eq. (53), we used the 
relation H z = i(q x A y — q y A x ) and took the limit q — » at the final stage of calculation 
[44] . Several kinds of Ward identities have to be correctly applied to maintain the gauge 
invariance [43]. 

In order to derive the conductivity 0)^,(0;), we have to perform the analytic 
continuations of eqs. (51) and (53); lo\ (I > 0) — > u + iS. Then, the e n summation 
is replaced with the integrations along the three cut lines in the complex plane in 
Fig. 9, together with the thermal factor (47riT) _1 th(e/2T). In region 1, g£ becomes 
(?k (e; a;) = G^(e)G^(e + u). In the same way, it becomes g^ (e; a;) = G^(e)G^(e + uo) 



And «-» a n\l 



tfivz 
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Figure 9. The analytic regions 1-3 as a function of a complex variable e (or e'). Here, 



. n:l 



n;l 



we put Im u> > 0. From each region, g k ' and A kx are analytically continued to become 
g£'(e;u>) and J^l{t;w) (i=l-3), respectively. 



and <7k (e;a>) = G^(e)G^(e + c<j) in region 2 and 3, respectively. Here, we stress the 
relation g^(e; 0) = 7rp k (e)/7^ ps 7r5(i?k)/7k! see eqs. (19)-(21). As a result, lf xx (a;H-i5) 
is given by [40] 

f°° de 
K xx {uj + iS)= - 2e 2 



+ th 



— oo 

e + uj 
IT 



Aiti 



th— Kjg{e,u) 

th^)^ 2 J(e;c)-thi±^^)(e^) 



«2(c;w) 



k 



vLoi ■ (e;^)^2(e;w), 



(54) 
(55) 



rW 



»,:( 



where i — 1,2, 3. ./^(e; w ) ^ s gi ven by the analytic continuation of A£ x from region i. It 
is expressed as 



■>£('. «0 



k', j=l,2,3 *' -t 



47T2 



r k,k'( e > e '^)^k' V^K'^ 



(56) 



where ^,(e, e'; a;) is given by the analytic continuation of the full four-point vertex V 
in eq. (18) from region (i,j) for (e, e'). Thermal factors [such as (47ri) _1 th(a;/2T) and 
(47ri) _1 ch(x/2T)] are included in the definition of T^ 3 V . The explicit expression for 7^'£, 
is given in eq. (12) of Ref. [40]. 

The DC-conductivity a xx is given by dhn.K xx /du\ UJ= Q. For example, if we take the 
^-derivative of the thermal factor (th^j^ — th^) associated with g^ in eq. (54), we 
obtain 



e 2 Myz k (-df/de) Ei vlJ^(0;0)/^, 



Z^ 



(57) 



where /°(e) = (e e / T + l)" 1 ; note that z k (-df°/de) E ^ = (-df°/de) Ek = 5(E k ) at 
sufficiently low temperatures. Since J^ x (e; 0) is not singular with respect to 7 _1 , eq. 
(57) is proportional to 7 _1 , which diverges when 7 — > 0. On the other hand, if we take 
the cu-derivative of the thermal factor in front of g( 1,3 ) in eq. (54), the obtained term is 
(7 ). To derive the exact expression for o~ xx with respect to 0(7 _1 ), we also have to take 



Anomalous transport phenomena in Fermi liquids with strong magnetic fluctuations 26 

the cu-derivative of the thermal factor in front of g^ 2 ' in K xx (e; a;). As a result, the exact 
expression for a xx is given by eq. (57) by replacing v^ x with v\^ x = i> ka . + <9ReSk(0)/(9A; :r 

[40]. 

(2) 

J^J is called the total current since it contains the CVC discussed in previous 
sections. Hereafter, we denote J k;r ;(e, 0) and T k k ,(e, e'; 0) as J^ x (e) and T k y(e,e'). The 
analytic continuation of K xy (iuJi) had been performed in Refs. [43, 44]. To summarize, 
the general expressions for a xx and o~ xy , which are exact within the most divergent term 
with respect to 0(7 _1 ), are given by 



0~xy 



df°\ A 

e3 TT ^ ( 9f°\ J k _ ^ ^ A 



o- xx = e 2 ^2 ^k ( — 7T- ) vk x — , (58) 



- e i H H z -{- d -f) -(^ x ^- ( 59 ) 

2 k v 9e ^ 7k 7k 

Equation (58) was derived by Eliashberg [40]. Equation (59) was derived by 
Fukuyama et al within the Born approximation [44], and it was proved to be correct in 
Fermi liquids [43]. Note that we dropped the [6 point VC term] in eq. (53) since they 
are less singular with respect to 7^ 1 [43]. In particular, its contribution to a xy vanishes 
in the FLEX approximation. In the same way, the present author has derived the exact 
expressions for the magneto-conductivity Aa xx = a xx (H z ) — o~ xx (0) [38] and the Peltier 
coefficient a xy = E y /{—V x T) in the presence of H z [39]. 

Apparently, eqs. (58) and (59) become equal to the RTA results given in eqs. (38) 
and (39) if we replace the total current J^ with the quasiparticle velocity v^, that is, if 
we drop the CVC. The RTA had been frequently used in analyzing transport anomaly 
in high-T c cuprates. However, we will explain that the neglect of the CVC frequently 
causes various unphysical results. 

According to Eliashberg, the total current J^ x can be rewritten as [40] 

/OO J / 
t-O, Ofl£V;o)./k<*(o), (60) 

K -oo 47n 

where v^ x = v^ x + <9Re£k(0)/<9/c x , and we put e = in Jk x (e) for simplicity. 7^J,(0, e') 
is the "irreducible" vertex with respect to g( l \ which is given by 

O ' e ') = « ( cth ^ - th ^) 2Imr k,k'(0, e' - iS). (61) 

Equation (60) is expressed in Fig. 10 (i). The CVC in the microscopic Fermi liquid 
theory is also called the backflow in the phenomenological Fermi liquid theory. 

Here, we analyze the Bethe-Salpeter equation (60): The solution of eq. (60) is real 
since 7^ ^,(0, e') in eq. (61) is pure imaginary. Since Imr k k ,(0, 0) = 0, we have to extract 
the e'-linear term of Imr k k ,(0, e'), which is derived from the cut of the particle- hole pair 
or that of the particle-particle pair in r kk ,(0,e') [31]. By this procedure, T 1 is divided 
into two T's. Therefore, e'-linear term of lmr(0, e') is given by [42] 

Imr kk , (0, e' - i8) = Tie J^ ^ ( k > k'; k' - q, k + q) 
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Figure 10. (i) Bethe-Salpeter equation (eq. (60)) for the total current Jk- (ii) 
Expression for 7^^/(0, e') in the microscopic Fermi liquid theory. The factor —1/2 is 
necessary to avoid double counting. 



x (^ k+q (0)p k ' +q (0) - -p k+q (0)p k ,_ q (0) 1 , (62) 

where T(k, k'; k' — q, k + q) is the full four-point vertex at the Fermi level, which is 
a real function of (k, k', q). The factor | in front of the second term in the curly 
bracket in eq. (62) is necessary to avoid double counting. Then, 7^ k ,(0,e') is given by 
i[cth(e/2T) — th(e/2T)]2ImT(0, e' — iS); it is schematically shown in Fig. 10 (ii). 

After changing momentum variables, the Bethe-Salpeter equation (60) is 
transformed into [42] 



AJ k 



t7 k + AJ k , 

-(0a) P k +q(0) 

27lc+q 



E^-(Oa) , 
^k.lc+q" 



J, 



k+q 



E^-(Ob) Pk'-q(O) 
J k,k' 



-q 



k'q 



k'q 



2 7l 



Jk' 



k'-q 



k'-q 



E^-(0c) Pk'(0) 
k , k ' k ' 2 7k ' 

k 'q 



Jk', 



(63) 



(64) 



where we used the relation j™ oo de[cth{e/2T) - th(e/2T)]e = (ttT) 2 . T^ {a = a,b,c) 
are functions of (k, k', q), and they represent the forward scattering amplitude [31, 40]. 
The expressions for T^ 0a ^ (a = a, b, c) at sufficiently low temperatures are given by [42] 



TKOa) 
J k,k+q 

rfiOb) 
^k,k'-q 

7-(0c) 
J k,k' 



^^ir 2 (k, k'; k' - q, k + q)p k '(0)p k '-q(0), 
^p-h 2 (k, k'; k' - q, k + q)pk'(0)p k+q (0), 
- ^2!ir 2 (k, k'; k' - q, k + q)p k+ q(0)pk'-q(0). 



(65) 
(66) 
(67) 



Since the CVC due to T^ represents the hole current, the minus sign appears in eq. 
(67). In eqs. (65)- (67), we dropped spin indices in F to simplify equations: If spin indices 
are taken into account, \Y 2 in each equation is replaced with \Y 2 *^ ^ + rf j_ ; j_ -f + Hl-LiTd 
as explained in Ref. [42]. Equations (65)-(67) are expressed in Fig. 11 (i). Note that 
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Figure 11. (i) Diagrammatic representation of eqs. (65)-(67). Each hatched part 
represents the real part of the full antisymmetrized four point vertex T(k, k';k + 
q, k' — q). Each line with arrow represents the imaginary part of the Green function. 
Note that (a) becomes different from (b) if we violate the antisymmetric nature of 
r(k, k';k + q,k' — q) in the course of approximation, (ii) Diagrammatic representation 
ofeq. (68). 



eqs. (65) and (66) are equivalent since r(k, k'; k' — q, k + q) is full antisymmetrized as 
a consequence of the Pauli principle. 

In the same way, ImEk(— iS) = 7k at sufficiently low temperatures can be expressed 
as [31, 42] 

de' 



7k 



k' 

jE^V +q (o), 



^kW)-w(e') 



(6* 



k'q 



which is proportional to T 2 at the zero-temperature limit, if the dimension is slightly 
higher than two. [In pure 2D systems, 7k oc — T 2 lnT.] Equation (68) is shown in Fig. 11 

(ii). We note that 7 k is also given by 7 k = \ E k ' q ^k,k'- q Pk'-q( ) = ~\ E k ' q ^JW( )- 
Here, we calculate the CVC in a free dispersion model in the absence of Umklapp 
scattering according to Yamada and Yosida [42]. In this case, both 7k and Jk on the 
Fermi surface are isotropic, that is, 7k = 7& F and Jk = J ■ k/fep. By noticing the 
relationship Jk+q + Jk'-q — J\J — Jk and using eqs. (65)-(68), it is easy to verify that 
the CVC in eq. (64) is exactly given by 

AJ k = J k , (69) 

— * — * 

in a free dispersion model. If we put AJk = cJk (c is a constant), the solution of eq. (63) 
is given by Jk = Vk/(1— c), which diverges when c — > 1. Therefore, the conductivity given 
by eq. (58) diverges, which is a natural consequence of the momentum conservation of 
the system [42] . Yamada and Yosida also showed that Jk remains finite in the presence 
of the Umklapp processes [42]. 

Next, we dicuss the case where e£ is anisotropic. Then, the formal solution of eq. 
(60) is given by 



j k = £(i - £)i£tk 



(70) 
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Figure 12. (i) Diagrammatic representation of eq. (71). The wavy lines represent 
the spin fluctuations, Xq(0). (ii) Diagrammatic representation of U in eq. (71). U is 
the bare Coulomb interaction in the Hubbard model and A(k + q, q) is the three-point 
vertex that is irreducible with respect to U. 



where Ck,k' = / ^^^(0, e') • pk'(0)/2 r yy s t. When e k is anisotropic, quasiparticle current 
is not conserved in the normal scattering process; v k + v k ' 7^ v k+q + v k /_ q . Nonetheless 
of the fact, one can show that det{l — C} = as a result of the momentum conservation 
[46, 155]. Therefore, J k given in eq. (70) and o~ xx diverge in any anisotropic model if 
the Umklapp processes are absent. [Exactly speaking, p xx = l/a xx is proportional to 
T 2N ~ 2 when the iV-particle Umklapp scattering processes are present [155].] 

We stress that the FLEX approximation can reproduce the divergence of the 
conductivity in the absence of the Umklapp scatterings, if one consider the CVC 
correctly. This result is assured by the fact that the Ward identity T 1 = 5Y,/5G is 
consistently satisfied in the FLEX approximation [139, 140]. Therefore, the FLEX 
approximation will be appropriate for the study of transport phenomena in correlated 
electron systems. 

Next, we consider the role of the CVC in nearly AF Fermi liquids in the presence 
of Umklapp scatterings, where the electron-electron scatterings are mainly given by the 
strong AF fluctuations with q ~ Q. In this case, we can approximate the full four-point 
vertex in eqs. (65)- (67) as 

U 2 

r s 1)S2 ;s 3 , S 4( k > k '' k ' ~ *1> k + <*) ~ "2" {Xq(°)^l.»4 ' ^3 

-Xk-k'+q(0Kl,«3 ' #82,84} , ( 71 ) 

which satisfies the antisymmetric nature of T that is the consequence of the Pauli 
principle. The diagrammatic representation of eq. (71) is shown in Fig. 12 (i). Here, 
Si (i — 1 ~ 4) represents the spin index and a is the Pauli matrix vector. U in eq. 
(71) is the effective interaction between the electrons and spin fluctuations, which is the 
renormalized Coulomb interaction due to the irreducible three-point vertex shown in 
Fig. 12 (ii). 

In the present subsection, however, we analyze the CVC using the following more 
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Figure 13. (i) According to the approximation in eq. (72), 7^ k °[ 7^ k /_ q , and 7^ £; 
in eqs. (65)-(67) are approximated as (A)-(C), respectively. (A) corresponds to the 
Maki-Thompson term in the FLEX approximation, which represents the CVC (MT) 
in Fig. 8. (B) and (C) correspond to the Aslamazov-Larkin terms, which represent the 
(AL1) and (AL2) in Fig. 8. (ii) The expression for 7k according to the approximation 
in eq. (72). 



simplified approximation for T: 

U 2 
r B1 , aa ;. SlW (k, k'; k' - q, k + q) « _ x ^(0)a Sl , S4 • a S2 , S3 . (72) 

Although eq. (72) violates the Pauli principle, this approximation produces the 
dominant CVC (Maki-Thompson term) near AF-QCP. In §3.4, we will analyze the CVC 
using eq. (71), and show that the analysis based on eq. (72) is justified. According to 
the present approximation, T(° a > (a = a, b, c) in eqs. (65)-(67) are expressed by (A)-(C) 
in Fig. 13 (i), respectively, where wavy lines represent the spin fluctuations. The CVCs 
expressed in (A) and (B) correspond to the current due to the quasiparticles at k + q 
and k' — q, respectively, and the CVC in (C) corresponds to the current due to the 
quasi-hole at k' in Fig. 8. In the field theory, (A) is called the Maki-Thompson term, 
and (B) and (C) are called the Aslamazov-Larkin terms. As explained in Fig. 8, the 
Aslamazov-Larkin terms approximately disappears [33] whereas the Maki-Thompson 
term plays an important role when £af 3> 1 and Q ~ (re, ir). 

According to eq. (72), \Y 2 in T^ becomes (3U 4 /A)[x s q (0)] 2 after taking the 
summation of spin indices. However, we replace |T 2 with (3U 4 /2)[x s „{0)} 2 since the 
term (A) in Fig. 13 appears twice if we use eq. (71). Then, the CVC in eq. (64) and 
7k given in eq. (68) near the AF QCP are obtained as 

AJ k = ^ T ^ix s ^ 2 l^(0) P -^J^ (73) 

q Z 7k+q 

Q/T4 

7 k = E( 7rT ) 2 ^-^(°)] 2lm ^( )^ q (°)' ( 74 ) 

q 

where we used the relation lm%q(0) = (tt/2) ^k/°k(0)pk+q(0). Since lm%g(0) = 
ImXq(0)[£7x q (C))] 2 , the Bethe-Salpeter equation (73) and 7k in eq. (74) in a nearly 
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AF Fermi liquid are given by 

& = * + £ ^(^) 2 lm^(0)^±^ J k+q . (75) 

q 4 7k+q 






3C/ 



2 



7k = ^ — (7rT) 2 Imx q (0) Pk+q (0). (76) 

q 
Note that eq. (76) is equivalent to 7k in the spin fluctuation theory in eq. (30). 

In this subsection, we have neglected the energy-dependence of the four-point 
vertices [in eqs. (65)- (67)] and that of the spin susceptibilities [in eqs. (73)- (76)] to 
simplify the discussion. For this reason, eqs. (75) and (76) are appropriate only for 
^sf ^ T where ^-dependence of Xq( u ) can be ignored. In the case of o; s f <C T, ijcT) 2 in 
eqs. (75) and (76) should be replaced with P^ sf [cth(e/2T) — th(e/2T)]ede = 4Tw s f. For 
any value of ou s f/T, they are expressed by eqs. (36) and (33) in Ref. [33], respectively. 

Here, we approximately solve eq. (75) in the case of £af ~^> 1 an d ^ F Ak c ~ 0(1). 
Here, we assume k is close to point A in Fig. 4 (i). Because x q (0) takes a large 
value only when |q — Q| < £^p, the CVC term in eq. (75) can be expressed as 
AJ k w (J k */ 7 k*)E q (3£/74)(7rT) 2 Imx q (0), where (fc*^*) = (-k y ,-k x ) for k x k y > 0, 

and (k*,k*) = (k y ,k x ) for k x k y < as shown in Fig. 4 (i). Considering eq. (76), 
A J k sa J k » in the case of £af ^ 1- The more detailed expression for the CVC term A J k 
is given by [33] 

AJ k S3 ( J k / )|k'-k*|<l/£ A F 

ss J k * • ( cos(0 k ,-0 k *) )|k'-k*|<i/? AF « a k J k *, (77) 

where a k ~ (1 — c/£ AF ) < 1 and c ~ 0(1) is a constant. The k-dependence of a k will 
be moderate if £afA/c c > 1. By this simplification, eq. (75) becomes 

J k = 4 + a k J k *, (78) 

and the solution is given by [33] 

Jk = Z 2"(Vk + «k^k*)- ( 79 ) 

l-a k 

Here, we examine the k-dependence of J k given in eq. (79): (a) near points A and B, 
J k is parallel to -u k by symmetry. At point A, J k = Vk/(1 + a k ) ~ \vy_ since -u k = — ik*. 
(b) Near point C, J k ~ (^ F /2c)(-u k + w k *), which is nearly parallel to the AFBZ. These 
behaviors of J k together with v k are schematically shown in Fig. 14 (i). In the present 
discussion, we assumed that £afAA; c ~ 0(1), where Ak c represents the deviation from 
the nesting condition at the cold spot in Fig. 4 (i). This condition seems to be satisfied 
in slightly under-doped YBCO and LSCO above T c , as discussed in §2.3. In Fig. 4 (ii), 
we show an "effective Fermi surface" obtained by bending the true Fermi surface such 
that it is orthogonal to J k around the cold spot [20] . we discuss the role of the CVC in 
the Hall coefficient based on the concept of the effective Fermi surface in §4.2. 

When £afAA; c ^> 1, one may think that the CVC is small around the cold spot, and 
therefore the increment in Ru due to the CVC is also small. However, this expectation 
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Figure 14. (i) Schematic behavior of the total current Jk and the quasiparticle 
velocity v^ on the Fermi surface. The cold spot is located around the point A (B) 
in hole-doped (electron-doped) HTSCs. (ii) Effective Fermi surface for <r xy . The 
large curvature of the effective Fermi surface around the cold spot is the origin of the 
enhancement (and the sign-change) of i?n- 



is not true unless Xq( w ) is a step function like 9(1 — |Q — q|^AF)- I n fact, we verified 
that the CVC around the cold spot is significant even if £afA/c c is much larger than 
unity, by assuming a phenomeno logical Xqi 00 ) m odel given in eq. (8) [34]. In later 
sections, we demonstrate that i?n is prominently enhanced by the CVC using the FLEX 
approximation, which yields a realistic functional form of xtt( u )- We will show that the 
CVC is important in NCCO, although £afAA; c ^> 1 is realized at low temperatures. 

Finally, we comment that the Aslamazov-Larkin terms ((B) and (C) in Fig. 13) 
become important when the AF fluctuations are weak. In the FLEX approximation, 
one can reproduce the divergence of the conductivity in the absence of the Umklapp 
scatterings, only if all the CVCs ((A)-(C) in Fig. 13) are taken into consideration. 



3.4- CVC beyond one loop approximation 

In the previous subsection, we analyzed the CVC by assuming eq. (72), which 
corresponds to the one-loop approximation for the self-energy like the FLEX 
approximation. In this approximation, eq. (73) and eq. (74) contain the same Kernel 
function [Xq(0)] 2 ImXg(0)pk+q(0)] as a consequence of the Ward identity between the 
CVC and the imaginary part of the self-energy. For this reason, the factor Ok in eq. 
(79) approaches unity near AF-QCP, which assures the anomalous k-dependence of Jk 
near AF-QCP that shown in Fig. 14. However, one may be afraid that this result is 
an artifact due to the violation of the Pauli principle. In this subsection, we analyze 
the self-energy and the CVC beyond the one-loop approximation by assuming eq. (71). 
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Figure 15. (i) According to the approximation in eq. (71), 7^ k °[ , 7^ k'_qj anc ^ ^k k' 
in eqs. (65)-(67) are approximated as (A)-(F), respectively. (A), (B), (D), and (E) 
come from T( 0a ) and T^ b \ (C) and (F) come from T^ 0c \ (ii) The expression for 7 k 
according to the approximation in eq. (71). 



Here, we confirm that the CVC near the AF QCP is important even if the higher-loop 
diagrams are taken into account. 

According to eq. (71), 7k in eq. (68) is approximately given by 

7k =7^+7 k 2) , (80) 

7k 1} = eq. (74), 

7r(7rT) 2 5£/ 4 



(2) 

7k 



£ 



-xyo)\-'.(o)p k+<1 (oK_ q .(oK +q _ q ,(o), 



q,q' 

'" /-. m .„..! .,(2) 



where 7 k (> 0) and 7 k (< 0) are shown in Fig. 15 (ii). The integrand in eq. 
(81) takes a large value when |q — Q|, |q' — Q| < £^ F , under the condition that 
£k+ q = ek-q' = ek+q-q' = /i. When Q = (tt, 7r) , 7 k can be significant at the hot 
spot since k + q^k — q'«k* and k + q — q' « k in modulo (2tc, 2tt) . 

In the same way, the CVC in eq. (64) is given by the six diagrams shown in Fig. 
15. As discussed in §3.3, terms (B) and (C) are negligible in nearly AF Fermi liquids. 
Other terms are given by 

(A) = eq. (73), 



;d,e) 



/ t\2 cry-4 
E^ ± ^Xq(0)Xq'(0)p k -q'(0)p k+ q-q'(0) 



q,q 



x Pk ± q(o)r 



(F) 



E 



q,q' 



7k+q 

2 cr74 



k+qj 



(82) 



n(nT) 2 5U 



Aq(0)Xq'(0)pk-q'(0)Pk + q(0) 



p k+q _ q /(0) 7 

x ^k+q-q' 

7k+q-q' 



(83) 
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Figure 16. Diagrammatic representation of T(g) in eq. (87). 

Here, X q (0)x q '(0) takes a large value only when |q — Q|, |q' — Q| < ^f- In (D)-(F), the 
set of momenta (k, k', q) is changed to (k, q, q') due to the transformation of variables. 
Applying the same approximation that was used in eq. (77), the total CVC is given by 

(1) (2) _ (2) 

AJ k = a ^- J k . + 2a' k ^- J k * - a£^_J k , (84) 

7k 7k 7k 

where (1— a k ), (1— CK k ) and (1— a k ) are positive, and they approach zero in proportion to 
£^f near the AF-QCP. Then, the solution of the Bethe-Salpeter equation J k = w k + AJ k 
is given by 

7k(7k + <7 k 2) ) 



Jv 



[7k - « k7k 1} - (2< - <) 7k 2) ][ 7k + « k7k 1} + (2< + <) 7k 2) ] 



.,(!) i <W.,( 2 ) 



7k + <7 k 
When ^af ^ 1, eq. (85) is approximately given by 

</k = rn r5T i^k + #k* } • (86) 

2[7k - OkT^ - (2< " <)7 k 2) ] 

Near the AF-QCP, 7 k -a k 7 k - (2a k -a k )7 k 2) oc 7 k ^AF since 7k = 7 k + 7 k • Therefore, 
J k in eq.(86) shows a similar behavior to that given in eq. (79). As a result, the 
significance of the CVC in nearly AF Fermi liquid is confirmed even if we assume eq. 
(71), which satisfy the Pauli principle. 

Finally, we discuss the CVC composed of the higher loop terms with respect to 
Xq(0), by considering the following approximation for the full four-point vertex: 

r si)S2;S3iS4 (k, k'; k' - q, k + q) « T(q)a SltS4 ■ a S2 , S3 

- f (k - k' + q)<r sliS3 -<7 S2;S4 , (87) 

by considering that T mainly depends on the momentum transfer in the electron- 
electron scattering. f(q) is given by Fig. 16. Using the full four-point vertex, the 
spin susceptibility is given by x q (0) = xj(0) + T 2 ^ kk , )f£ , G k (e)G k+q (e)r(k, k'; k' - 
q,k + q)G k ,(e , )G' k '- q (e') ) where Xq (0) = -r£ k>e G k (e)G k+q (e). According to the 
approximation in eq. (87), X q (0) is given by 

X q (0) + 2[ Xq (0)] 2 f(q) 

-2T 2 ]T G , k (e)G k+q (e)f(k-k')G' k , +q ( e / )G' k ,( e '). (88) 

k,k',e,e' 
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The third term of eq. (88) is expected to be smaller than the second term since the 
integrand takes large values only when k — k' ~ Q. Since the q-dependence of x q (0) 
is moderate, the present analysis suggests that the approximate relation x q (0) oc f (q) 
holds for q ~ Q. For this reason, the solution of J^ obtained in this subsection, eq. (85) 
or (86), will be valid even in higher-loop approximations. 

As a result, the anomalous k-dependence of J^ in Fig. 14, which is derived in the 
previous section within the one-loop approximation, is expected to be realized even in 
the close vicinity of the AF QCP where the higher-loop diagrams becomes important. 
In later sections, we present a numerical study of the transport phenomena based on 
the FLEX+CVC approximation. Although it is a one-loop approximation with respect 
to the spin fluctuation, the obtained numerical results will be qualitatively reliable even 
in the under-doped region, since the CVC due to the higher-loop diagrams do not alter 
the one-loop approximation results qualitatively. 

4. Transport phenomena in HTSCs above T* 

4-1. Total current J^ 

Here, we discuss the transport phenomena in HTSCs above the pseudo-gap temperature 
T* where the SC fluctuations are negligibly small. In the FLEX+CVC approximation, 
the self-energy is given by eq. (14), and the vertex function in the Bethe-Salpeter 
equation, T^,, is given by eq. (61). The total current Jk is obtained by solving the 
Bethe-Salpeter equation [33] . According to the Ward identity, the irreducible four-point 
vertex is given by T 1 = 5T,/5G. In the FLEX approximation, the self-energy is given 
by the convolution of G and V as in eq. (14). Then, T 1 contains one Maki-Thompson 
term that is given by taking derivative of G, and two Aslamazov-Larkin terms that are 
given by taking derivative of V, which is composed of infinite G"s. These three terms 
corresponds to (A)-(C) in Fig. 13 [33, 57]. 

Considering the relation <5x q /5Gk' = — (Gk'+q + GW-q), we can derive r£ k , = 
5Sk/5Gk' i n the FLEX approximation as 

r k k'( e n' e n') = H:-k'( e n ~ £ n ') 



tJ^u 2 



\{ux s M) + iy + -{u x c M-iy 



^ ~qv ., , 2 



1,1 

x G k '- q (e n / ~~ ^/)Ck- q (e„ — LUi) 



-tJ^u 2 



^XqM) + I) 2 + ^XqM - I) 2 



1,1 

x Gk'+ q (en' + ^)G k _ q (e n - wj). (89) 

The first term V^iuoi) corresponds to the Maki-Thompson vertex correction. The last 
two terms in eq. (89), both of which contain [x q ] 2 , are the Aslamazov-Larkin vertex 
corrections. Hereafter, we drop the Aslamazov-Larkin terms since they are negligible in 
nearly AF Fermi liquids (xh ~^> Xqi Xq)i as we have discussed in previous sections. 
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Figure 17. (i) 8^ = taxi -1 (Jk x / Jk y ) and (ii) mean free path with CVC obtained by 
eq. (90); |Zk| = |>/k/7k|- —d0^/dk\\ at the cold spot (A) increases as T decreases. |£k| 
attains its local minimum at the cold spot for T < 0.02, which induces increases in the 
Nernst coefficient and magnetoresistance. We also plot the mean free path without 
CVC; |f k | = |« k /7k|. 



After the analytic continuation of the Maki-Thompson term iu>i — > u + iS from 
the region e n , e n t < and e n + u>i, e n i + u>i > 0, the Bethe-Salpeter equation in the 
FLEX+CVC theory for u = is given by [33], 



J\ 



k e 



W 



E/£ 



cth 



e' -e 
IT 



th4 

2T 



x ImV q (e 



\S) ■ |G k+q (e + e' + id)\ 2 ■ J k+q (e + e' 



)'k 



+E 



^)! Iml > q( o)^(2) /k+qi 



7k^ 



(90) 



where v^e) 
J[cth(e/2T) 



Vk(e k + ReSk(e)). In deriving eq. (90), we used the relation 
th(e/2T)]ede = (vrT) 2 , and assumed the relation co> sf ^> T. If we 
approximate Vq ~ 3f/ 2 x q /2, eq. (90) is equivalent to eq. (75). [Note that U = U 
in the FLEX approximation.] The numerical solution of eq. (90) is shown in Fig. 17, 
whose schematic behavior is shown in Fig. 14 (i). This singular k-dependence of J k is 
realized because a k in eq. (77) approaches unity for £af — ► oo. As discussed in §3.4, 
this result is not specific to the FLEX approximation since the same vertex function 
Pick'] appears in both (60) and (68) in the microscopic Fermi liquid theory, which is 
one of the Ward identity. For this reason, singular k-dependence of Jk in the case of 
£af ^ 1, which is shown in Fig. 17, is not specific to the FLEX approximation, but is 
a universal behavior in Fermi liquids near the AF QCP. 

In hole-doped systems, 7k attains its minimum (maximum) at point A (point B) as 
shown in Fig. 7 [24, 25]. At low T, the local minimum of | J k | is located at points A and 



Anomalous transport phenomena in Fermi liquids with strong magnetic fluctuations 37 

B, and its maximum is located at point C: this is schematically shown in Fig. 14 (i). 
Further, Fig. 17 (ii) shows that the mean free path with CVC, |Lk| = |Jk/7k|> attains its 
local minimum at point A for T < 0.02, since | Jk| rapidly increases as k deviates from 
the point A. In fact, eq. (79) suggests that |J k | oc ^af!^ + ^k*|- Therefore, |Jk| takes 
a large value near the AF QCP except for point A and B. In later sections, we explain 
that the a local minimum structure of |Lk| at point A becomes prominent below T*, 
which causes the significant increases in the Nernst coefficient and magnetoresistance. 
In Ref. [33], we solved the Bethe-Salpeter equation by considering both Maki- 
Thompson term and the Aslamazov-Larkin terms, and calculated both p and Rr. As a 
result, we verified numerically that the Aslamazov-Larkin terms are negligible in nearly 
AF metals. 

4-2. Resistivity and Hall coefficient 

First, we analyze the conductivity a xx using eq. (58). According to the approximate 
expression for J^ given in eq. (79), 



i^kx^fccJcold 



v L 



(91) 

- cold 



1 + «k 

at the cold spot of YBCO [point A]. We have used the relationship k = — k* at the 
cold spot. Since 1/(1 + «k) ~ 1/2 when £af ~^> 1, then a xx ~ ct^J a /2. Therefore, the 
resistivity p is slightly increased by the CVC [33]. Since the effect of the CVC on p is 
not large, the RTA result for p can be qualitatively justified. 

Figure 18 (i) shows the numerical results of p for LSCO (n = 0.92 and U = 4.5) 
obtained from the FLEX approximation and the FLEX+T-matrix approximation. In 
the FLEX approximation without CVC, p shows an approximate T-linear behavior for 
T > 0.02 (~ 80 K). Because of the CVC, p is slightly enhanced for a wide range of 
temperatures. Moreover, p obtained by the FLEX+CVC approximation shows a tiny 
"kink" structure at ~ To, below which AF fluctuations grows prominently. This result 
is consistent with the experimental results [167]. The kink becomes more prominent 
in the FLEX+T-matrix approximation, which will be discussed in §5.2. In optimally- 
doped YBCO and LSCO, the resistivity at 300 K is 200 — 300 /ificm, and it increases 
to 400 — 600 yuffcm in slightly under-doped compounds (n ~ 0.1) [156]. In Fig. 18 
(i), p = 1 ~ 250 /ificm at T = 0.08 ~ 320 K: The resistivity increases with U, and 
p ~ 450 /idem for YBCO with U — 8 and n = 0.9 [33], which is consistent with 
experimental values. However, the FLEX+CVC method cannot reproduce the huge 
resistivity in heavily under-doped compounds, since the self-consistency condition in 
the FLEX approximation tends to suppress 7k too strongly; see §9.2. In under-doped 
compounds, one has to take account of the effect of residual disorder since the residual 
resistivity grows prominently. We will explain the reason in §7.2 based on the Fermi 
liquid theory with strong AF fluctuations. 

Figure 18 (ii) shows the p for NCCO (U = 5.5) obtained from the FLEX 
approximation (within the RTA) and the FLEX+CVC approximation. In the case 
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Figure 18. (i) p for LSCO obtained by the FLEX (without CVC) and FLEX+CVC 
approximations, p is slightly enhanced because of the CVC. We also plot p obtained by 
the FLEX+T-matrix (without CVC) and [FLEX+T-matrix] +CVC approximations. 
The FLEX+T-matrix approximation coincides with the FLEX approximation for 
T ^> T* , where SC fluctuations disappear. T = 0.2 corresponds to 800 K. (ii) p 
for NCCO obtained by the FLEX (without CVC) and FLEX+CVC approximations. 



of n — 1.20 (over-doped), p shows a T 2 -like behavior below T = 0.1 ~ 400 K, 
which is consistent with experiments. In the case of n — 1.10 (under-doped), p 
shows an approximate T-linear behavior when the CVC is included in the calculation. 
Interestingly, pflex+cvc < Pflex(rta) for n = 1.10 below T = 0.06, since |J k | > |u k | 
due to the CVC in NCCO at low temperatures. 

Next, we discuss the Hall coefficient. Using the Onsager's relation a xy = —a, 
the general expression for a xy in eq. (59) and o^y 1 ^ in eq. (39) in 2D systems ( 
rewritten as [33] 
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7 k 
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(92) 
(93) 



where L k = J k /7 k , 9£ = tan _1 (J kx /J k?/ ) and 9£ = tan _1 (i^ CB /vk V ). dk\\ is the momentum 
tangent to the Fermi surface, which is depicted in Fig. 14 (i). Note that |Z k | = | J k /7 k | 
represents the mean free path by considering the CVC. In deriving above equations, we 
used the relation |f k |<9/<9A;|| = (z x u k ) ■ V k = (v k x V k ) z in 2D systems. Using the 
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relation \ v^ + v%_ d/dk\\ = (v^ x Vk) 2 , o xy in a 3D system is given by 



,2 i „,2 



'- ,h '-tJ„^t' ki |Jt| l«rj'i' (94) 

where (iSk represents the Fermi surface element, and dfcii is the momentum tangent to 
the Fermi surface and parallel to the xy-plane. We stress that the /c-derivative of 7k 
does not enter into the expression of a xy , while it exists in expressions of Aa xx and a xy . 

In the expression of o",^ A , {—d9^/dk\\) represents the curvature of the Fermi surface 
[168, 169]; in both hole-doped and electron-doped HTSCs, (—d9^/dk\\) ~ l/kp (> 0) 
on the Fermi surface. On the other hand, (—dd^/dku) exhibits strong k-dependence 
in nearly AF metals, as shown in Fig. 14 (i); (—dO^/dku) is positive around point A, 
whereas it is negative around point B. Since the cold spot is located around point A in 
hole-doped systems [33], the present study predicts that Rn > in hole-doped systems 
when the AF fluctuations are strong. [Note that the charge of electron is — e (e > 0) in 
the present study] On the other hand, Ru < in electron-doped systems since the cold 
spot in electron- doped systems is B [33]. 

For an intuitive understanding of the CVC, we introduce an "effective Fermi 
surface" obtained by bending the true Fermi surface such that it is orthogonal to Jk 
around the cold spot [20]. The effective Fermi surface in HTSCs is shown in Fig. 14 
(ii). It can be seen that the curvature of the effective Fermi surface, which is equal to 
(— d0£/dk\\) by definition, takes a large positive (negative) value around the cold spot 
in hole-doped (electron-doped) systems. Therefore, Ru in hole-doped (electron-doped) 
systems exhibits a large positive (negative) value at low temperatures [33]. 

Figures 19 (i), (iii), and (iv) show the numerical results for i?n obtained from 
the FLEX approximation by including the CVC. In hole-doped systems (n < 1), i?n 
increases as the doping S = [1 — n\ decreases. On the other hand, i?H for electron- 
doped systems (n > 1) becomes negative below T = 0.09 ~400 K. These results are 
consistent with the experimental results. Previously, negative Ru in electron-doped 
systems had been considered to be strong evidence of the breakdown of the Fermi liquid 
state since the curvature of the Fermi surface in HTSCs is positive everywhere. However, 
it is naturally explained in terms of the Fermi liquid picture since the curvature of the 
effective Fermi surface becomes negative around point B, as shown in Fig. 14 (ii) [33]. 

As shown in Fig. 19 (i), (iii), and (iv), i?n given by the RTA (without CVC) 
decreases with decreasing T, since the curvature of the true interacting Fermi surface 
becomes smaller in the FLEX approximation at low temperatures [33]: The deformation 
of the Fermi surface is caused by the strong k-dependence of the self-energy, which 
becomes prominent when AF fluctuations are strong. A tiny increment in R^ TA in 
YBCO below T = 0.02 is caused by the strong anisotropy of 7k- In the present 
calculation, 7hot/7coid is approximately 3 at T = 0.02 

Figure 19 (ii) shows the /cy-dependence of Aa xx (k\\) = v^ ■ Lk/I^kl and Aa xy (k\\) = 
|-^k| 2 ( _ 90^/dk\\), where a xx = <f FS dk\\Aa xx (k\\) and a xy = § ¥ ^dk\\Aa xy {k\\). In both 
o xx and a xy , the quasiparticles around the cold spot gives the dominant contributions. 
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Figure 19. (i) Rn for LSCO obtained by the FLEX+CVC approximation. In contrast 
to its prominent T-dependence, i?§ TA given by the RTA (without CVC) shows only 
a weak T-dependence. Rn obtained using the FLEX+T-matrix approximation is also 
shown; it starts to decrease below T* ~ 0.04 (~ 160) K (pseudo-gap behavior), (ii) k\\- 
dependence of Acr xx (ku) and Aa xy (ku) at T = 0.02. Note that a^ v = J ps dk\\ A(jfi V {k\\). 
(iii) R R for YBCO obtained using the FLEX+CVC method, (iv) i? H for NCCO 
obtained using the FLEX+CVC method. 



We can see that a xx is slightly reduced by including the CVC. On the other hand, 
o~„, obtained from the FLEX+CVC theory is considerably larger than al 1 ™, since 



■xy 



xy 



\d9^/dk\\\ 3> \d9^/dk\\\ at the cold spot. Aa xy (k\\) is large only around the cold spot, 
and it becomes negative around the hot spot. 

Here, we discuss the temperature dependence of Ru in detail using the approximate 
expression for Jk in eq. (79). Since dak/dk\\ = at the cold spot [point A] because of 
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the symmetry, 

dfk _ 1 f dih, dv k * \ 

dk\\ ~ l-al \dk\\ k dkJ [ } 

at the cold spot. By noticing the relationships (vk* x dvk*/dk\\) z = — (vk x dvk/dk\\) z 
and (vk* x dvk/dk») z = — (vk x dvk*/dkn) z , we obtain that 

: <x ^] \u(~ deJ A 

dk\\ ) V dk\\ I ,, 

II / n II ' cold 

= 1 2^ h^ » ( 96 ) 

l-«k V dk\\ J cold 

which is proportional to (1 — «k) x oc £^f [33]. In fact, Fig. 17 (i) shows that eq. (96) 
increases as the temperature decreases. As a result, Ru behaves as [33] 

|i? H |oc£i F . (97) 

Therefore, both the sign and the T-dependence of Ru in hole-doped HTSC are 
successfully reproduced in the present approach. 

Finally, we discuss electron-doped systems. At the point B in Fig. 4 (i), 
|k B — k B ' — Q| is equal to |k B — k B » — Q| for Q = (it, —it). Then, the CVC at k B , which 
is given by the second term of eq. (75), is approximately proportional to Jb> + Jb" = 0. 
Therefore, a^ = in the simplified Bethe-Salpeter equation (78) and J^ ~ Vk at point 
B. Since ak rapidly increases if k deviates from point B, (—dd£/dk\\) attains a large 
negative value around the cold spot of NCCO. Therefore, the negative sign of i?H is 
realized by considering the CVC. 

Finally, we stress that we have also studied the CVC by using a widely used 
phenomenological AF fluctuation model given in eq. (8). Assuming a reasonable set of 
parameters, we find that i?H is prominently enhanced due to the CVC [34]. Thus, the 
enhancement of i? B i n nearly AF metals is not an artifact in the FLEX approximation, 
but a universal phenomena near the AF QCP. 

4-3. Magnetoresistance 

Next, we study the orbital magnetoresistance Ap/p in HTSCs by considering the CVC. 
According to the linear response theory [1], the magnetoresistance is give by 

Ap/p = - Aa xx /a° xx - (a xy /a° xx ) 2 , (98) 

where a xx denotes the conductivity without a magnetic field, and Aa xx = cr xx (H z ) — a xx 
is the magneto-conductivity, which is always negative. 

To derive the magnetoresistance, we have to calculate the magneto-conductivity, 
which is given in eq. (40) in the RTA. By performing the partial integration, eq. (40) 
for 2D system is rewritten as 

A RTA _ tt2 e J. Uft ll |„-? 

Aa xx - - n- — <p 7— r- \Vk 



■ ° r ra (2tt)= 




— , (99) 

7k 
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where l^ = i>k/7k- Here, we used the relation \v^\d/dk\\ — (z x v^) ■ Vk = (v^ x Vk) z 
in 2D systems. When the Fermi surface is spherical, the second term in the bracket in 
eq. (99) vanishes identically. In this case, a xx = e 2 /cFfF/47T7, a xy = (ev-F/2j)a xx and 
Aa xx = — (evp/2'j) 2 a xx according to the RTA, where kp and v-p are the Fermi momentum 
and the Fermi velocity, respectively. Therefore, Ap/po given in eq. (98) becomes zero. 
Except for this special case, the orbital magnetoresistance is always positive [1]. 

The general expression for the magneto-conductivity in Fermi liquids are derived 
in Ref. [38], which is exact of order 0(7~ 3 ). It is derived by performing the analytic 
continuation of eq. (48) for m — 2. This work had enabled us to calculate the magneto- 
conductivity along with satisfying the conservation laws. At low temperatures, the 
magneto-conductivity can be expressed by a simple form: Ao~ xx = Ao~ xx + Ao~ xx ; 

A< x = - -H? • J 5> {-%) 0U 2 ^ ( 10 °) 

AaL =-flJ-jE z * (-jPl ^AcA ( 101 ) 

4 k V fcJ K 7k 

dkx = \v k \-^-, (102) 

/OO J / 
4^ r k,k'( ' e')Pk'(0)rf k ',/7k', (103) 

where L k = Jk/7k and ^k k /(0, e') is the irreducible four-point vertex in the particle- 
hole channel introduced in eq. (60) [40]. A<j^J A in eq. (40) is given by eqs. (100) 
and (102), by replacing J^ with v^. In the FLEX approximation, f p7jj k ,(0, e') ~ 
(3f/ 2 /2)(7rT) 2 Im9x k _k'(^)/^U=o, as shown in eq. (75). 
Here, we can rewrite Ac^ as 







1 

7k : 



(104) 



where Lk = Jk/7k- According to eq. (96), the first term in the bracket in eq. (104) is 
proportional to (^k/^||)coid7coki ^ ^AF7ra?d [35]- Further, remaining terms (the second 
term in eq. (104) and Aa xx ) are also proportional to ^AF7raid f° r a wide range of 
temperatures, since the fey-derivative of 7^ k , yields a factor proportional to £ AF [33, 35]. 
In §5.2, we will explain that |Lk| becomes highly anisotropic below T* due to AF+SC 
fluctuations. For this reason, due to the second term in eq. (104), Ap/p is prominently 
enhanced below T* in under-doped HTSCs. 

Therefore, the magnetoresistance in nearly AF metals behaves as 

Ap/po oc £ AF • pf oc T- 4 (105) 

above T* due to the CVC since po oc T and £ AF oc T _1 . Equation (105) is consistent with 
experimental results in LSCO [12, 71, 72], YBCO [71] and TBCO [73]. On the other 
hand, in a conventional Fermi liquid where £af is constant, eq. (105) gives Kohler's rule 
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Figure 20. ( = (Ap/p ) cot 2 6 H for (i) YBCO and (ii) LSCO calculated by the 
FLEX+CVC method for 0.02 < T < 0.1 (80K - 400K). The modified Kohlcr's rule 
(( =const.) holds well by considering the CVC. On the other hand, Crta obtained by 
the RTA exhibits a strong T-dependence. [Ref. [35]] 



Ap/po oc p 2 . Equations (97) and (105) directly imply the "modified Kohler's rule" [35] 



Ap/po oc tan 2 #H, 



(106) 



where tan# H = °~xy/o~x X = Rn/Po- Experimentally, the modified Kohler's rule is well 
satisfied in various HTSC compounds [12, 71, 73] and in CeMIn 5 (M=Co or Rh) [19, 20], 
whereas the conventional Kohler's rule is completely broken. This long-standing problem 
had been naturally explained in terms of the Fermi liquid theory. The modified Kohler's 
rule implies that the enhancements of i?n and Ap/p are induced by the same origin; 
the CVC due to AF fluctuations. 

Figure 20 shows the numerical results for (i) YBCO and (ii) LSCO by using the 
FLEX+CVC method. The coefficient 



C = (Ap/po) cot 2 H 



(107) 



takes an almost constant value for a wide range of temperatures, only when the CVCs 
are correctly considered: In Fig. 20, we show that Crta = (Ap/po)RTA(cot 2 #h)rta and 
Cmc-rta = (Ap/po)RTACot 2 #h: Crta increases as T decreases, since the anisotropy of 
\h\ = |^k/7k| increases at low temperatures. In contrast, Cmc-rta decreases at low 
temperatures since the i?H is strikingly enhanced by the CVC. Only when both o~ xy 
and Aa xx are calculated using the FLEX+CVC method, the modified Kohler's rule 
C ~const. is satisfied. In YBCO, the dominant contribution to the total magneto- 
conductivity is given by Au^, particularly by the first term in the bracket in eq. (104) 
that includes (d6^/dk\\) 2 . In the case of LSCO, the second term in eq. (104) is very 
important. Further, |Aa"J ~ |A<t£ x |. In both these systems, the modified Kohler's rule 
(106) is well reproduced by the FLEX+CVC method. 
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The experimental value is ( ~ 1.5 ~ 1.7 for YBCO [71] and ( ~ 2 ~ 3 for 
Tl 2 Ba 2 Cu06+<5 [73]. In contrast, ( takes much larger value in LSCO: ( ~ 13.6 for a 
slightly over-doped sample (x = 0.17) [72], and it increases (decreases) as the doping 
increases (decreases) [12]. These experimental results are reproduced by the present 
study and they are shown in Fig. 20. In LSCO, the Fermi surface is very close to the 
van-Hove singularity point (n, 0) since both \t'/t\ and \t" /t\ in eq. (12) are smaller than 
those in other systems. Since v^ = at (ir, 0), the anisotropy of the velocity \v^\ on the 
Fermi surface is larger in LSCO. Hence, the second term in eq. (104) [or eq. (99)] takes a 
large value in LSCO. In heavily over-doped LSCO at 5 = 0.225, where CVC is expected 
to be unimportant, ( exceeds 100 [12]. This result is consistent with recent ARPES 
measurement [136], which shows that the Fermi surface in LSCO passes through (it, 0) 
in the over-doped region 5 = 0.2 ~ 0.22. 

Therefore, experimental value of ( gives us a useful measure of the anisotropy of 
-^k — Jk/lk (^k — ^k/7k m weakly correlated systems). In CeMIn 5 (M=Rh,Co), ( is of 
order O(100) [20], which indicates that the anisotropy of Lk is large in the main Fermi 
surface. 

4-4- Thermoelectric power 

We also discuss the thermoelectric power, S. In HTSCs, S takes a large value in under- 
doped systems, and it increases as T decreases from the room temperature. Except for 
over-doped compounds, S > in hole-doped systems [2, 157, 158, 159, 160, 161, 162] 
whereas S < in electron-doped systems [7, 163] below the room temperature. 
Interestingly, the peak temperature of S is nearly equal to the pseudo-gap temperature 
T* in many hole-doped compounds; in HgBa 2 Cu0 4+( 5 [159, 160], LSCO, YBCO [161] 
and Bi 2 Sr 2 RCu 2 8 (R = Ca, Y, Pr, Dy or Er) [162]. By neglecting the CVC, Hildebrand 
et al. calculated the thermoelectric power for YBCO using the FLEX approximation 
[164]. Here, we study the thermoelectric power both for hole-doped and electron- doped 
systems using the FLEX+CVC method. 

According to the linear response theory [165, 166, 39], the thermoelectric power in 
a Fermi liquid system is given by 

S = a xx /a xx , (108) 

where a xx is the diagonal Peltier conductivity; j x = a xx {—V x T). It is given by 
[165, 166, 39], 

a xx = K2 x (u + i0)/iu>l_ , (109) 

where K^ x (u + iO) is given by the analytic continuation of [39] 

KZA^x) = fj o dre-^(T T £(0)j x (r)), (110) 

where j and j® are charge current operator and heat current one, respectively. The heat 
current is defined by the following energy conservation law; d t (h(r)— /xn(r)) + V-j^(r) = 
0. Then, j contains complex two-body terms in the presence of electron-electron 
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Figure 21. (i) S 1 for YBCO and NCCO given by the FLEX+CVC approximation, 
by including the CVC for Jy. Reflecting the reduction in AF fluctuations, S starts 
to decrease below T* s=a 0.04 (<~ 160) K in the FLEX+T-matrix method, which is 



consistent with the experimental pseudo-gap behavior in S. (ii) S; 
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thermoelectric power given by the inside (outside) of the AFBZ. Stotai = ^in 
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Oo 



correlation. Fortunately, the Ward identity allows us to rewrite the heat current as 
ik ( e ) = e ^k 5 where v^ = Vk(ek + R-e^k) [39]. By performing the analytic continuation 
of eq. (110), we can derive the exact expression of a xx of order 0{E-p/ r i) as [39] 

e V fde(_df°^ 
T 



a. 



£ 



7T 



de 



tv kx {e 



x 



(112) 



(Jkx(e)|G k (e)| 2 - 2^(e)ReG k (e) 2 ) , (111) 

where Jkx(e) is the total current given in eq. (60). At sufficiently low temperatures, a xx 
is approximately simplified as 

tt 2 t r dk\\ i d ( v kx (K)J^(KY 
axx ~ e 3 L{^?z^W{Ei)\dhAW{K)\iAK) 

where — e (e > 0) is the electron charge. dk± represents the momentum vertical to the 
Fermi surface. When 7 is energy- in dependent, eq. (112) yields S = —eTiT/Qz r )a xx = 
—7iTm*/3en (m* = m/z is the effective mass) in the 2D isotropic system with 
£ k = k 2 /2m - ii. 

In Fig. 21 (i), we show the numerical results for S of YBCO and NCCO given by the 
FLEX+CVC method. In LSCO, S > in the optimally-doped case (n = 0.85) below the 
room temperature, and it increases in the under-doped case. In NCCO (n = 1.10), S < 
and it takes a peak around 200K. These results are consistent with the experimental 
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results [2, 157, 158, 159, 160, 161]. The origin of the non-Fermi-liquid-like behavior is 
the prominent k— and e-dependences of 7k (e) due to strong AF fluctuations. According 
to the spin fluctuation theory, d^^E^j/dk^ > at point A in Fig. 4 (i) because 7k 
becomes large near the AFBZ. Since point A is the cold spot of YBCO, eq. (112) 
becomes positive in YBCO. Therefore, S > in under-doped YBCO and LSCO. On the 
other hand, S < in NCCO because d^ k (E^)/dk ± < at point B. In YBCO, S slightly 
decreases if the CVC is included, since | Jk| < \vk\ around the cold spot. In NCCO, in 
contrast, absolute value of S is enhanced by the CVC. In both cases, the effect of the 
CVC on S is smaller than that on R^. Figure 21 (ii) shows the S for YBCO (n = 0.90), 
together with Si n (S^ut) which denotes the contribution from the inside (outside) of the 
AFBZ. S[ n is positive and it deviates from a conventional behavior (S oc T) since the 
e-dependence of 7k (e) becomes large inside of the AFBZ at lower temperatures. On the 
other hand, S ont is negative and it approaches zero as T decreases. Therefore, the sign 
of S — Si n + Sout becomes positive below T = 0.075 (~ 300 K). 

4-5. Summary of this section and comments on other transport coefficients 

In HTSCs, anomalous transport phenomena (such as the violation of Kohler's rule) had 
been frequently considered as strong evidence for the breakdown of the quasiparticle 
picture. For example, the RTA for the highly anisotropic t^ model cannot reproduce 
relationships (l)-(4) for hole-doped systems at the same time for a wide range of 
temperatures. Furthermore, the RTA cannot explain the negative Hall coefficient 
in electron-doped systems, since the curvature of the true Fermi surface is positive 
everywhere. The explanation for the nearly symmetric behavior of Rn and S with 
respect to the type of carrier doping (shown in Fig. 3) was desirable for a long time. 
These highly nontrivial transport phenomena had been one of the central issues in 
HTSCs, which should serve to elucidate the actual electric ground state of HTSCs. 

To resolve this long-standing problem, we developed a method to calculate various 
transport coefficients based on the microscopic Fermi liquid theory. In the RTA, 
the momentum and energy transfers between the quasiparticles by scattering are not 
treated correctly. Therefore, the RTA results frequently yield erroneous results. To 
overcome this defect, we study the role of the CVC in nearly AF Fermi liquids. 
We find that the total current J^ shows an anomalous k-dependence because of the 
CVC in nearly AF metals, which gives rise to the non-Fermi liquid behaviors such as 
i?H oc £af c* r ^~ 1 an d Ap/p oc P^ 2 £af ^ T~ 4 - Consequently, the modified Kohler's rule 
Ap/po oc (i?n/p) 2 oc T~ 4 is realized in HTSCs. We also studied the S for HTSCs, 
which increases as the temperature decreases above T*. This experimental fact is 
reproduced by considering the strong k— and e-dependences of 7k(e) [36]. Therefore, 
various anomalous transport phenomena in HTSCs above T* can be well explained by 
the FLEX+CVC method. 

Recently, Tsukada et al. measured the Ru in clean heavily over-doped LSCO 
samples [170], and found that Rn becomes almost temperature-independent for 5 > 0.24; 
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further, its sign smoothly changes from positive to negative between S = 0.28 and 0.32, 
corresponding to the change in the curvature of the Fermi surface. Hussey et al. also 
showed that the RTA analysis is successful in heavily over-doped T12201 [149] and LSCO 
[150]. These experimental results are consistent with the present theory: According to 
the present analysis, the RTA works well in heavily over-doped systems since the CVC 
is unimportant when the AF fluctuations are very weak. This fact is well explained in 
Ref. [167]. 

Finally, we comment on other important theoretical study involving the CVC in 
strongly correlated systems. In the interacting electron gas model, the electron cyclotron 
frequency u c is given by eH/mc, where m is the bare electron mass. In 1961, Kohn 
proved that u c is unchanged by electron-electron interaction due to the consequence of 
the angular momentum conservation law, which is called the "Kohn's theorem" [171]. 
The Kohn's theorem is also proved by the phenomenological Fermi liquid theory, by 
correctly considering the CVC [172]. Kanki and Yamada studied this problem based 
on the microscopic Fermi liquid theory, and found that uj c is influenced by the electron- 
electron correlation if the Umklapp processes are present [173]. The Umklapp processes 
also reduces the Drude weight in o~{u) [174] and the penetration depth in the SC state 
[175]. Moreover, the effect of the vertex correction on the Raman spectroscopy in HTSC 
was studied [176]. 

5. Transport phenomena in HTSCs below T* 

5.1. Mechanism of pseudo-gap phenomena 

In under-doped HTSCs, a deep pseudo-gap in the DOS emerges at the chemical potential 
below T* ~200 K. The origin of the pseudo-gap has been a central issue with regard 
to HTSCs. The pseudo-gap does not originate from the spin fluctuations since the 
strength of the spin fluctuations decrease below T* [78, 79, 80, 81, 82]. According 
to recent ARPES measurements [83, 84, 85, 86], the k-dependence of the pseudo-gap 
coincides with that of the d x 2_ y 2 -wave SC gap function. This pseudo-gap structure in 
the quasiparticle spectrum p(k, u) = ImG(k, uj — i5)/ir starts to appear around u = 
below T*. Below T c , the quasiparticle spectrum shows sharp peaks at the edge of the 
SC gap {uj ~ A) since the inelastic scattering is reduced by the SC gap. At the same 
time, a dip-hump structure is induced around uud (> A) by the strong resonance peak 
in \mx s {uo) at uo r ~ 40 meV in YBCO and Bi 2 Sr 2 CaCu 2 8 (BSCCO). We note that 
uid ~ A + uj r and uj t < 2A. However, the overall gap-like structure in the quasiparticle 
spectrum remains unchanged around T c . 

Motivated by these experimental facts, various strong-coupling theories of SC 
fluctuations have been studied to reproduce the pseudo-gap phenomena [87, 89, 51, 88]. 
Several groups developed the "FLEX+T-matrix theory" , where the normal self-energy 
correction induced by the strong SC fluctuations, S^ CF (e), has been self-consistently 
included into the FLEX approximation. The SC "amplitude" fluctuations are given by 
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Figure 22. (i) T-matrix induced by the AF fluctuations, (i) Self-energy Sk(w) 
in the FLEX+T-matrix approximation, (iii) DOS given by the FLEX+T-matrix 
approximation. The pseudo-gap appears below T* ~ 0.04 due to the strong SC 
fluctuations, (iv) Im££ FF (-i<5) and ImS^ CF (-i(5) given by the FLEX+T-matrix 
approximation at T = 0.02. For comparison, ImE^—iS) given by the FLEX 
approximation at T = 0.02 is also shown. 



the T-matrix (particle-particle scattering amplitude) induced by the AF fluctuations, 
which prominently increase below T* [51, 52, 88, 89, 37]. The T-matrix with respect 
to the AF fluctuations, which is responsible for the Thouless instability for the d x 2_ y 2- 
channel, is given by 

^k,k';q(en, ^n r , ^l) = V k _ p ( e n — e n') + T ^ V k _ k , (e n — € m ) 

x G , p(e m )G_ p+q (— e m + u;i)T Pi k' ; q(e m , e n i\ ui), (113) 

which is shown in Fig. 22 (i). V k AFF is given in eq. (15). In the FLEX+T-matrix 
approximation, the Green function and the self-energy are given by 



G k (e„) 

S AFF (eJ = eq. (14) 



K + /i-e k -Sf F (e r 



^SCF, 



en)}' 1 , 



k \ e n) 

^k ( 6 n) = T £ , ( J-k+q(~ e n + ^z)Tt,k;q(£n; ^n] ^>l), 
1,1 



(114) 
(115) 
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where S^ FF and S| CF are shown in Fig. 22 (ii). In the self-consistent FLEX+T- 
matrix approximation, we have to solve eqs. (14), (15), and (113)-(115) self- 
consistently. Unfortunately, it is very difficult to calculate the T-matrix in eq. (113) 
since a lot of memory in computer is required. Fortunately, eq. (113) can be 
simplified by considering only the d^.^-channel and drop other pairing channels as 
Tk >k '-q(e n ,e n >;uji) « ^(e n )^>(e n ')tq(u}i), where ^ k oc cosfc x - cosA;^ [51, 52, 88, 89, 37]. 
This approximation is expected to be reasonable when the temperature is close to T c in 
HTSCs. Then, the k-dependence of £ k CF becomes £ k CF oc T kk;q oc ?/> k - 

In the present approximation, T* is defined as the temperature below which 
S| CF takes considerable values. T* is slightly greater than the T C FLEX in the FLEX 
approximation. In a pure 2D system, T c = since the FLEX+T-matrix theory 
satisfies the Mermin-Wagner-Hohenberg theorem. T c becomes finite when weak three- 
dimensionality is assumed. Sufficiently above T*, where S| CF m 0, the FLEX+T-matrix 
theory is equivalent to the FLEX approximation. The obtained DOS and the self-energy 
are shown in Fig. 22 (iii) and (iv), respectively. Below T* ~ 0.04, a large pseudo-gap 
structure appears in the FLEX+T-matrix approximation, due to the cooperation of 
the real part and the imaginary part of S^ CF [37]. The AF fluctuations are suppressed 
below T* in accordance with the pseudo-gap in the DOS [88, 52, 51, 89]. For this reason, 
ImE k vFF (0) in the FLEX+T-matrix approximation is much smaller than ImE k (0) in the 
FLEX approximation, as recognized in Fig. 22 (iv). Below T c , the pseudo-gap in the 
FLEX+T-matrix approximation smoothly changes to the SC gap [87]. 

Here, we discuss transport phenomena in the pseudo-gap region. In YBCO [9, 74], 
Bi2201 and Bi2212 [177], i?H shows a maximum around the pseudo-gap temperature 
T*, at which \jT{T shows the maximum value. Also, the peak temperature of S is 
nearly equal to T* in HgBa 2 Cu0 4+( 5 [159, 160], LSCO, YBCO [161] and Bi 2 Sr 2 RCu 2 8 
(R = Ca, Y, Pr, Dy and Er) [162]. As discussed in Ref. [33], these behaviors are 
naturally understood since both of them are strikingly enhanced by the AF fluctuations 
above T*. Therefore, both Rn and S should decrease below T* in accordance with 
the reduction in the AF fluctuation. This behavior is easily reproduced by the 
FLEX+T-matrix method [33, 51]. The obtained p, i?n and S for hole-doped systems 
by using the [FLEX+T-matrix] +CVC method are shown in Figs. 18 (i), 19 (i), 
and 21 (i), respectively. In this approximation, T k l k ,(e, e') in eq. (60) is given by 
i [cth^=r + th^p] Im{Vk_ k /(e — e') — T kik;k _ k /(e, e; e — e')}, where T kk;k _ k / represents the 
T-matrix in eq. (113). The detailed method of calculation is explained in Ref. [37]. 

As shown in Figs. 19 (i) and 21 (i), both i?n and S start to decrease below T* ~ 0.04 
(~160 K) in the FLEX+T-matrix approximation. These results are consistent with 
experiments [159, 160, 161]. Figure 18 (i) shows the resistivity given by the FLEX 
approximation and the FLEX+T-matrix approximation, both of which are given by 
including CVCs. As pointed out in Ref. [167], p given by the FLEX+CVC method 
shows a tiny "kink" structure at ~ T because of the CVC. This kink structure at 
To ~ 0.12 is experimentally observed in LSCO [185] and YBCO [186]. The kink becomes 
more prominent in the [FLEX+T-matrix] +CVC method since the inelastic scattering 
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is reduced by the formation of the pseudo-gap. As shown in Fig. 22 (iv), ImS^ FF (0) 
decreases prominently below T*, in accordance with the emergence of ImS| CF (0). As a 
result, p is approximately proportional to T 2 below T*. To summarize, p, Rn and S are 
suppressed in the pseudo-gap region in the [FLEX+T-matrix]+CVC method. 

According to the FLEX+T-matrix theory, the SC "amplitude" fluctuations is the 
origin of the pseudo-gap formation. Other than this scenario, various mechanisms of 
pseudo-gap formation have been proposed. For example, Kivelson and Emery proposed 
the SC "phase" fluctuation scenario [178]: they considered that the amplitude of 
the SC order-parameter develops in the pseudo-gap region, whereas the global phase 
coherence is absent. Moreover, several hidden-order scenarios, such as d-densitj wave 
formation [179], have been proposed. To elucidate the correct scenario, anomalous 
transport phenomena in the pseudo-gap region are significant and they severely constrain 
the theories. In the following, we show that the scenario involving SC "amplitude" 
fluctuations enables us to understand the anomalous transport phenomena below T* in 
a unified way. 

5.2. Enhancement of Nernst coefficient and magnetoresistance 

In contrast to Ru and S, the Nernst coefficient v [y = S yx /H z = —E y /H Z V X T; 
off-diagonal thermoelectric power under a magnetic field) prominently increases in 
the pseudo-gap region by approximately 100 times than that in conventional metals 
[180, 181, 182]. Moreover, in the pseudo-gap region, the magnetoresistance Ap/p 
increases faster than T~ 4 [75, 76]. Therefore, the modified Kohler's rule in eq. (106) 
violates below T*. Such nontrivial behaviors of v and Ap/p have attracted considerable 
attention as a key phenomenon closely related to the origin of the pseudo-gap. 

Here, we study the Nernst coefficient. According to the linear response theory 
[165, 166, 39], v is given by 

v = [a xy /cr xx - Sa xy /cr xx ] /H z , (116) 

where a xx = j x /(—'V x T), and a xy = j x /(—V y T) is the off-diagonal Peltier conductivity 
under the magnetic field H z . In conventional metals with simple Fermi surfaces, v is 
small because of an approximate cancellation between the first and the second terms in 
Eq.(116), which is known as the Sondheimer cancellation. (Sondheimer cancellation is 
exact only when e£ = k 2 /2m.) In HTSCs, however, Sondheimer cancellation is totally 
violated since a xy is considerably enhanced. That is, a xy is the origin of giant Nernst 
coefficient in LSCO below T* and in NCCO below T [180, 181]. 

In this study, we investigate a xy due to the quasiparticle transport, and find that 
a xy is prominently enhanced below T* if we include the CVC caused by SC and AF 
fluctuations. As a result, v ~ ot xy ja xx takes a large value below T*. Moreover, we show 
that all the transport anomalies below T* are understood as the quasiparticle transport 
phenomena in a unified way. Note that Yip showed that the Maki-Thompson-type 
CVC due to SC fluctuations disappears in <i-wave superconductors when the inelastic 
scatterings are negligible [184]. However, this is inappropriate for HTSC since inelastic 



Anomalous transport phenomena in Fermi liquids with strong magnetic fluctuations 51 




---,.Z2 



Nemst Coefficient by FLEX-T-matiix 




60 



40 



(ui) 
Nernst Coefficient for MCCO 



in-L \5T 






Stanth/B ^ 




in-L15} 



iM 



10» 



200 30J 

T[K] 



400 



Figure 23. (i) Ap/p Q for LSCO obtained by the [FLEX+T-matrix]+CVC; The 
modified Kohler's rule Ap/p Q ex tan 2 6 H holds above T*. Below T* ~ 0.04 (~ 160) 
K, modified Kohler's rule is violated since Ap/po increases quite drastically that is 
caused by the CVC due to AF+SC fluctuations, (ii) v for LSCO obtained by the 
[FLEX+T-matrix]+CVC; the abrupt increase in v below T* is caused by the CVC 
due to AF+SC fluctuations. Experimental data ^ cxp is cited from [180]. (iii) v for 
NCCO is obtained by using the FLEX+CVC approximation. Experimental data j/ cxp 
is cited from [7]. [Rcf. [32]] 



scatterings are much larger than elastic scatterings. We find that the Maki-Thompson 
type CVC is significant in HTSC. 

We note that Ussishkin et al. studied v induced by short-lived Cooper pairs, which 
is independent of the quasiparticle lifetime and order 0(7°) [183]. It will be observed in 
the close vicinity of T c . However, quasiparticle transport is dominant in in good metals 
where k-pl ~ E-p/'y ^> 1. (Here, / = Vk F r is the mean free path.) In HTSCs, kpl = 1 
corresponds to 1700 fiQcm if we omit the k-dependence of l^. Therefore, quasiparticle 
transport will be dominant in slightly under-doped HTSC (p < 200 /iflcm above T c ) 
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for a wide range of temperatures. Hereafter, we show that a xy due to quasiparticle 
transport is strongly enhanced by the CVC. 

According to the microscopic Fermi liquid theory, a xy is given by the correlation 
function between the heat current and the charge current in the presence of H z ^ 0. 
Therefore, a xy is given by eq. (50) by replacing j M with j®. After the analytic 
continuation, the exact expression for a xy of the order 7 -1 is given by [39] 

e 2 ^ fdef d/°V T ~ R ,, M _ R ,, |2 



k 

x |4(e)| 7 k(e)A k (e), (117) 

M4 = (q k (e) x ^0\ , (118) 

where L k (e) = J k (e)/7 k (e), <7k(e) = e • Vk is the quasiparticle heat velocity, and Qk(e) is 
the total heat current [39, 37]. We stress that the CVC term in eq. (119) vanishes if we 
omit the energy dependence of 7^ k ,(e,e'). This fact means that the heat CVC is small 
and thus Qk(e) ~ <fk(e) in general cases [39, 37]. In the FLEX+T-matrix approximation, 
T$,(e, e') is given by i [cth^=r + th^] Im{Vk_ k /(e - e') - T k , k;k _ k /(e, e; e - e')}. 

In contrast to Ru and S, the Nernst coefficient v [180, 181] and Ap/p [75, 76, 12] in 
LSCO rapidly increase below T*. These experimental facts are also reproduced by the 
[FLEX+T-matrix] +CVC method as shown in Fig. 23 (i) and (ii). Since i?H decreases 
whereas Ap/p increases drastically below T*, the plot of Ap/p as a function of tan 2 #h 
forms an "inverse S-shape". Therefore, the modified Kohler's rule is completely violated 
below T*. This result is very similar to the experimental results provided in Ref. [12]. 
Next, we discuss v in NCCO using the FLEX (not the FLEX+T-matrix) approximation 
since SC fluctuations is considered to be absent in NCCO. The obtained v is shown in 
Fig. 23 (iii); it takes a large value only when the CVC due to the AF fluctuations 
are considered. In NCCO, v increases gradually as T decreases, and it starts to 
decrease below the maximum temperature of 120 K. The obtained behaviors of v are 
semiquantitatively consistent with the experimental results [7]. Therefore, the present 
numerical study can explain the experimental behaviors of v and Ap/p in HTSCs both 
above and below T*. 

Here, we explain a theoretical reason why v is enhanced in both LSCO and NCCO: 
A k in eq. (118) can be rewritten as 

A k = Q k L k (dei/dk\\) cos(0j[ - 9?) + Q k (<9L k /<%||) sin(^ - #g), (120) 

where 9£ = tan _1 (J k?; /J ka .), 9^ = t&n. -1 (Qk y / Qk x ) , an d L k = J k /7 k is the mean free 
path with CVC. In the RTA, the heat current is given by q\ = e ■ £T k . Since #£ = k , 
the second term is absent in the RTA. If one includes the CVC, Q k is not parallel to 
J k because the CVC for q k is usually small (i.e., Q k ~ q k ) as we will discuss below 
[39, 37]. Therefore, we obtain 9^ 7^ 0^. In this case, the second term gives rise to an 
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enhancement of v if L^ is highly anisotropic around the cold spot. Figure 17 (ii) shows 
Lk in LSCO given by the FLEX+CVC approximation. We can see that dL^/dkn takes 
a large value near the cold spot below T = 0.02. (Note that dLk/dk\\ = just at point 
A.) This is the origin of increment in z^flex i n Fig- 23 (ii). As shown in Fig. 23 (iii), 
z^lex for NCCO is much larger than z/flex for LSCO: One reason is that the anisotropy 
of L^ LEX in NCCO is greater since £af ^ 1 in NCCO. Another reason is that 7 k 1 (oc v) 
in NCCO is larger than that in LSCO. For these reasons, FLEX+CVC method can 
explain the huge v observed in NCCO. 

Next, we consider the increment of v in the pseudo-gap region in LSCO using the 
FLEX+T-matrix method, which is shown as ^flex+t in Fig. 23 (ii). The CVC due 
to SC fluctuations (Maki-Thompson term) represents the acceleration of quasiparticles 
caused by the short-lived Cooper pairs. Since HTSCs are J^.^-wave superconductors, 
the CVC due to the SC fluctuations magnifies Jk except at the nodal point (point A in 
Fig. 14). In fact, when the Ginzburg-Landau correlation length is longer than £af> the 
total current in the FLEX+T-matrix approximation is approximately given by [37] 

t ti , u \ ^k + «kW k * , . 

J k « (1 + hk) 2 , (121) 

K = ImS| CF (0)/ImSf F (0), (122) 

where ak = (1 — c^f) has been introduced in eq. (77). h^ represents the enhancement 
factor due to SC fluctuations. The anisotropy of L k = | Jk/7k| becomes considerably 
prominent below T* since ImS^ CF (0) oc ■?/>£ increases whereas Im£k FF (0) decreases below 
T* as shown in Fig. 22 (iv). For this reason, the second term in eq. (120) takes a large 
value in the neighborhood of the point A, as explained in Ref. [37]. Thus, [FLEX+T- 
matrix] +CVC method can explain the rapid increment in v in LSCO below T*. 

If the AF fluctuations are absent, 6£ = 9^ = #£ even if the CVC due to SC 
fluctuations are taken into account. Therefore, to explain the enhancement in v in the 
pseudo-gap region, we have to include CVCs due to both the AF fluctuations and SC 
fluctuations. In contrast, i?n decreases below T* in proportion to £ AF since the second 
term in eq. (120) is absent if we replace Qk with Jk. 

Finally, we explain why the heat CVC is usually small. As we have explained, 
the CVC represents the current of the Fermi sea that is transfered from the excited 
quasiparticles by the electron-electron scattering. The charge CVC can be large because 
of the momentum conservation law. However, the heat current q\ = (e k — /x)vk is not 
conserved even in the free dispersion model, that is, q\ + q^ ^ gk +q + <?k'-q even if 
e k + e k' = e k+ q + e k'- q - F° r this reason, the heat CVC is small in general, and therefore 
the thermal conductivity k, is finite even if the Umklapp process is absent: In a 3D 
free dispersion model, k, = (9/8)k rta within the second-order perturbation theory with 
respect to U [39] . The heat CVC is also small in nearly AF metals: Since the direciton 
of the the heat CVC due to the quasiparticle at k + q in Fig. 8 (the MT term) depends 
on the sign of 6k+q — t 1 -, the heat CVC becomes small due to cancellation after the 
q-summation for |q — Q| < ^f- 
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5.3. Summary of this section 

In conclusion, we studied p, Ru, Ap/p , S and v in HTSCs using the FLEX+T-matrix 
approximation. The results are shown in Figs. 18 (i), 19 (i), 21 (i), and 23 (i)-(ii). 
We can explain that (i) R^ and S start to decrease below T*, (ii) p starts to deviate 
from the T-linear behavior at To, and it is approximately proportional to T 2 below T*. 
Moreover, (iii) v and Ap/p further increase below T*. Therefore, the present study 
gives us a unified understanding of the various anomalous transport phenomena in both 
below and above T*. 

This study provides strong evidence that the SC fluctuations are predominant in the 
pseudo-gap region. The striking increment in v below T* is a cooperative phenomenon 
between the d-wave SC fluctuations and AF fluctuations. We stress that the RTA cannot 
explain the enhancement of v since the second term in eq. (120) vanishes identically 
in the RTA. Further, a very large Nernst signal appears also in electron-doped systems, 
which starts to increase in proportion to T^ 1 below To ~ 600 K. This experimental 
fact is reproduced by the CVC due to the AF-fluctuations, even in the absence of SC 
fluctuations. In §5.2, we found that the origin of rapid increment in v is the second term 
of eq. (120), which emerges only when the CVC is included. In NCCO (and CeMIns), v 
starts to increase below T because of the CVC due to strong AF fluctuations; £af ~^> 1- 
In LSCO, on the other hand, increment in v below To is small since £af is of order 1; 
v in LSCO starts to increase rapidly below T* <C T , with the aid of the strong SC 
fluctuations. 

Here, we comment on the transport phenomena in Nd-doped LSCO, 
Lai.4_ a; Ndo.6Sr- r Cu04, which shows the static stripe ordered phase at T st ~ 80 K [187]. 
In this compound, p increases whereas i?H decreases monotonically below T st , which 
indicates that the system becomes one- dimensional due to the stripe order. This ex- 
perimental fact might tempt us to consider that the stripe order is the origin of the 
reduction in R^ in other HTSCs in the pseudo-gap region. However, the stripe order 
scenario cannot explain the reduction in p nor increment in v below T*. On the other 
hand, the overall transport phenomena in HTSCs can be explained by considering the 
CVCs due to AF+SC fluctuations. 

In this section, we explained that the pseudo-gap behaviors in p and Rn originate 
from the reduction of AF fluctuations, and the pseudo-gap behavior of v reflects the 
increment of SC fluctuations. The pseudo-gap temperature determined by p, T* is 
larger than T^ H and T* determined by i?n and u, respectively. Recently, impurity (Zn) 
effect on the pseudo-gap behavior in YBCO was studied [188], and found that T* and 
Tft H are independent of impurities, whereas T* decreases with Zn-doping in parallel to 
T c . It is an important future problem to show whether these experimental facts can be 
explained by the FLEX+T-matrix approximation or not. In §7, we will show that the 
CVC due to fluctuations is usually sensitive to the impurities. Therefore, to reconcile 
this problem, we have to study the impurity effect on the CVC seriously. 

It should be noted that v takes a very large value in the vortex-liquid state above 
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Figure 24. i?H(w) in optimally-doped YBCO at 95 K, derived from Faraday rotation. 
The closed circles and squares represent experimental values of RcRh(w) an d Imi?H(w), 
respectively. [Ref. [192]] 



H c i in a clean 2D sample, reflecting the high mobility of the vortices. Therefore, v is 
frequently used as a sensitive probe for the mixed state. Based on his observation, Ong et 
al. proposed that spontaneous vortex-antivortex pairs emerge in under-doped systems 
below T*, and they govern the transport phenomena in the pseudo-gap region [180]. 
However, this assumption seems to contradict with the other transport coefficients; for 
example, the flux-flow resistance does not appear below T*. 

6. AC transport phenomena in HTSCs 

In previous sections, we discussed the DC transport phenomena in nearly AF metals, 
and found that the CVC produces various non-Fermi-liquid-like behaviors. In principle, 
the AC transport phenomena can yield further useful and decisive information about 
the electronic status. Unfortunately, the measurements of the AC transport coefficients 
are not common because of the difficulty in their observations, except for the optical 
conductivity cr xx {u)) measurements. 

Fortunately, Drew's group has performed intensive measurements of the AC Hall 
coefficient R H (u) = a xy {oj) / <j 2 xx {u) in YBCO [189, 190, 191, 192], BSCCO [193], LSCO 
[194], and PCCO [195]. They found that the w-dependence of Rr(uj) in HTSC shows 
amazing non-Fermi-liquid-like behaviors, which have been a big challenge for researchers 
for a long time. Here, we show that this crucial experimental constraint is well satisfied 
by the numerical study using the FLEX+CVC method. 

In the RTA, both a xx (oj) and o~ xy (u) in a single-band model will follow the following 
"extended Drude forms": 



o 



RTA, 



a;) = Q xx (2-f (uj) 



iz 1 uj) 1 , 



(123) 



Anomalous transport phenomena in Fermi liquids with strong magnetic fluctuations 56 



(") 




LSCO ( n=0. m 



Figure 25. (i) <r xy {uj) and (ii) R u (u>) for LSCO (n = 0.9, U = 5) given by the 



FLEX+CVC method. 



0.1 corresponds to ~ 300cm . Because of the CVC, 



a xy (u) >• a^y A (uj) at low frequencies. [Rcf. [196]] 






a 



^(2 7o (a;) - mT^)" 2 , (124) 

where z _1 is the mass-enhancement factor and 70 (u;) is the cu-dependent damping rate 
in the optical conductivity, which is approximately given by 70 (w) ~ (7coid(w/2) + 
7coid(— ^/2))/2 for small ui. According to the spin fluctuation theory [56], 70(0;) oc 
max{cj/2, 7rT}, which is observed by the optical conductivity measurements. The in- 
dependence of z is important in heavy-fermion systems {1/z ^> 1 at u — 0), whereas 
it will not be so important in HTSC since 1/z is rather small. Expressions (123) and 
(124) are called the "extended Drude form". Within the RTA, the AC-Hall coefficient 
is independent of u even if the cu-dependence of z is considered: 



/C>) = V xy /n 2 xx ~ Ijne. 



(125) 



Very interestingly, Drew's group has revealed that Rn(w) in HTSC decreases 
drastically with uo: as shown in Fig. 24, \mR^{uj) shows a peak at ~ ujq = 50 cm -1 



in 



optimally-doped YBCO [192]. Moreover, lmRn(ou) is as large as ReRu(uj) for u > u , 
as a consequence of the Kramers- Kronig relation between ReR}i(uj) and IhlRh (<*>). Such 
a large u;-dependence of i?n cannot be explained by the RTA, even if one assume an 
arbitrary (k,u;)-dependence of the quasiparticle damping rate 7k(^)- Therefore, the 
AC-Hall effect severely constrains the theories involving the normal state of HTSCs. 

Recently, we studied both cr xx (u) and a xy (uj) in HTSC using the FLEX+CVC 
method, by performing the analytic continuation of eqs. (51) and (53) using the Pade 
approximation [196, 197]. Since the 10- dependence of the CVC is correctly considered, 
the obtained a xx {oS) and o~ xy {uj) satisfy the /-sum rule very well - - / °° Rea xx (u)du = 
Tie 2 J2k(d 2 e k/ dkx) n k and J °° Rea xy {ui)dui = 0; the relative error is less than 2.5% [196]. 
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Note that a useful /-sum rule for the Hall angle has been derived in Ref. [198]. The 
obtained numerical results are shown in Fig. 25. We find that (a) o~ xx {uj) m cr^^iuj) 
follows eq. (123), by substituting 70 (u;) = (7 co id(^/2) + 7coid( — u>/2))/2. Therefore, 
the reduction in the Drude weight in a xx (yj) due to the CVC is small within the 
FLEX+CVC method. However, we stress that the /-sum rules for both <r^J A (a;) and 
<t^J a (cj) do not hold since the CVC is ignored [196]. (b) R^ fa (uj) is a nearly real 
constant of the order 1/ne since [o"^J A (w)] 2 oc o~ x x £ a {uj) is well satisfied, (c) o~ xy {uS) 
obtained using the FLEX+CVC method completely deviates from the extended Drude 
forms in eq. (124). The effect of CVC on o~ xy (u) is very large for the far-infrared region 
uj < 0.03 (~ 100cm -1 ), since the CVC due to AF fluctuations is strong for uj < uj s { ~ T. 
Therefore, a xy (u) and Rh_(ou) obtained using the FLEX+CVC method show striking 
^-dependences as shown in Fig. 25. In Fig. 25 (i), the peak frequency of lma xy (uj), u xy , 
is six times less than that of Imcr^u;), uj xx , which is consistent with experiments. This 
fact cannot be explained by using the extended Drude form. 

Here, we analyze a qualitative cu-dependence of o~ xy {uj) by considering the in- 
dependence of the CVC. The Bethe-Salpeter equation for uj 7^ is given by 






k'q 

x G£ Hl (c' + a;/2)G^ q (c'-a;/2)J k+q (e , ;a;), (126) 

which is equivalent to eq. (60) when uj = 0. A simplified Bethe-Salpeter equation for 
uj = is given in eq. (78). Here, we extend eq. (78) to the case of uj 7^ 0. By noticing 
the relationship Gf +q (e' + cj/2)G£ +q (e' - uj/2) « 7rp k (e') • 2/(2 7k - iujz' 1 ), eq. (126) is 
simplified for \uj\ <C 7 as 

JkM = 4 + 9 ak " 27k r J k *M, (127) 

27 k — iujz 1 
where (1 — ctk) -1 oc £^f- In deriving eq. (127), we assumed that the e-dependence of 
X q (e + iS) is small for |e| < 7. The solution of eq. (127) is given by 



Jv 



(27 — iujz 



-l*- 



«k • 27 _ 

u k + 7. : TT^k* 

27 — IUJZ 1 



:i28) 



(2 7 - iujz- 1 ) 2 - (a k • 2 7 ) 2 
which is equivalent to eq. (79) when 00 — 0. Similar to the derivation of eq. (96), we 
get 

r <9Jk (27 -iuz- 1 ) 2 w dv k 

J* x — = \ / , 2 • ^k x — -. 129 

ok\\ (27 — iujz l ) z — (a ■ l^y ok\\ 

As a result, o~ xy {oj) for uj <C 7 is approximately given by 

a xy (uj)/H z = e ^ J^ ^—ZTy (* * ^ 

((l-a k )2 7o -^2 ; - 1 )((l + a k )27o-^2 ; - 1 )- v ; 

If we omit the CVC (i.e., a k = 0), eq. (130) is equivalent to the extended Drude 
form in eq. (124). According to eq. (130), the solution of Re<j xy (uj ) = is given by 
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Figure 26. (i) w-dependence of \J B ja xy in BSCCO at 300K. Open squares 
and diamonds represent real and imaginary part of theoretical result given by the 
FLEX+CVC method, (ii) y/B/a xy in optimally-doped BSCCO at 950 cm" 1 . The 
solid dots represent the result of the FLEX+CVC method. [Ref. [193]] 



ujq = 2v/l — 0:^.-27, which is much smaller than the solution obtained from the RTA; 
UJ o' TA — 2-27- This analytical result can be recognized in the numerical result shown in 
Fig. 25. As a result, the striking u;-dependence of Ra(w) in HTSCs for \u\ < 70 can be 
explained by including the CVC. 

Furthermore, the effect of CVC is still large even for the infrared region u < 0.3 (~ 
1000cm -1 ), which is consistent with the experimental results. This fact requires the 
satisfaction of the /-sum rule, L dua xy (u) = 0, since the DC a xy is significantly 
enhanced. Schmadel et al. [193] observed a xy {uj) in optimally doped BSCCO for a 
wide range of u (< 0.3 ~ 1000cm -1 ), and found that a xy {oj) follows a "simple Drude 
form"; a xy oc (2-f xy — izu})~ 2 , where w-dependence of ^ xy is much smaller than 7o(cu) and 
Ixy "C 7o(0). The experimental u;-dependence of y/B/a xy as well as theoretical result 
are shown in Fig. 26 (i). Moreover, Fig. 26 (ii) shows that Ke{ax y } oc j xy depends 
on T sensitively. Therefore, ^ xy is independent of u in the infrared region. In terms of 
the RTA, this experimental result highly contradicts the fact that 70 is independent of 
T for cu ^> T; 70 ex max{w, ttT}. This mysterious behavior is well reproduced by the 
FLEX+CVC method, as indicated by the solid dots in Fig. 26 (ii). This numerical result 
suggests that the effect of the ^-dependence of 70 (o>) on the AC-Hall conductivity, by 
which a^y (u>) deviates from the simple Drude form, approximately cancels with that 
of the CVC. 

Interestingly, the Hall angle 6r(uj) = a xy {uj)/a xx {uj) follows an approximate simple 
Drude form in both YBCO [190] and LSCO [194] - - 6 H (cu) oc (2 7h - izuj)~ l with 
an cj-independent constant 7h oc T~ n (n = 1.5 ~ 2). [Exactly speaking, both a xy 
and #h cannot follow the simple Drude forms at the same time when a xx exhibits an 
extended Drude form.] This experimental fact has also been theoretically reproduced 
[196]. Therefore, the anomalous (a;, Independences of cr xy (uj) and Rr(u>) in HTSCs can 
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be semiquantitatively explained by the FLEX+CVC method --a microscopic theory 
without any fitting parameters except for U. At the present stage, the FLEX+CVC 
method is the only theory with the capability to explain the anomalous AC and DC 
transport coefficients in a unified way. We briefly comment on the carrier-doping 
dependence of Qr = —co/lra{9 }l 1 (uj)}: The experimental value of Qn/H z for optimally- 
doped YBCO (n ~ 0.85) is 0.15 [cm _1 T" 1 ] and it increases to 0.3 [cm~ 1 T~ 1 ] in slightly 
under-doped YBCO and BSCCO (n ~ 0.9). This experimental doping dependence 
can be quantitatively reproduced by using the FLEX+CVC method for n < 0.9 
[196]. However, this method cannot explain the large experimental value VLyi/H z rs 0.5 
[cm _1 T _1 ] in heavily under-doped YBCO (T c = 0) [189]. This experimental fact will 
give us an important hint to understand the essential electronic properties in heavily 
under-doped HTSCs. 

Recently, the infrared Hall coefficient in PCCO (5 = 0.12 ~ 0.18) was measured 
for io = 1000 ~ 3000 cm -1 [195]. The observed ^-dependences for 5 < 0.15 cannot be 
explained by simple extended Drude models. They show that the experimental results 
can be understood by assuming a SDW gap Asdw ~ 1000 cm -1 . It is an important 
future problem to study the role of the CVC in the infrared Hall coefficient in PCCO. 

7. Impurity effects in nearly AF metals 

In earlier sections, we studied the transport phenomena in nearly AF metals without 
randomness. By including CVCs, we succeeded in explaining the anomalous DC and 
AC transport phenomena in a unified way, in both hole- and the electron-doped HTSCs. 
However, we have neglected the impurity or disorder effects on the transport phenomena, 
although they are prominent in real under-doped HTSCs. In fact, the residual resistivity 
due to disorder increases drastically as the system approaches the half-filling [156]. 
Moreover, STM/STS measurements revealed that the electronic states in under-doped 
HTSCs are highly inhomogeneous at the nanoscale, reflecting the random potential 
induced by the disordered atoms outside of the Cu0 2 plane [199]. These experimental 
facts are reproduced by assuming that the SC pairing potential is strongly influenced 
by the impurity potential [200, 201] 

Further, the anomalous transport phenomena near the AF QCP are sensitive to 
randomness. For example, within the Born approximation, the CVC due to electron- 
electron interaction vanishes at zero temperature since only elastic scattering exists at 
T = 0. In fact, the CVC term in eq. (75) vanishes at T = if we replace 7k+ q 
with 7 imp ^ 0. However, the Born approximation is applicable only when the impurity 
potential is weak. In fact, we will show that a "strong" local impurity potential in 
under-doped HTSCs induces drastic and widespread changes in the electron-electron 
correlation around the impurity site. Thus, the impurity effect in HTSCs depends on 
the strength of the impurity potential. Hereafter, we study the nontrivial impurity 
effects in nearly AF Fermi liquids. 
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7.1. Hall coefficient in the presence of weak local impurities 

First, we study the effect of "weak local impurities" on the transport phenomena 
within the Born approximation, where the quasiparticle damping rate due to impurity 
scattering is given by 

7 imp = n imp I 2 J2 ImG k (-a) = 7m imp / 2 iV(0), (131) 

k 

where n- irav is the impurity density, / is the impurity potential, and iV(0) is the DOS at 
the chemical potential. Then, the total quasiparticle damping rate 7k is given by 

7k = 7k + 7 imp , (132) 

where 7k represents the quasiparticle damping rate due to inelastic scattering, which is 
given by eq. (68) in the Fermi liquid theory or by eq. (24) in the FLEX approximation. 
(Here, we ignored the self-energy correction given by the cross terms between U and I. 
This simplification is not allowed for a large J, as we will show in §7.2.) Also, the CVC 
due to local impurities vanishes identically within the Born approximation. Therefore, 
in the case of 7 imp 7^ 0, the Bethe-Salpeter equation (75) is changed to become 

~4 + «kJk*, (133) 

where a^ = a^ ■ (7k/7k)- Thus, an approximate solution of eq. (133) is 

1 

Jk=z zo (^k + «kWk*)- (I 34 ) 

l-a k 

In the absence of impurities, Jk exhibits singular k-dependence because «k ~ (1 — c /£af) 
approaches one. In the case of 7k ~ 7 imp ) on the other hand, cx^ <ti a^ < I and therefore 
J\s. ~ ^k- As a result, the CVC is strongly suppressed by high density "weak local 
impurities" , when the elastic scattering is comparable to the inelastic scattering due to 
AF fluctuations. 

In Ref. [34] , we have pointed out that the enhancement of i?H due to the CVC is 
easily suppressed by the weak impurities. Figure 27 shows the impurity effects on p and 
i?H obtained by the FLEX+CVC method for electron-doped systems. As expected, Rn 
is suppressed by a small amount of impurities, although the induced residual resistivity 
is small. (Note that p = 1 corresponds to 250 /zffcm.) As a result, i?H becomes 
positive at low temperatures in the presence of impurities, which is consistent with the 
experimental behavior of Rn in PCCO for S = 0.16 ~ 0.18 [8]. Finally, we consider the 
case where the impurity potential is widespread (nonlocal). In this case, the impurity 
potential /(q) is momentum- dependent, and the kernel of the Bethe-Salpeter equation 
in eq. (133), (3f/ 2 /2)(7rT) 2 ImXq(0), is replaced with n imp / 2 (q) + (3f/ 2 /2)(7rT) 2 Im^(0). 
When J(q) is large only for q ~ (forward scattering), both the residual resistivity and 
the reduction in Rn due to the impurities will be small. In HTSCs, it is considered that 
impurities outside of the Cu0 2 plane causes the forward impurity scattering [202] . 
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Figure 27. Theoretically obtained T-dependence of (i) p and (ii) i?H for NCCO in the 
presence of elastic scattering due to weak local impurities. In the case of 7i mp ^ 0, i?jj 
becomes positive at low temperatures since the CVC is reduced by elastic scattering. 
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Figure 28. Obtained T-dependence of cotfe = p/Rb. for LSCO for 0.01 < T < 0.08. 
cot#H oc T ° + c is realized well, and c is proportional to 7; mp . 



Figure 28 shows the impurity effects on cot#n = p/Rr obtained by the FLEX+CVC 
method for LSCO. The relationship cot#H oc T 1 - 6 + c holds well, and c oc 7i mp . Such 
a parallel shift of cot#n by impurity doping is observed in various hole-doped HTSCs 
[169]. 



7.2. Effect of strong local impurities near AF QCP 

Now, we discuss the effect of "strong local impurities" in HTSCs. According to a recent 
LDA study [203], a Zn atom introduced in the Cu02 plane of HTSCs induces a large 
positive potential (less than lOeV), and the potential radius is only ~ 1A. In HTSCs, 
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however, Zn doping causes a nontrivial widespread change of the electronic states. In 
Zn-doped YBCO compounds, site-selective 89 Y NMR measurements revealed that both 
local spin susceptibility [204, 205] and staggered susceptibility [206] are prominently 
enhanced around the Zn site, within a radius of the AF correlation length £af- The 
same result was obtained by the 7 Li Knight shift measurement in Li-doped YBCO 
compounds [207], and by the 63 Cu NMR measurement in Zn-doped YBCO compounds 
[208]. Moreover, a small concentration of Zn induces a huge residual resistivity, which 
is significantly greater than the s-wave unitary scattering limit [156]. These nontrivial 
impurity effects were frequently considered as evidence for the breakdown of the Fermi 
liquid state in under-doped HTSCs. 

Up to now, many theorists have studied this important issue. The single-impurity 
problem in a cluster t-J model has been studied by using the exact diagonalization 
method [209, 210]. When the number of holes is two, both AF correlations and electron 
density increase near the impurity site. This problem was also studied by using the 
extended Gutzwiller approximation [211]. Although these methods of calculation are 
founded, the upper limit of the cluster size is rather small. Moreover, these studies are 
restricted to T = 0. 

In this section, we study a single-impurity problem in a large size square lattice 
(say 64 x 64) Hubbard model based on the nearly AF Fermi liquid theory [130]: 

H = 5Z *«'4t c j» + U 5Z n ^ n ^ + j Kt + noi), (135) 

i,j,cr i 

where I is the local impurity potential at site % — 0. It is a difficult problem since 
we have to consider two different types of strong interactions (U and I) on the same 
footing. Moreover, the absence of translational symmetry severely complicates the 
numerical analysis. To overcome these difficulties, we developed the GV 1 method - 
a powerful method for calculating the electronic states in real space in the presence of 
impurities [130]. The GV 1 method is applicable for finite temperatures since the thermal 
fluctuation effect is taken into account appropriately. Based on the GV 1 method, we 
successfully explain the nontrivial impurity effects in HTSCs in a unified way without 
assuming any exotic mechanisms. 

In the GV 1 method, the real-space spin susceptibility x Is (ui) is given by 



X Is = T1 I (1- UT1 1 ) , (136) 

n 7 (ri, iv,-; ui) = ~T^2 G \ r ^ r i5 e n + ^i) Gl ( r j, r i', e n)- ( 13 7) 

Here, G 1 is given by solving the following Dyson equation; 

G\e n ) = G°(e n ) + G°(e n )/G J (e n ), (138) 

where (I)ij = Io~i,oo~j,o, an d G° is the real-space Green function in the FLEX 
approximation without the impurity potential /. The solution of eq. (138) is 

Gi^Gl + I^f, (139) 

r o,o 



'" ~ 1 ' 3 ' '1-JG5, 
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Figure 29. (i) The uniform susceptibility given by the GV 1 method in the presence of 
impurities (n lmp = 0.01) for the impurity potential / = oo. A single impurity induces a 
large Curie-term Ax", which is consistent with the experimental results, (ii) The local 
susceptibility % s (r, r) and staggered one x s ( r i r + (1, 0)) given by the GV 1 method. 
The inset describes the state in which both local and staggered susceptibilities are 
enhanced around the impurity site. 



where the position of the impurity potential is i = 0. 

In Fig. 29, we show the numerical results for spin susceptibilities, x Is , for I = oo. 
Surprisingly, a nonmagnetic impurity induces a huge Curie-like component in the 



uniform susceptibility; Ax 



n 



imp 



/Xg ff /3T. This result explains the long-standing 



experimental problem [212, 213, 214]. The obtained value of /^ e ff is 0.74/xb, which is close 
to the experimental value /i e s ~ 1/xb in YBa 2 Cu30 6 . 6 6 (T c w 60K) [212] and in LSCO 
(5 = 0.1) [214]. We also emphasize that both the local and staggered susceptibilities 
get enhanced around the impurity site within a radius of about 3 (~ £af) at T = 0.02. 
Here, we discuss the physical reason for this drastic impurity effect. In the FLEX 
approximation, the AF order (in the mean-field level) is suppressed by the self-energy 
correction (see §2), which represents the destruction of the long-range order due to 
thermal and quantum fluctuations. In the FLEX approximation, the reduction in the 
DOS due to the large quasiparticle damping rate 7k = ImE^ renormalizes the spin 
susceptibility. Around the impurity site, the self-energy effect due to electron-electron 
correlation is expected to be smaller. In fact, quantum fluctuation is reduced near 
the vacant site in the s = 1/2 Heisenberg model [215]. For this reason, in the GV 1 
method, the AF correlations are prominently enhanced around the impurity site due to 
the reduction of thermal and quantum fluctuations. The examples of the cross terms 
between U and I for susceptibility are shown in Fig. 30 (i). The same physics occurs 
in quantum spin systems with vacancies; AF spin correlations are enhanced around a 
vacancy, due to the reduction in quantum fluctuations [215, 216, 217, 218, 219]. 
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Figure 30. (i) Examples of diagrams for \ in the presence of a single-impurity 
potential I in the Hubbard model. They are cross terms between / and U. The cross 
represents /, and the large circle with the shadow represents the self-energy given by 
the FLEX approximation for 1 = 0. (ii) The local DOS given by the GV 1 method: it 
decreases around the impurity site within the radius of about 3, due to the nonlocal 
"effective complex impurity potential", <5E(r,r'; e = 0). 



Prior to the GV 1 method, the single impurity problem in the Hubbard model had 
been studied using the RPA [220, 221, 222], assuming a nonlocal impurity potential; Vo 
for the onsite and V\ for the site adjacent to the impurity atom. The obtained result 
strongly depends on the value of V\. When V\ = 0, the enhancement of susceptibility is 
tiny or absent although V\ ~ for the Zn impurity atom in the CuCVplane according 
to the LDA study [203]. On the other hand, in the GV 1 method, the considerable 
enhancement of the susceptibility is realized even when V± — (local impurity potential 
case), since xq is renormalized by the self-energy correction in the FLEX approximation: 
Then, \Q is easily enhanced when the self-energy effect is reduced by introducing an 
impurity. Therefore, we have to use the GV 1 method for investigating the impurity 
problem in HTSCs. 

Next, we discuss the transport phenomena in the presence of dilute impurities 
according to the GV 1 method [130]. The enhanced susceptibilities due to the impurity 
induce the "additional self-energy correction <5X = X 7 — S /=0 " around the impurity 
site. If the area of 5Y, ^ is large, a large residual resistivity will be induced by the 
non-s-wave scattering channels. In the GV 1 method, <5E is given by 



5E(rj, r,-; e„)=TV G(r h r.,-; io x + t n )V\r u r^; 



i 

E°f 



UJl 



^ j) ^-n}i 



(140) 
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(141) 

where S° is the self-energy given by the FLEX approximation without the impurity 
potential (/ = 0). The Green function G(rj,rj-;e n ) is given by solving the following 
Dyson equation in real space: 

G(e n ) = G\e n ) + G\e n )5t{t n )G{e n ). (142) 

In the GV 1 method, we solve eqs. (140) and (142) self-consistently The local DOS 
given by the GV 1 method in real space -- p(r,uj) = ImG(r, uo — i8)/n -- is shown in 
Fig. 30 (ii). The local DOS decreases around the impurity site within the radius of 
approximately 3a (a is the lattice spacing), since the quasiparticle lifetime is very short 
due to the large Imo"S. We verified that the radius of Imo"£ increases as the filling 
number n approaches unity, in proportion to £af- As explained in Ref. [130], we should 
not solve V 1 in eq. (141) self-consistently, since the feedback effect on V 1 introduced 
by iteration is canceled by the vertex correction for V 1 that is absent in the FLEX 
approximation. 

To derive the resistivity in the presence of impurities, we have to obtain the t- 
matrix, t(e), which is defined as G = G° + G°tG°. The expression of t(e) in the case of 
^imp <C 1 is derived in Ref. [130]. Using the t-matrix, the quasiparticle damping rate 
due to the impurity is given by [40, 41] 

"imp V r 
~Jj2 ' 1^ 



7r P (e) = ^VlmT,(e-ioy kr <, (143) 



i 

where T}(e) = ^2 m t m , m +i(c) ■ n- mip is the density of the impurities. In the case of 
^imp 7^ 0, the total quasiparticle damping rate is 7k(e) = 7k( e ) + 7k mP ( e )> where 



cold' 



7k (e) = ImE£(e — iS). The resistivity for n- imx> ^ is approximately given by p oc 7 

Figure 31 (i) shows k-dependences of 7™ p (n imp = 0.05) and 7^. We emphasize that 
7™ p exhibits strong k-dependence that is similar to the k-dependence of 7^. This result 
suggests that the effective impurity potential SH is nonlocal. As a result, the structure 
of the "hot spot" and the "cold spot" , is not smeared out by the strong non-magnetic 
impurities, although it will be smeared out by the weak local impurities. This finding 
strongly suggests that the enhancement of the Hall coefficient near the AF QCP, which 
is induced by the strong CVC around the cold spot [33, 32], does not decrease due to the 
strong impurities. Moreover, the enhancement of xaf due to n imp will also enhance i?H- 
Experimentally, i?n for under-doped YBCO and Bi 2 Sr2_ x La x Cu06+«s are approximately 
independent of the doping of Zn and other non-magnetic impurities [9, 10, 5], whereas 
i?H for LSCO decreases with Zn doping [11, 12]. 

Figure 31 (ii) shows the resistivity p for YBCO with n imp =0.02 for / = 0,4,6, 
and 00. The obtained impurity effect is the most prominent when 1 = 6. As T 
decreases, nonmagnetic impurities cause a "Kondo-like upturn" of p below T x , reflecting 
an extremely short quasiparticle lifetime near the impurities. In other words, the 
considerable residual resistivity is caused by the large scattering cross section of the 
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Figure 31. (i) 7™ p for n imp = 0.05 and 7° for YBCO and LSCO at T = 0.02. 
k-dependent 7J c mp obtained using the GV 1 method will be verified by the ARPES 
measurements for Zn-doped compounds, (ii) p in YBCO (n = 0.9) for n; mp = 0.02 
given by the GV 1 method. At lower temperatures, an insulating behavior is observed 
in the close vicinity of the AF QCP. In HTSCs, T = 0.1 and p = 1 correspond to 400K 
and 250/i$7-cm, respectively, (iii) The filling dependence of Ap = p lmp — po per ni mp 
at T = 0.02 given by the GV 1 method. Note that Ap = (4n imp /n) [h/e 2 ] in a 2D 
electron-gas model due to strong local impurities. 



"effective impurity potential", 511, as shown in Fig. 30 (ii). Note that T x increases 
with n- imp . This amazing result strongly suggests that the insulating behavior of p 
observed in under-doped LSCO [223, 224] and NCCO [225] is caused by the residual 
disorder in the Cu02 plane. As shown in Fig. 31 (iii), the parallel shift of the resistivity 
at finite temperatures due to impurities (Ap) grows drastically in the under-doped 
region; Ap exceeds the s-wave unitary scattering limit in the 2D electron gas model, 
Ap = (ft/e 2 )(4ni mP/ M). This result effectively explains experimental carrier-doping 
dependence of Ap [156]. 

The enhancement of the residual resistivity is also observed near the AF QCP in 
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heavy-fermion systems such as CeAl 3 [226] and CeCusAu [227], and in the organic 
superconductor /t-(BEDT-TTF) 4 Hg 2 .89Br 8 [239]. In these compounds, Ap quickly 
decreases with pressure, as the distance from the AF QCP increases. In k-(BEDT- 
TTF) 4 Hg2.89Br 8 , Ap under 2 GPa is six times smaller than the value at 0.5 GPa. 
Such a large change in Ap is difficult to be explained by the pressure dependence of 
the DOS. In fact, according to the t-matrix approximation, 7 imp = (irI 2 N(0)/2)/(l + 
(ttIN(0)/2) 2 ). The residual resistivity in 2D free dispersion model (e£ = k 2 /2m) is 
Ap = 27riV(0)7i m p/e 2 n. Therefore, Ap is given by 

H 4(«IN(0)/2y n imp 
1 e 2 1 + (7r/iV(0)/2) 2 n ' K } 

in 2D free dispersion model. Note that the renormalization factor z does not appear 
in the expression of Ap. In case of IN(0) ^> 1, eq. (144) gives the s-wave unitary 
scattering value; Ap = {%/ e 2 )4(ni mp /n) . In case of weak impurity scattering where 
IN(0) <C 1 (Born limit), we obtain the relation Ap oc J 2 iV 2 (0), which will decrease 
under pressure since N(0) oc 1/Wband- However, pressure dependence of Ap in k- 
(BEDT-TTF) 4 Hg 2 .89Br 8 seems too strong to be explained by the Born approximation. 
We comment that the increment in Ap due to charge fluctuations is discussed in some 
heavy-fermion systems [228, 229]. 

In summary, the present study revealed that a single impurity strongly influence 
the electronic states in a wide area around the impurity site near the AF QCP. Using 
the GV 1 method, the characteristic impurity effects in under-doped HTSCs are well 
explained in a unified way, without assuming any exotic non- Fermi-liquid ground states. 
We successfully explain the nontrivial impurity effects in HTSCs in a unified way in 
terms of a spin fluctuation theory, which strongly suggests that the ground state of 
HTSCs is a Fermi liquid. We expect that the novel impurity effects in other metals near 
the AF QCP, such as heavy-fermion systems and organic metals, will be explained by 
the GV 1 method. The validity of the GV 1 method is verified in Ref. [130] based on the 
microscopic Fermi liquid theory. 

Recently, we extended the GV^FLEX method to be applicable even in below T c , 
and analyzed the Hubbard cluster model given by eq. (135). We found that both local 
and staggered susceptibilities are enhanced around the impurity site even in the SC 
state [230]. This result is consistent with 17 NMR measurement for YBCO with Zn 
impurities [231]. Interestingly, it is also reproduced by analyzing the Hubbard model 
together with the rf-wave pairing interaction using the mean- field approximation [232]. 
In the normal state, on the other hand, the enhancement of the local and staggered 
susceptibilities cannot be explained qualitatively within the mean-field approximation 
or RPA [220, 221, 222]. It is an important future issue to explain the sensitive impurity 
effects on the SC state in HTSCs [199] based on the repulsive Hubbard model, using the 
GV 7 -FLEX method. 
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Figure 32. Temperature dependence of i?H in (a) CeRhIn5 for P > P c 
and in (b) CeCoIn 5 for P > 0. [Rcf. [20]] 



2.01 GPa, 



8. Anomalous transport behaviors in CeMIn 5 (M=Co or Rh) and 
k-(BEDT-TTF) 

Thus far, we showed that the CVC due to strong AF+SC fluctuations induces various 
anomalous transport phenomena in HTSCs. To validate this idea, we have to study the 
various nearly AF systems other than HTSCs. In general, the electronic structure of 
heavy-fermion systems and organic metals are very sensitive to the pressure. Therefore, 
the distance from the AF QCP can be easily changed by applying the pressure, without 
introducing disorders in the compounds. This is a great advantage with respect to 
investigating the intrinsic electronic states near the AF QCP, free from the disorder 
effects. A useful theoretical review for heavy-fermion systems near the AF QCP is 
given in Ref . [233] . Recently, detailed measurements of the transport phenomena under 
pressure have been performed in the heavy-fermion superconductor CeMIns (M=Co 
or Rh) and in the organic superconductor K-(BEDT-TTF). They exhibit striking non- 
Fermi-liquid-like behaviors as observed in HTSCs — eqs. (l)-(4). Hereafter, we explain 
the experimental and theoretical studies on the transport phenomena in these systems. 



8.1. CeMIn 5 (M=Co or Rh) 

CeMIns is a quasi 2D heavy-fermion superconductor with T c = 2.3K. According to the 
angle-resolved measurements of thermal conductivity [102] and specific heat [105], the 
symmetry of the SC gap is rf-wave. Figure 32 shows the temperature dependence of 
i?H in CeCoIns and CeRhhis in the limit of H z — 0. In CeRhlns, the critical pressure 
under which the AF QCP is realized is P c « 2.01GPa. At P = P c , J R H (2.3K)/ J R H (300K) 
reaches 50, whereas the magnitude of -Rh rapidly decreases as the pressure is increased. 
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Figure 33. Ap/po plotted as a function of tan 2 Oh for CeRhlns. [Rcf. 



A similar pressure dependence of -Rr- is observed in CeCoIns, where P c is slightly below 
the ambient pressure. In both these compounds, Ru is inversely proportional to T 
at higher temperatures. Moreover, as shown in Fig. 33, the modified Kohler's rule 
given in eq. (106) is well satisfied in CeRhIn 5 [20] for over four orders of magnitude. 
Furthermore, this is well satisfied in CeCoIn 5 [19, 20]. Therefore, both Rn and Ap/p 
in CeMIns shows anomalous behaviors which are similar to HTSCs. This experimental 
fact strongly suggests that the CVC due to the strong AF fluctuations is the origin of 
the anomalous transport phenomena. 

Now, we discuss the magnetic field dependence of the transport coefficients. 
Interestingly, Rn = dpu(H z )/dH z in CeMIn 5 near P c is easily suppressed only by a 
small magnetic field, as shown in the inset of Fig. 32. At the same time, (Ap/po)H~ 2 
is also significantly suppressed. Therefore, the relationships a xy oc H z and Ao~ xx oc H 2 
are satisfied only below ~ 0.1 Tesla near the QCP [19, 20]. These behaviors cannot be 
attributed to the orbital effect (i.e., the cyclotron motion of conduction electrons) since 
u*t* = (eH z /m*c)T* > 1 is satisfied only when H z > H c2 = 5 Tesla and T< IK 
[234]. The condition uj*t* <C 1 is also recognized from the relationship Ap/p i$ 0.1 for 
T > 2K and H < 3 Tesla, as shown in Fig. 33. 

When uj*t* <^ 1 is satisfied, we can safely expand Nakano-Kubo formula given in 
eq. (42) with respect to the vector potential, as we did in §4 - 7. Then, the following 
relationships derived in §4 

p H oc H z e AF , (145) 

Ap/p oc Hl&v/pl (146) 

are expected to be valid for CeMIn 5 . Here, the factors £^ F and £^f in eqs. (145) 
and (146), respectively, come from the CVCs. Therefore, experimental striking non- 
linear behaviors of -Rr- and Ap/po with respect to H z should originate from the field- 
dependence of Xq(0) oc £af : The suppression of £af due to the magnetic field results 
in reducing both eqs. (145) and (146). Since the correlation length is sensitive to the 
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outer parameters in the vicinity of the QCP [235] , the anomalous sensitivity of pn and 
Ap/pcj to the magnetic field in CeMIn 5 originates from the field-dependence of the CVC. 
Recently, pu and Ap/p in PCCO were measured in magnetic field up to 60 T, and it is 
found that both of them shows striking non- linear behaviors with respect to H z [236]. 
Their behaviors will be explained by the field-dependence of the CVC. 

Even if the field dependence of CVC is prominent, the modified Kohler's rule in 
eq. (106) should be satisfied for a wide range of the magnetic field strength, since both 
Ap/po and cot 2 #H are proportional to £afP -2 if ^e preset theory of CVC is correct. In 
fact, the modified Kohler's rule is well satisfied for T = 2.5 — 30K and < H z < 3 
Tesla in CeRhhis as shown in Fig. 33, regardless of the fact that the conventional 
Kohler's rule is violated only for 0.1 Tesla. This fact strongly suggests that both a xy 
and Aa xx are enhanced by the same origin, namely, the CVC due to the AF fluctuations. 
Therefore, the anomalous transport phenomena in CeCoIns are consistently described 
by the theory of CVC in nearly AF Fermi liquids. 

The Nernst signal v in CeCoIns is also very anomalous [107]. Below 20K, v starts 
to increase approximately in proportion to T" 1 , exhibiting anomalously large values 
[y ~ lpV/KT) below 4K. This behavior is very similar to that of v in electron-doped 
HTSC, whose T-dependence and the magnitude are well reproduced by the CVC due to 
strong AF fluctuations. It reaches ~ O.lpV/KT in optimally-doped NCCO. Considering 
that v is proportional to r oc p^ 1 , the experimental relations & , c e coin 6 (4K)/* / Ncco(iooK) ~ 
l(pV/KT)/0.1(pV/KT) = 10 and PccComWPNccoaooK) ~ 5(pfi)/50(pfi) = 1/10 are 
naturally understood as the quasiparticle transport phenomena. Thus, the giant Nernst 
signal in CeCoIn 5 is expected to be caused by the CVC due to strong AF fluctuations. 

Although the spin susceptibility x s (q) i n CeMIn 5 has a quasi 2D structure, the 
anisotropy of resistivity is only two (p c /p a b ~ 2) at low temperatures. In the dynamical- 
mean-field-theory (DMFT) [109, 110], which is believed to effectively describe the 
electronic properties in 3D systems, the CVC vanishes identically. Therefore, it is 
highly desired to confirm whether the CVC can be really significant in 3D systems 
or not. Recently, we performed the numerical study for 3D Hubbard model with the 
conduction electron spectrum is e^ D = e£ + 2t z cos k z , where e£ is given in eq. (12). The 
obtained p, Ru and v in quasi 2D Hubbard model are shown in Fig. 34 (i)-(iii). Since 
the unit cell lengths in CeMIns are a a = a^ = 4.6 Aand a c = 7.6 A, p = 1 and Ru = 1 
correspond to 300 pf2cm and 1.0 x 10 -9 m 3 /C, respectively. We find that both i?n and 
v are strikingly enhanced due to the CVC even if t z /to ~ 1. In CeCoIns, R^ starts to 
increase below ~ 40K, which corresponds to T = 0,08 in Fig. 34 (ii). At T = 0.02 
(~8K for CeMIn 5 ), Ru ~ 5 for t z /t = 0.8. Then, the extrapolated value of Rn at 2K 
is 20. 

As shown in Fig. 34 (iii), the obtained v is drastically enhanced by the CVC, which 
is consistent with experiments. To derive the correct absolute value of z/, we have to 
take account of the fact that v is proportional to the mass-enhancement factor 1/z. 
According to the de Haas-van Alphen measurement, 1/z = m*/m^, an d ~ 50 in CeCoIns. 
Since the mass-enhancement factor in the present FLEX approximation is 1/^flex ~ 3, 
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we describe the experimental mass-enhancement factor as 1/z = 1/(^flex-2*)i where 
1/z* = 50/3 is the mass enhancement factor which cannot be described by the FLEX 
approximation in the Hubbard model. [Even in the FLEX approximation, relatively 
large 1/-2flex ~ 10 is obtained in the periodic Anderson model, which is an effective 
model for heavy fermion systems [237]. We will comment on this fact in §9.2.] We present 
v — ^flex/- 2 * in Fig. 34 (iii), where z^flex is given by the FLEX approximation, by 
using the relation h^a\/h = 28 nV/KT. Recently, we extended the FLEX approximation 
to reproduce appropriate results under a finite magnetic field, and calculated the 
field dependence of Ru and v. It was suggested that both these quantities were 
rapidly suppressed by the magnetic field near AF QCP, reflecting the reduction of AF 
fluctuations. The obtained result is in good agreement with the experimental results. 

In CeCoIns, Rn exhibits a peak at T^ ~ 4K at ambient pressure, and T^ increases 
with pressure [see Fig. 32]. This is different from the pseudo-gap behavior in under- 
doped HTSCs since T^ H in CeCoIn 5 increases as the system goes away from the AF-QCP. 
In particular, i?H of CeRhlns at P = P c maintains its increasing trend just above T c . 
These behaviors can be understood as the effect of the weak (local) residual disorders, 
as we have discussed in §7.1: as shown in Fig. 27 (ii), \Ru — -Rh TA | starts to decrease 
at lower temperatures when 7 imp > 0. The peak temperature T^ increases with 7i mp . 
In CeCoIn 5 , the resistivity at T* Rn is p(T = T£j ~ 6//fi cm [20] for P = ~ 2.5 GPa. 
On the other hand, in CeRhlns at P = 2 GPa, p(T > T c ) f« 10/xfi cm because of the 
large inelastic scattering. Since 7 imp in CeCoIn 5 and that in CeRhIn 5 are expected to 
be similar, <Sk = «k • 7k/ (7k + 7im P ) in eq. (134) will be smaller in CeCoIn 5 . That is, 
reduction of CVC due to impurities is more prominent in CeCoIn 5 . Therefore, we can 
explain the different behaviors of Ru in CeRhlns and CeCoIn5 at low temperatures as 
the effect of residual disorders. 

One may ascribe the temperature dependence of i?H i n CeM In 5 to the multiband 
effect. For example, a sign change in i?n can occur if a hole-like Fermi surface and 
an electron-like Fermi surface coexist and their mean free paths have different In- 
dependences. However, it is very difficult to explain the relation \Rn\ ^> 1/ne in a 
multiband model based on the RTA. In fact, in order to explain the pressure dependence 
of i?H in CeCoIns, one has to assume that a small Fermi surface governs the transport 
phenomena at OGPa (near AF QCP), whereas a large Fermi surface should consequently 
govern the transport phenomena under 2.5GPa. The same drastic change in the 
electronic states should occur by applying a magnetic field H ~ ITesla. This unnatural 
assumption is not true since the other transport coefficients cannot be explained at all, 
as discussed in Ref. [20] in detail. In particular, the elegant modified Kohler's rule 
plot for CeRhIn 5 shown in Fig. 33 for over four orders of magnitude is strong evidence 
that the CVC is the origin of the anomalous transport phenomena, and the anomalous 
transport phenomena are mainly caused by a single large Fermi surface with heavy 
quasiparticles. Since the enhancements of both i?n and Ap/p originate from a small 
portion of the Fermi surface (i.e., the cold spot) as shown in Fig. 19 (ii), the multiband 
effect is unimportant in CeMIn 5 . Therefore, the modified Kohler's rule will be realized 
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Figure 34. (i) p, (ii) i? H and (iii) v for the 3D Hubbard model (t z = 0.4, 0.8) 
given by the FLEX+CVC method. Broken lines represent the numerical results for 
U = 9.4 (t z = 0), U = 6.2 (t z = 0.4) and U = 5.4 (t z = 0.8). In these parameters, the 
Stoner factors as is 0.995 at T — 0.02. Full lines represents the numerical results for Us 
with which a s = 0.998 is satisfied at T = 0.02. (iv) Fermi surfaces for t z = 0, 0.4, 0.8. 
[Ref. [108]] 



near the AF QCP even in multiband systems like CeMIn 5 . 

In usual heavy-fermion compounds, i?H exhibits a Curie-like behavior above the 
coherent temperature T due to the large anomalous Hall effect, which is caused by the 
angular momenta of the /-orbitals as we will discuss in §9.4. For T <C T , -Rr- due to 
the AHE is proportional to p 2 . In heavy-fermion compounds near the AF QCP, on the 
other hand, both the ordinary Hall coefficient K^ and the anomalous Hall coefficient 
Rft UE can show large temperature dependences. Paschen et al. extracted i?g from the 
experimental Hall coefficient in YbRh 2 Si2, and discussed the critical behavior of K^ near 
AF QCP, which is realized under the magnetic field [22]. In contrast, i?n in CeMIns 
(M=Co,Rh) is almost constant above 50 K as shown in Fig. 32. Therefore, anomalous 
Hall effect in CeMIn 5 is very small, and therefore _R^ F ~ R^ as explained in detail in 
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Figure 35. Temperature dependence of (i) i?H and (ii) cot6>H for k-(BEDT- 
TTF) 2 Cu(NCS) 2 under pressure. [Ref. [17]] 



Ref. [19, 20]. This fact is a great advantage to study the anomalous T-dependence of 
the ordinary Hall effect in CeMIns. 

In CeCoIn 5 , novel kinds of critical behaviors are observed near H c2 ~ 5 Tesla 
for T ^C 1 K, where prominent increment in the effective mass was observed. This 
phenomenon is referred to as a "field-induced QCP". One possible origin will be the 
field-induced SDW state that is hidden in the SC state. We will discuss this important 
future problem in §9.2 in more detail. 

8.2. k-(BEDT-TTF) 

The measurements of i?H under pressure have been intensively performed in k- 
(BEDT-TTF) 2 X. Figure 35 shows the temperature dependence of Ru and cot# H for 
X=Cu(NCS)2 under homogeneous pressure. Here, the AF correlations are the strongest 
at the lowest experimental pressure (0.19GPa). As the pressure is increased, the AF 
fluctuations get reduced and a conventional Fermi liquid state is realized. At the same 
time, i?H is reduced and exhibits a constant value. The observed Rn is independent of 
pressure at room temperature. Therefore, the origin of the increment in i?n below 100 K 
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cannot be related to deformation of the Fermi surface under pressure. Note that similar 
increment in i?n is observed in X=Cu[N(CN) 2 ]Cl, which has a single elliptical Fermi 
surface [92]. Because the observed T dependences of Ru and cot6>H are very similar 
to the observations involving HTSCs and in CeMIn 5 , the CVC is expected to play a 
significant role in k-(BEDT-TTF) 2 X. 

As discussed in §2.3, the resistivity in 2D systems in the presence of AF fluctuations 
is [147], 

P oc T 2 e AF - (147) 

This relationship is reliable when cj s f > T, which is satisfied in optimally- or over- 
doped HTSCs. The T-dependence of o; s f is given in eq. (10). We consider eq. (147) is 
realized in X=Cu(NCS)2 since the AF fluctuations are not so prominent. Since -Rr- is 
proportional to £j[ F as shown in eq. (97), we obtain 

cot# H = p/Rh ocT 2 . (148) 

Figure 35 (ii) shows that eq. (148) is well satisfied below 80K, except at P = 0.19 
GPa. [The thermal contraction of the sample might modify this relationship at the 
lowest pressure.] The success of the scaling relationships (97), (147) and (148) is strong 
evidence that the enhancement of Rn in k-(BEDT-TTF) 2 X is caused by the CVC due 
to AF fluctuations. 

In many /t-(BEDT-TTF)2X compounds, the phase transition between the Mott 
insulating phase and the metallic (SC) phase induced by pressure is weak first order. 
Therefore, the AF fluctuations are not so strong even in the vicinity of the AF insulating 
phase. As a result, the enhancement of Ru in /t-(BEDT-TTF)2X is much smaller than 
that in CeMIn 5 . On the contrary, carrier doped (11% doping) K-type superconductor 
K-(BEDT-TTF)4Hg2.89Brg (T c = 4 K at ambient pressure) exhibits very strong AF 
fluctuations, and therefore the critical pressures is expected to be slightly below kbar. 
According to NMR measurement, \jT{T [oc £^f] increases with decreasing temperature 
above ~ 10 K, in proportion to (T + @) _1 with = 13 K [238]. The P-T phase 
diagram and transport properties of /T>(BEDT-TTF) 4 Hg 2 .89Br 8 are carefully measured 
by Taniguchi et al. [239]: It was found that T c shows two-peak structure under pressure. 
Also, the power n in p = p + AT n increases from 1 to 2 with increasing pressure, which 
indicates that the AF fluctuations are suppressed under the pressure. 

Taniguchi et al. also measured the Ru in K-(BEDT-TTF)4Hg2.89Br8 under 
pressures, and found that Ru exhibits a Curie- Weiss temperature dependence; Ru oc 
(T — @rh) _1 where Grh < [18]. At 0.19GPa where AF fluctuations are strong, 
/?h(10K)/_R h (300K) reaches 10, whereas the enhancement of Ru is totally suppressed 
by 1 GPa. The behavior of Ru in K-(BEDT-TTF)4Hg2.89Br8 is very similar to the 
observation in CeMIn 5 . 6 RH decreases with pressure, and its extrapolated value to 
GPa (—13 K) coincides with the Weiss temperature of \jT{T at ambient pressure. This 
is strong evidence that the Hall coefficient in K-(BEDT-TTF)4Hg2.89Br8 is proportional 
to £^ F , due to the CVC induced by AF fluctuations. It is noteworthy that the residual 
resistivity Ap in /t-(BEDT-TTF) 4 Hg2.89Br 8 drastically decreases with pressure [239]. As 
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Figure 36. (i) U-T phase diagram for the anisotropic triangular lattice Hubbard 
model [t\ — 0.3 ~ 0.7] given by the FLEX approximation. The large (small) dots 
represent the d-wave T c (T N ). (ii) R R given by the FLEX+CVC method. [Refs. [93] 
and 



we have discussed in §7.2, the large residual resistivity near GPa is expected to be 
given by the enlarged effective impurity potential due to the many-body effect. 

Note that Ru in K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br changes to be negative under high 
pressures and low temperatures, when the electron-electron correlation becomes very 
weak. The origin of this experimental fact is considered to be the lattice deformation 
with the long-period structure inherent in this compound [13]. 

Here, we analyzed the anisotropic triangular lattice Hubbard model at half-filling, 
which is a theoretical effective model for /t-(BEDT-TTF) 2 X. Figure 36 (i) shows the 
numerical results for the U-T phase diagram given by the FLEX approximation. Here, 
to (> 0) and t\ (> 0) are the nearest neighbor integrals of the triangular lattice; to is 
for two of three axes and t\ is for the remaining axis, respectively. In many k-(BEDT- 
TTF) 2 X compounds, ti/to ~ 0.7. Then, the corresponding Fermi surface is hole-like 
and ellipsoidal around the T-point. Hereafter, we substitute t = 1 ~ 600 K. Then, 
the AF ordered state appears at Tn < 0.04 ~24K when U is as large as the bandwidth 
(~ 10). As U decreases, Tn decreases and the d^-^-wave SC appears next to the AF 
phase. The obtained U-T phase diagram well explains the experimental P-T phase 
diagram (Kanoda diagram [91]), since the increment in P corresponds to the reduction 
in £//Hband- Figure 36 (ii) also shows the T-dependence of Ru for ti/t = 0.7. We 
see that Ru increases at low temperatures due to the CVC and it becomes smaller 
as U decreases. At half-filling, obtained enhancement of i?n is n °t so large since the 
AF fluctuations in this model are rather weak due to the strong geometric frustration 
(ti/to ~ 0.7). Therefore, a consistent understanding of the P-T phase diagram as well 
as the P-dependence of Ru are achieved by using the FLEX+CVC method. 
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9. Discussions 

9.1. Summary of the present study 

In the present study, we investigated the mechanism of non-Fermi-liquid-like transport 
phenomena in strongly correlated Fermi liquids in the presence of strong magnetic 
fluctuations. This problem was first realized in the study of HTSC, which has been one 
of the central issues in HTSC. For example, the relaxation time approximation (RTA) 
for the highly anisotropic t^ model cannot reproduce relationships (l)-(4) at the same 
time, even for a narrow range of temperatures. Recent experiments have been revealed 
that these anomalous transport properties are not a special phenomena in HTSC, but a 
universal phenomena in Fermi liquids near AF QCP. We discussed two typical examples 
- CeMIn 5 (M=Co or Rh) and /t-(BEDT-TTF) 2 X. The ratio R H (2.5K)/R H (300K) in 
CeRhIn 5 reaches ~ 50, which is much larger than the ratio in HTSCs. 

Here, we studied this long-standing problem by developing the microscopic Fermi 
liquid theory. In the RTA, the momentum and energy conservation laws in quasiparticles 
scattering are violated, as explained in §1 and §3. To overcome this defect, we take 
account of the CVC to satisfy the conservation laws. The total current Jk, which is 
responsible for the transport coefficients, is given by the summation of the quasiparticle 
velocity v^ and the CVC AJk. In interacting electron systems, an excited quasiparticle 
induces other particle-hole excitations by collisions. The CVC represents the current due 
to these particle-hole excitations. Since the CVC is caused by quasiparticle interactions, 
it can be significant in strongly correlated Fermi liquids. However, the effects of CVC in 
HTSCs have not been studied in detail until recently. Since the CVC is totally ignored 
in the RTA, the RTA frequently yields unphysical results in strongly correlated systems. 

In the presence of strong AF fluctuations, we find that (a) quasiparticle damping 
rate 7k = l/2rk becomes anisotropic, and the portion of the Fermi surface with small 
7k (cold spot) governs the transport phenomena [24, 25]. At the same time, (b) the 
total current Jk becomes highly anisotropic due to the prominent CVC [33]. (a) and (b) 
should occur simultaneously since both of them are caused by the same origin; highly 
anisotropic incoherent scattering due to AF fluctuations. That is, both (a) and (b) are 
induced by the strong backward scattering between k and k* (k — k* m Q), as shown in 
Fig. 4 (i). Mathematically, 7k = ImSk(— iS) and the CVC are closely connected by the 
Ward identity near AF QCP, that is, the same vertex function 7^ k , appears in eqs. (68) 
and (60). The facts (a) and (b) naturally explain the enhancement of Rn, Ap/p , and 
v in nearly AF metals in a unified way. As shown in Fig. 14 (ii), the large curvature 
of the effective Fermi surface around the cold spot gives rise to the enhancement of i?n- 
In §3.4, we proved the significance of the CVC beyond the one-loop (such as the FLEX) 
approximation. 

It is also very important to demonstrate to what extent the CVC can reproduce the 
aspects of anomalous transport properties based on a standard spin fluctuations theory. 
For this purpose, we study -Rh, Ap/po, S and v based on the FLEX (or FLEX+T- 
matrix) approximation by including CVCs. The obtained results semiquantitatively 
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reproduce the various experimental facts in slightly under-doped HTSCs. The present 
study strongly suggests that the striking deviation from the Fermi liquid behaviors (such 
as eqs. (l)-(J,.)) are ubiquitous in strongly correlated metals near the AF QCP, not 
specific to HTSCs. Note that the increment in i?n is also observed in non-SC metals 
near AF phase, such as V2-j / Os[240] and i^-xBirRx^Oy (R = Sm or Eu)[241]. 
The main results in the present study are as follows: 

(i) Hall coefficient: Due to the CVC, i?n oc £^ F oc T _1 . The large curvature of the 
effective Fermi surface around the cold spot is the origin of the enhancement of R-^. 
In particular, we can explain that the /?h for NCCO and PCCO are negative even 
above Tn ~ 150 K, since the curvature of the "effective Fermi surface" at point B 
in Fig. 14 is negative, although the curvature of the true Fermi surface is positive. 

(ii) Magnetoresistance: Ap/p oc £af/Po ^ ^~ 4 due t° the CVC, which means that 
the conventional Kohler's rule Ap/p oc p$ 2 is violated. On the other hand, the 
modified Kohler's rule Ap/po oc (Ru/po) 2 is realized since the same factor £^f 
appears in both sides. 

(iii) Thermoelectric power: For T > T*, \S\ increases as T decreases since the 
e-dependence of Tk(e) becomes prominent when the AF fluctuations are strong. 
Further, S > (S < 0) in many hole-doped (electron-doped) systems. 

(iv) Nernst coefficient: Due to the CVC, v gradually increases in NCCO below the 
room temperature. With regard to LSCO, a rapid increase in v in the pseudo-gap 
region is well explained by the FLEX+T-matrix approximation, due to the CVC 
induced by the strong AF+SC fluctuations. 

(v) AC Hall effect: The strong ^-dependence of Ru(uj) appears due to the frequency 
dependence of the CVC. Both T- and ^-dependences of Ru(w) are well reproduced 
by the FLEX+CVC method. 

(vi) Impurity effect: The strong CVC near the AF QCP is suppressed by the weak 
local impurities. For strong local impurities, the residual resistivity exceeds the 
s-wave unitary scattering value since the local and staggered susceptibilities are 
strongly enhanced around the impurity sites. 

Results (i)-(vi) are given by the same mechanism - - the hot/ cold spot structure 
and the singular ^-dependence of the total current J^ in nearly AF metals, which is 
shown in Fig. 14. We emphasize that the abovementioned results are not derived from 
a special defect of the FLEX approximation. They are also reproduced by another 
spin fluctuation theory. In fact, Kanki and Kontani [34] calculated J^ based on the 
Millis-Monien-Pines model in eq. (8), by assuming a realistic set of parameters for 
optimally-doped YBCO. The obtained Jk represents the characteristic k-dependence in 
nearly AF metals as shown in Fig. 14. Moreover, we proved the significance of the 
CVC near AF QCP beyond the one-loop (such as the FLEX) approximation in §3.4. 
As a result, it is general that Ru is prominently enhanced by the CVC when the AF 
fluctuations are strong, independent of the types of the spin fluctuation theories. 
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9.2. Applicable scope of the present study and future problems 

In this study, we explained that the CVC plays a significant roles in nearly AF metals, 
based on the microscopic Fermi liquid theory. To illustrate this theoretical idea, we 
performed numerical studies using the FLEX+CVC method in §4 - 7. Since the 
FLEX is a conserving approximation, we can perform a reliable study of the transport 
phenomena by including CVCs. As we discussed in §2.2, the FLEX approximation 
can explain characteristic electronic properties in optimally-doped HTSC. For example, 
the appropriate spin fluctuation Xqi 00 ) an d the hot/cold spot structure of quasiparticle 
damping rate 7k are reproduced satisfactorily. Although the FLEX approximation 
cannot explain the "strong pseudo-gap behavior" below T* ~ 200 K, the FLEX+T- 
matrix approximation can reproduce various transport anomalies below T* in slightly 
under-doped systems as discussed in §5. In the FLEX+T-matrix approximation, we 
take account of the strong SC fluctuations that are induced by AF fluctuations. The 
success of the FLEX+T-matrix method is strong evidence that SC fluctuations are the 
origin of strong pseudo-gap. 

Here, we discuss the applicable scope of the FLEX approximation. First, we discuss 
the "weak pseudo-gap behavior" below T ~ 600 K, where the DOS [242, 243, 244, 245], 
the Knight shift [80, 246, 247] and the uniform susceptibility [248, 249] decrease inversely 
with the growth of AF fluctuations. The weak pseudo-gap in the DOS is shallow and 
wide in energy. Although weak pseudo-gap behaviors are reproduced by the FLEX 
approximation, the obtained gap behaviors are too moderate. [Since the weak pseudo- 
gap affects transport phenomena only slightly, FLEX+CVC approximation can explain 
anomalous transport phenomena in HTSC] This failure of the FLEX will originate from 
the fully self-consistent determination of the self-energy and the Green function. In 
fact, although a large weak pseudo-gap in the DOS is reproduced at the first iteration 
stage, it vanishes in the course of iteration. Mathematically, additional self-energy 
correction introduced by iteration (i.e., the feedback effect) should be canceled by the 
vertex correction in the self-energy that is absent in the FLEX approximation (e.g., 7^ 
in Fig. 15 (ii)). The same result is reported in the GW approximation: This is a first 
principle calculation for the self-energy, which is given by the convolution of the Green 
function G and the screened interaction W within the RPA. Although the descriptions 
of the bandwidth reduction and satellite structure in the quasiparticle spectrum are 
satisfactory in a partially self-consistent GW method, they are smeared out in the fully 
self-consistent GW due to the feedback effect by iteration [250, 251, 252]. 

To produce a weak pseudo-gap successfully, it is better to derive the self-energy by 
performing (i) fully self-consistent calculation with vertex corrections or (ii) partially 
(or no) self-consistent calculation. Along the lines of method (i), Schmalian and Pines 
calculated the self-energy with all the vertex corrections by applying a high-temperature 
approximation [253]. Using a similar technique, Fujimoto calculated the pseudogap 
phenomena due to quasistatic SC fluctuations [254]. It is noteworthy that the fourth- 
order-perturbation theory with respect to U well reproduces the weak pseudo-gap 
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since the vertex corrections for the self-energy are considered [255]. This method also 
reproduces the difference of T c between YBCO and LSCO [256], and the p-wave SC 
state in Sr 2 Ru04 [257]. Along the lines of method (ii), Vilk and Tremblay proposed 
a two-particle self-consistent (TPSC) method [258], where full self-consistency is not 
imposed on the self-energy The renormalization of the spin susceptibility is described 
by introducing the effective Coulomb interaction U e s (< U) that is determined so as to 
satisfy the sum rule for Xq(w). This philosophy has been applied to extend the dynamical 
mean field theory (DMFT) by including self-energy correction due to spin fluctuations 
[259, 260]. These theories can explain the weak pseudo-gap in the DOS. However, they 
are not suitable for the study of the transport phenomena since they are not conserving 
approximations. For this reason, we use the FLEX approximation in the numerical study 
for the transport phenomena. Since the FLEX approximation tends to underestimate 
the anisotropy of t^ due to the overestimated feedback effect, the enhancements of Ru 
and Ap imp given by the FLEX+CVC method might be underestimated quantitatively. 

Finally, we present important future problems. Since the FLEX (or FLEX+T- 
matrix) approximation is a one-loop approximation, it does not work satisfactorily in 
under-doped systems. However, the analysis in §3.4 have shown that the enhancement 
of i?H due to the CVC is not an artifact of the FLEX approximation, but will be 
universal near the AF-QCP beyond the one-loop approximation. Therefore, relations 
Ru oc £^f an d Ap/Po oc cot 2 #H are expected to hold in under-doped HTSCs. It is an 
important challenge to perform a numerical study of the CVC using the higher-loop 
approximation, by taking the vertex corrections for the self-energy Then, we will be 
able to judge clearly whether the striking transport anomaly in heavily under-doped 
systems can be understood based on the Fermi liquid theory, or the idea of non-Fermi 
liquid theory explained in §1.4 has to be employed. 

In under-doped systems, small amount of impurities drastically changes the 
electronic states of the system as discussed in §7.2. Therefore, to understand transport 
phenomena in under-doped systems, we have to develop the theory of transport 
phenomena in the presence of disorder. 

In §8.1, we presented a numerical study of the transport coefficients in the quasi 
2D and 3D Hubbard models using the FLEX+CVC approximation, and explained the 
experimental results in CeMIns. However, in real heavy fermion systems, the Fermi 
surfaces is much more complicated. As discussed in §8.1, the huge Ru (^> 1/ne) in 
CeMIns cannot be accounted for a simple RTA in the multiband system. However, 
we have to study the multiband system for a more realistic study. In fact, we recently 
studied the 2D and 3D periodic Anderson model, which is a two-band model for a heavy 
fermion system, and found that a relatively large mass-enhancement factor z~ l ~ 10 is 
obtained in the FLEX approximation [237]. Since p, R^ and Ap/p are independent 
of z -1 , we could discuss their critical behaviors using the Hubbard model as we have 
explained in §8.1. Since S and v are proportional to z _1 , on the other hand, we have 
to take account of the experimental z~ x if one compare the theoretical results with 
experimental values. 
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Finally, we comment on the "field-induced QCP" in CeCoIn 5 near H c2 ~ 5 
Tesla, where prominent increment in the effective mass m* was inferred experimentally. 
Paglione et al. [234] measured the in-plane resistivity near H c2 , and derived the relation 
A oc (H — i? c2 )~ 4 / 3 , where A is the coefficient of the T 2 -term in the resistivity. According 
to the Kadowaki- Woods relation yA oc m* (see eq. (164)), m* diverges at H = H c2 . 
This fact strongly suggests that the field-induced QCP occurs at H = H c2 = 5 Tesla 
in CeCoIn 5 . The prominent increment of m* around the field-induced QCP is also 
confirmed by thermoelectric transport measurements (S and v) [261]. Interestingly, 
a similar field-induced QCP is also realized in heavily over-doped Bi2212 (T c = 10 K) 
[262] . Up to now, the origin of the critical behavior is unknown. One possible origin will 
be the (field-induced) SDW state that is hidden in the SC state. In CeCoIn 5 , Tanatar et 
al. [263] measured the anisotropy of the resistivity near H c2 , and found that the in-plane 
resistivity (p a ) is linear-in-T above 4 K, whereas it decreases more quickly below 4 K. 
The observed T-dependence of p a looks similar to the theoretical result in Fig. 34 (i). 
In contrast, inter-plane resistivity (p c ) shows a complete T-linear dependence from 25 
mK to 16 K, which might suggest that the CVC is unimportant for p c . Moreover, the 
Wiedemann- Frantz law (n/aT = (7r 2 /3)(/cB/e) 2 ; k being the thermal conductivity) is 
strongly violated only for the inter-plane direction. This strong anisotropy in transport 
may be difficult to explain by the current theory. Since the transport phenomena under 
high magnetic field is outside of the scope of the present study, this is an important 
future problem of transport phenomena. 

Up to now, electronic properties in strongly correlated metals under high magnetic 
field have not been studied sufficiently. This issue is an important future challenge. We 
should note that the if C 2-line in CeCoIns is a weak first-order line for T <C T c q [106]. 
Therefore, exactly speaking, the field-induced QCP in CeCoIn 5 cannot be a true QCP. 

9. 3. Fermi arc picture and transport phenomena 

According to the ARPES measurements for under-doped Bi2212 [85, 264, 265, 266] and 
LSCO [267], the intensity of the spectrum Pk(w) = ImG^uj — i5)/n on the Fermi surface 
at uj = takes very small value around the hot spot below the pseudo-gap temperature 
T*. As a result, pk(0) has finite intensity only around the cold spot, which is frequently 
called the "Fermi arc". The length of the Fermi arc diminishes as the temperature 
decreases, particularly below T*. References [85, 264, 265] have shown that the pseudo- 
gap is a precursor of the SC gap. Based on this observation, T-matrix theory [87] 
and the FLEX+T-matrix theory [51, 52, 88, 89] have been proposed. Moreover, recent 
ARPES measurements [268] suggests that the reduction in Pk(0) around the hot spot 
in the Fermi arc state is not due to a "density order" . 

The temperature dependence of the Fermi arc in under-doped systems is 
qualitatively reproduced by the FLEX+T-matrix approximation, where no hidden 
order is assumed. Figure 37 shows pk(0) given by the FLEX and FLEX+T-matrix 
approximations for LSCO (n = 0.9) at T = 0.02 ~ 80 K. The intensity of Pk(0) at around 
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Figure 37. (i) Quasiparticle spectrum pk(0) = ImGk( — iS)/w for LSCO at T = 0.02 
given by the FLEX approximation. pk(0) is shown only when pk(0) > 0.1. The contour 
is shown on the basal plane, (ii) Pk(0) given by the FLEX+T-matrix approximation. 
The contour shows a "Fermi arc structure" . 



(it, 0) is suppressed in the FLEX approximation. The Fermi arc structure becomes more 
prominent in the FLEX+T-matrix approximation, because the quasiparticle damping 
rate due to strong <i x 2_ J/ 2-wave SC fluctuations, which is given in eq. (115) in the 
FLEX+T-matrix approximation, takes large values around (tt, 0). In §5, we explained 
that anomalous transport phenomena below T* are well explained by the FLEX+T- 
matrix approximation, by considering the CVC due to AF+SC fluctuations. Even 
above T*, the weak pseudo-gap due to the AF fluctuations is widely observed by ARPES 
measurements, which reflects the hot/cold-spot structure of r^. Therefore, a kind of the 
Fermi arc structure exists even above T*. 

The Fermi arc structure is also observed in under-doped NCCO: In the ARPES 
measurements [70], the intensity of Pk(0) for a portion of the Fermi surface near the 
AFBZ (see Fig. 4) is strongly suppressed below Tn ~ 150 K by the SDW order. Below 
Tn, i?H is negative by reflecting the small electron- like Fermi surface around (tt, 0), 
which is caused by the reconstruction of the Fermi surface due to the SDW order. 
Surprisingly, R^ remains negative even above Tn, although SDW-induced Fermi surface 
reconstruction is absent. Moreover, no anomaly in i?H is observed at Tsf. This highly 
nontrivial fact is explained by the Fermi liquid theory by considering the CVC, as shown 
in Fig. 19 (iv). According to the ARPES measurements [269], pseudo-gap state due to 
AF fluctuations is observed even above Tn- Therefore, a Fermi arc state appears to be 
realized in electron-doped systems. 

When a distinct Fermi arc structure exists, the "effective carrier density at the Fermi 
level (n e ff)" reduces in proportion to the length of the Fermi arc. Since -Rh oc l/en e s 
in the RTA, it is sometimes claimed that the enhancement of i?n m HTSCs can be 
explained by the reduction in n e & due to the emergence of the Fermi arc. This idea 
seems to be valid in heavily under-doped compounds [270]. However, this idea is not 
true for slightly under-doped systems since Ru starts to decrease below T*, whereas 
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Figure 38. Schematic picture of the "Fermi arc structure" in a hole-dopes system. 
The "effective Fermi surface (FS)" above T*, which was explained in Fig. 14 (ii), 
is also described. The smallness of the effective Fermi surface is the origin of the 
enhancement of i?H above T* . Below T*, the effective Fermi surface approaches the 
true Fermi surface since the CVC due to AF fluctuations are reduces. Therefore, i?n 
decreases below T*. 



the length of the Fermi arc shrinks due to SC fluctuations. The reduction in Rn in the 
pseudo-gap region is explained by the reduction in the AF fluctuation, as explained in 
§5.1. 

Here, we discuss the Hall coefficient in the presence of the "Fermi arc structure due 
to strong AF fluctuations" , by taking the CVC into consideration. The cases we consider 
correspond to the hole-doped systems above T*, and the electron- doped systems above 
T N . Figure 38 shows the Fermi arc structure due to the weak pseudo-gap formation and 
the "effective Fermi surface" that is defined to be perpendicular to the total current 
Jk, as introduced in Fig. 14 (ii). (Note that the "effective Fermi surface" is not visible 
in the ARPES measurements.) Within the RTA, a^J A oc Z c /7coid and <x^J A oc l c /lcoidi 
where l c is the length of the cold spot (Fermi arc). Therefore, 

(149) 



i?l TA oc l/f c . 



As we have discussed, a xy is proportional to the curvature of the effective Fermi surface 
at the cold spot. In the presence of strong AF fluctuations above T*, 9£ rotates 
approximately by it around the cold spot, as shown in Fig. 38. Therefore, the curvature 
of the effective Fermi surface at the cold spot in Fig. 38 is proportional to l/l c . As a 
result, due to the CVC above T*, 

R R oc l/l 2 c (150) 

approximately. The large curvature of the effective Fermi surface around the cold spot 
is the origin of the enhancement of R^. In conclusion, when the Fermi arc emerges, the 
magnitude of Ru is much larger than the value obtained by the RTA that is given in 
eq. (149). 
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We note that if the distinct Fermi arc structure due to the highly anisotropic r^ 
is realized, it should lead to a huge magnetoresistance within the RTA: in fact, Aa xx 
diverges when the k-dependence of the mean free path /k = fk^k is not continuous on 
the Fermi surface; in this case, the second term in eq. (99) diverges as pointed out in 
Refs. [26] . Therefore, the modified Kohler's rule given by eq. (4) cannot be satisfied. In 
summary, a unified understanding of anomalous transport phenomena in HTSC cannot 
be achieved by the RTA even if the Fermi arc structure is taken into account. 

9.4- Unconventional transport phenomena in multiorbital systems: anomalous Hall 
effect and grand Kadowaki- Woods relation 

In previous sections, we studied single-band models, and showed that the CVC induces 
various striking non- Fermi-liquid-like behaviors in the presence of AF fluctuations. For 
example, Ru shows strong temperature dependence due to the CVC. In multiband 
systems, however, Rn can exhibit temperature dependence or sign change within the 
RTA, if a hole-like Fermi surface and an electron-like Fermi surface coexist and their 
relaxation times have different T-dependences. There are many such examples even in 
conventional metals. In this study, we did not discuss the multiband effect, since the 
magnitude \Rn\ will remain ~ 1/ne in this mechanism. Thus, the multiband mechanism 
is impossible to explain the huge \Rn\ observed in various systems near the AF QCP. 

However, a kind of multiorbital effect causes a huge Hall coefficient in d- and /- 
electron systems, which is known as the "anomalous Hall effect (AHE)". In the presence 
of the AHE, the Hall resistivity is given by 

p H = R n H B + R a H M, (151) 

where B is the magnetic field and R\ is the ordinary Hall coefficient. R\ is the 
anomalous Hall coefficient, which is (generally) proportional to the magnetization M. 
In general, R\ and R\ are even functions of B and M. In ferromagnets, the second 
term R & H M takes a finite value even if B — 0. Study of the AHE due to multiband 
effect was initiated by Karplus and Luttinger [271], and by Luttinger [272]. They found 
that the anomalous Hall conductivity (AHC) a% y (= R^M/p 2 ) is independent of the 
resistivity p. This Karplus-Luttinger-term is called the "intrinsic AHE" because it exists 
even in systems without impurities. Later, Smit presented a mechanism of "extrinsic 
AHE" [273]: he found that the spin-polarized electrons are scattered asymmetrically 
around an impurity site in the presence of spin-orbit interaction. The AHC due to this 
skew-scattering mechanism is linearly proportional to p. 

In paramagnetic heavy-fermion systems, the observed Hall coefficient is given by 

R* F =p H /B = R n H + R* I RE , (152) 

where i?^ HE = R^ ■ x and x = M/B is the uniform magnetic susceptibility. In heavy- 
fermion systems, the AHE due to the second term of eq. (152) takes large value since 
the uniform susceptibility M/B = x is widely enhanced due to the strong Coulomb 
interaction. Due to the AHE, R^f starts to increase with increasing temperature from 
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OK, and it turns to decrease after exhibiting its maximum value around the coherent 
temperature T co ^. The maximum value of \R^ F \ is more than one order of magnitude 
greater than l/\ne\ ~ 1.0 x 10~ 9 m 3 /C. The AHE cannot be derived from the RTA 
since the interband hopping of electrons is important. [In the AHE, the effect of CVC 
is not expected to be crucial.] Experimentally, R% is positive in usual Ce and U based 
heavy-fermion systems [274, 275]. 

Here, we discuss the AHE in heavy-fermion systems. Bandstructure of Ce- and Yb- 
based heavy-fermion system is described by the following orbitally degenerate periodic 
Anderson model [276, 277]: 

H= ^ekcl a c ka + ^2E f fl M f kM +^2(V£ Ma fl M c ka + h.c.) 

kcr kM kM<r 

+ ~2 ^ /k+qJVf/k'-qM'/k'M'/kM, (153) 

kk'qMM' 

where cj^ (/j[ M ) is the creation operator of the conduction electron (/-electron) 
with a — ±1 (M = J, J — 1, •••,— J). /-orbital degeneracy is Nf = 2J + 
1. ek (Ef) is the spectrum for conduction electrons (/ electrons). In the 
case of Ce-compound (J = 5/2), the complex c-f mixing potential is given 
by IW = o-\/ v / 4V3 V / (7/2 - 2Ma)/7Y l ^ 3 ' (7 (e k , <p k ), where Y^ 3 (6 k , Vk ) is the 
spherical harmonic function, In the case of Yb-compound (J = 7/2), V-^Mo — 
VoV^Vi 7 / 2 + 2Ma)/7YM-°(e k , <p k ). Note that the relation Em=-j |VkM*| 2 = V 2 
holds. 

The Green function for c electron is given by [276, 277] 

Gl(u) = ( uj + fj, - e k 



uj + /i — Ef — S(co>) 



u; + /x-e k {V J )2 , ) , (154) 

uj — Ef — ry* I 

where \/ Q * = v^F , #/ = z(E f + ReS(O) - //), 7* = zlmE(O), and s = (1 - 
(9ReS(a;)/(9c<j) _1 | w= o is the renormalization factor. Since the quasiparticle energy E£ 
satisfies Re{l/G k (-E^)} = 0, the equation for _E£ is given by 

^ + / ,_ 6k _ W) 2 [K-E f ) ^ Q (155) 

k (^ k -^/) 2 + (r) 2 

Here, we derive the quasiparticle energy by analyzing eq. (155): At zero temperature 
where 7* = 0, heavy quasiparticles band is formed due to c-f hybridization. Below T co h, 
the quasiparticles band is 



K± = Ue 1t -n + E f ± y^k - ~n - EfY + A(V *yj , (156) 

where ££_ (££ + ) represents the lower (upper) quasiparticles band. 7* grows 
monotonically as temperature increases, and \Ef\ ~ 7* at the coherent temperature 
Tcoh- When T ^> T coh , c-f hybridization is prohibited by 7* ^> \Ef\. 
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Figure 39. i?H in U3Ni3Sii4 plotted as functions of (i) T 2 and (ii) p 2 . The relation 
K& oc p 2 is satisfied well, as predicted by the theory of AHE [277]. [Rcf. [278]] 



RHF 



We summarize the electronic states in heavy-fermion systems: (i) Below T co h, 
Fermi liquid state with heavy quasiparticles is realized due to c-f hybridization. 
Mass enhancement factor z~ l and uniform susceptibility \ °^ z ~ 1 are constant. The 
heavy quasiparticle bandwidth W-^f is approximately given by Hhf ~ min{_Ej^ + — 
-££_} ~ (Vo) 2 /W c ~ \Ef\, where W c is the c-electron bandwidth, (ii) Above T coh , 
c-f hybridization ceases, and localized /-electrons causes Curie-Weiss susceptibility. 
Experimentally, the temperature of maximum resistivity T* is larger than T co h. 

In paramagnetic heavy-fermion systems, the Hall coefficients takes huge values due 
to the AHE [274, 275, 278]. In the early stage, Coleman et al [279] and Fert and Levy 
[280] developed theories of extrinsic AHE: they studied the extrinsic mechanism based 
on the /-electron impurity Anderson models with rf-orbital channels, and predicted the 
relation i?^ HE oc XP above T co h when the ci-orbital phase-shift is finite. 

On the other hand, Kontani and Yamada studied the intrinsic AHE due to 
multiband effect based on the J = 5/2 periodic Anderson model, which a the model 
in the strong limit of spin-orbit interaction [277, 281]. In this model, the anomalous 
velocity originates from the angular momentum of the localized /-electrons. That 
is, the AHE originates from the transfer of angular momentum of the /-electron to 
the conduction electron. Based on the linear-response theory, they derived a large 
intrinsic AHC: a% y /H[l/Q cm Gauss] = 5.8 x HT^cm]- 3 ) 1 / 3 / E f [K] where E f is 
the renormalized /-level measured from the Fermi energy. \Ef\ gives approximate 
quasiparticle bandwidth. The present study predicts that i?^ HE > (< 0) in Ce (Yb) 
based heavy-fermion systems since Ef > (< 0). This prediction is consistent with 
the experimental results [274, 275]. If we substitute Ef = 10K and n = 10 22 cm~ 3 , 
o~®/H = 1.3 x 10 2 [l/fi cm Gauss], which is a typical experimental value. The predicted 



temperature dependence of Ry\ 



AHE 



[cr 
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is R*m oc ( 7 *7(£ 2 f + l* 2 ) X [277]. 
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Therefore, Rjf LE shows the following cross-over behavior: 

^ HE oc XP 2 [oc p 2 ] : below T coh (7* < \E f \), (157) 

< HE oc X [oc 1/T] : above T coh (7* > \E f \). (158) 

In typical Ce-based heavy-fermion systems, \j z ~ 100 and \Ef\ ~ 10 K. Below T co h, eq. 
(157) is proportional to p 2 since x is constant for T < T co h. Figure 39 shows the i?n in 
U 3 Ni 3 Sn 4 (7 = 380mJ/K 2 ), where the relation R^ HE = R^ F -i2g oc p 2 holds well below 
~ 0.3T co h ~ 25K. Note that the extrinsic type AHE [279, 280] and the intrinsic type 
AHE [277] can coexist. However, the relation i?^ HE oc p 2 suggests that the intrinsic-type 
AHE is dominant at least below T co h. Above T co h, the AHE due to interband transition 
is suppressed when c-f mixing is prohibited. Therefore, eq. (158) is independent of p. 
This coherent-incoherent crossover of intrinsic AHE was first theoretically derived in 
Ref. [277]. This crossover behavior is not restricted to heavy-fermions systems, but also 
observed in various transition ferromagnets [282]. 

In heavy-fermion systems near the AF QCP, i?g shows strong temperature 
dependence due to the CVC. To extract i?g from the observed Hall coefficient, we have 
to seriously consider the T-dependence of R^ RE . Very fortunately, AHE is vanishingly 
small in CeMIns [19, 20]. Therefore, we could perform reliable analysis of the ordinary 
Hall effect in CeMIn 5 as discussed in §8.1, without the necessity of subtracting the AHE. 
Note that the AHE vanishes when the crystal-field splitting of the /-levels is much larger 
than T co h. This may be the reason for the small AHE in CeMIn 5 . 

Recently, theory of the intrinsic AHE in (ferromagnetic) <i-electron systems has 
been developed based on realistic multiorbital tight-binding models [283, 284]. It was 
revealed that a large anomalous velocity emerges in general multiorbital rf-electron 
systems because of the inter-orbital hopping. This is the origin of the large AHE in 
c?-electron systems, which is very similar to the origin of the AHE in /-electron systems 
[277]. The intrinsic AHEs for Fe [286] and in SrRuC>3 [287] were also calculated based 
on the LDA band calculations. 

Here, we present an intuitive explanation of the AHE due to the multiorbital 
effect [284, 288]. In a multiorbital system, a conduction electron acquires the "effective 
Aharonov-Bohm phase factor" due to rf-atomic angular momentum with the aid of the 
spin-orbit interaction and the inter-orbital hoppings, which is responsible for the Hall 
effect. The intuitive explanation based on the tight-binding model, which is given in 
Fig. 40, which represents a square lattice (d xz , G^-orbital tight-binding model with 
^-component of the spin-orbit interaction ^2^1 z ■ s 2 )j [284], which is a simplified model 
for (Ca,Sr) 2 Ru04: This compound shows large AHE under the magnetic field [285]. 
Its magnitude is comparable with the large AHE in /-electron systems (such as UPts). 
In Fig. 40, ±t' represents the interorbital hopping integral between next nearest sites. 
Now, let us consider the motion of a down-spin electron along a triangle of half unit 
cell: An electron in the cL^-orbital can transfer to c^-orbital and vise versa using the 
spin-orbit interaction ±h~\l z , where (yz\l z \xz) = — (xz\l z \yz) = i. By combination of 
angular dependence of interorbital hopping and the spin-orbit interaction, any clockwise 
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Figure 40. Effective magnetic flux for down-electrons in a square lattice (d xz ,d yz )- 
orbital tight-binding model with z-component of the spin-orbit interaction "Y^Alz ■ s z )i. 
The effective magnetic flux for up-electrons is negative. Both \xz) and \yz) are given 
by the linear combination of \l z = ±1). Since l 2 z = 1 in the present basis, the rotation 
operator R(ff) = e * = cos6 — il z s'm0 is equal to —il z for 9 = tt/2. Therefore, 
l z \xz) = iR(jr/2)\xz) = i\yz) and l z \yz) = iR(n/2)\yz) — —i\xz). 



(anti-clockwise) motion along any triangle path with spin-orbit interaction causes the 
factor +i (— i). This factor can be interpreted as the "Aharonov-Bohm phase factor " 
exp(27ri0/0 o ) [</>o = ^ C /M]; where 4> = 0o/4 represents the "effective magnetic flux" in 
the half unit cell. We revealed the fact that the effective magnetic flux, which is inherent 
in multiorbital <i-electron systems, causes huge AHE in various transition metals. 

In many <i-electron ferromagnets, intrinsic AHE (R a H oc p 2 ) seems to be observed 
experimentally. In particular, Asamitsu et al. measured AHE in various transition- 
metal ferromagnets, and found out a universal crossover behavior from R ^ oc p 2 to 
R a H oc p n (n = ~ 0.5) as the resistivity increases, around p ~ 100 pficm. Recent 
theoretical study in Ref. [284] has revealed that the experimental result is explained by 
the following coherent-incoherent crossover 

R a H ocp 2 : 7 < A, (159) 

R a H oc const. : 7 > A, (160) 

where A represents the minimum bandsplitting near the Fermi level. A > 0.1 eV in usual 
transition metal ferromagnets. It is apparent that eqs. (159) and (160) correspond to 
eqs. (157) and (158) for heavy-fermion systems. [Note that R^ nE = -Rh'X-] Therefore, a 
conventional behavior R a H oc p 2 is violated in high-resistivity metals since the interband 
particle-hole excitation, which is the origin of the AHE, is suppressed when 7 is larger 
than bandsplitting energy A. 

There is a simple intuitive explanation for this coherent-incoherent crossover of 
intrinsic AHE: When 7 is sufficiently small, the intrinsic Hall conductivity is proportional 
to the lifetime of the interband particle-hole excitation: h/A [271, 272, 277]. In the 
high-resistivity regime where 7 ^> A, the SHC decreases drastically with 7 since the 
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interband excitation is suppressed when the quasiparticle lifetime /1/7 is shorter than 
h/A. According to the above discussion, coherent-incoherent crossover of intrinsic AHE 
should be universal and is widely observed in various multiorbital p-, d-, and /-electron 
systems. 

It is noteworthy that the large anomalous velocity also induces a sizable spin Hall 
effect (SHE) in paramagnetic multiorbital systems [284, 288, 289]. The SHE is the 
phenomena that an applied electric field induces a spin current j s = j-\—Ji in a transverse 
direction. Recently, SHE attracts great attention due to its fundamental interest and its 
potential application in spintronics. Karplus and Luttinger [271] showed that an applied 
electric field induces a spin-dependent transverse current in the presence of spin-orbit 
interaction. This mechanism causes the AHE in ferromagnetic metals, and the SHE 
in paramagnetic metals. Recently, the spin Hall conductivity (SHC) in Pt had been 
observed [290], and found that the SHC in Pt is 10 4 times larger than that observed 
in n-type semiconductors. This experiment fact has attracted considerable attention to 
the study of the SHE in transition metals. 

Recently, we propose a new mechanism for the giant SHE originate from the d- 
orbital degrees of freedom, which is absent in semiconductors [288, 289, 291]. In a 
multiorbital system, a conduction electron acquires the "effective Aharonov-Bohm phase 
factor" due to d-atomic angular momentum, as explained in Fig. 40. We have studied 
the SHEs in Sr2Ru04 [288], which is the first theoretical study of the SHE in <i-electron 
systems, in Pt [289], and in various Ad and 5d transition metals [291]. It is found that the 
explanation in Fig. 40 seems to capture the characteristics of SHE in rf-electron systems. 
In contrast, the "Dirac monopole mechanism" [292] is appropriate when massless Dirac 
corn dispersion exists close to the Fermi level like in semiconductors. The present 
mechanism of SHE is also different from that in Rashba-type 2D electron-gas model due 
to momentum- dependent SOI [293]. Thus, giant SHE due to atomic orbital degrees of 
freedom is ubiquitous in various p-, d-, /-electron systems. 

Finally, we discuss the multiorbital effect on the Kadowaki- Woods (KW) ratio 
A'-/' 2 , where A is the coefficient of the T 2 term in the resistivity, and 7 is the 
coefficient of the T-linear term of the electric specific heat. Experimentally, the KW 
ratio in Ce- and U-based heavy-fermion compounds shows an approximate universal 
value A'-f^ 2 rs 10~ 5 [fiQ cm (mol-K/mJ) 2 ], which is known as the KW relation [294]. 
Although it was believed to be universal in heavy-fermion systems for a long time, recent 
experimental activities have revealed that the KW relation is strongly violated in many 
Yb-based compounds. Recently, the author had derived a generalized KW relation 
that is applicable for systems with general /-orbital degeneracy Nf for Ce- and Yb- 
based compounds [295] and for Sm- and Er-based compounds [296]. By considering the 
material dependence of Nf, the mystery of the failure of the KW relation was resolved. 

Here, we analyze the KW ratio in terms of the DMFT (d = oo-limit) [109, 110], 
which is believed to be useful in heavy-fermion systems that is not close to the AF- 
QCP. By using the DMFT, we can utilize the strong-coupling Fermi liquid theory for 
the impurity Anderson model developed by Yamada and Yosida [297] : Using the Ward 
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Figure 41. (a) Grand Kadowaki- Woods relation in Rcf. [299], where A = 
2A/N f (N f - 1) and 7 = 2^/N f (N f - 1). (b) A scaling relation between E mil 

and yu/770 [305]. -Emir corresponds to 2V^, and ya/TTo <x ^/zT^ (\/o is a constant of 
order 0(1)). [Ref. [306]] 



identity, they derived the exact Wilson ratio R = {x/x°)/{l/l°) = 2 for Nf = 2 in the 
strong coupling limit, that is, z^ 1 ^> 1 and the /-electron charge susceptibility is zero 
(Kondo model limit). Here, x° an d 7° represent the non- interacting values. Using the 
DMFT, 7 and ImE(0) in the PAM (in eq. (153) is given by [295, 298] 

7 = N A k 2 B ^-N f (N f - l)r loc (0, 0)p / (0) 2 , (161) 



ImS(0) 



N A k 2 B -N f (N f - l)T loc (0,0)p^0y 
^ B T) 2 /Ar ^ 2 , ftn , ,, 



(^-1)^(0,0^(0)^, 



(162) 



in the strong coupling limit. Here, N A = 6.02 x 10 23 is the Avogadro constant, and 
rioc(0,0) is the local four-point vertex. p-^(O) is the DOS for /-electron per channel; 
Nfpf(0) is the total DOS at the Fermi level. When ek in eq. (153) is a free dispersion, 
the conductivity is given by 

_e 2 (3/ 7 r) 1 AV/ 3 a 3 
° hNpf(0)-hnZ(0y [ } 

where n(= kp/3n 2 ) is the density of quasiparticles that form the conduction band, and 
a is the unit cell length. According to eqs. (161)-(163), we find the scaling properties 
7 oc Nf(Nf — l)ri oc p^ and A oc Nf(Nf — l)T 2 oc pf . Therefore, we obtain the "grand 
KW relation" that is valid for any N f (> 2) [295, 296, 298, 



A-f 



-2 



h 9(3 



TV 



1-1/3 



e 2 kl 4n 4 / 3 a 3 iV| \N f {N f - 1) 



'B 

1 X 10~ 5 



\Nf{N f 



[/j,fi cm (mol-K/mJ)" 



(164) 
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where we put h/e 2 = 2.6 x 10 4 fi, k B = 1.38 x 1CT 23 JK" 1 , and we assumed 
l/n 4//3 a 3 k4x 10~ 8 cm. 

Without crystal field splitting of the /-level, Nf = 2 J + 1 = 6 for Ce 3+ and Sm 3+ 
ions, and Nf = 8 for Yb 3+ and Er 2+ ions. Therefore, the previous KW relation turned 
out to be valid only when Nf = 2 (Kramers doublet case due to strong crystal field 
splitting). A similar universal relation lim^^o eS/T'-f f» ±1 (S is the Seebeck coefficient) 
was recently found [300, 301]. The characteristics of the electronic state in heavy- 
fermion systems are (i) large mass-enhancement and (ii) small charge susceptibility 
since the /-electron is almost localized. The grand KW relation (164) is derived only by 
imposing these constraints on the microscopic Fermi liquid theory. This fact illustrates 
a remarkable advantage of the Fermi liquid theory for the analysis of strongly correlated 
systems. 

Figure 41 (i) shows the experimental verification of eq. (164) for various heavy- 
fermion compounds, where Nf in each compound was determined by the temperature 
dependence of x an d the inelastic neutron scattering [299]. Tsujii's study confirmed 
that Nf ~ 2 in many Ce-based compounds, where crystal-field splitting is larger than 
the renormalized Fermi energy W^hf- On the other hand, Nf ~ 8 in many Yb-based 
ones, where crystal-field splitting is smaller than W^hf- Torikachvili et al. also found 
that other Yb-based heavy-fermion systems YbT 2 Zn 2 o (T = Fe, Co, Ru, Rh, Os, or 
Ir) follows the grand KW relation shown in eq. (164) [302]. When Nf is larger than 
the number of conduction band N c , almost unhybridized /-levels exist near the Fermi 
level, except when the crystal-field splitting is very large. This situation induces a huge 
anomalous Hall coefficient [277] as well as the large Van Vleck susceptibilities [303, 304] 
in Kondo insulators and in singlet SC heavy-fermion systems. 

We shortly discuss the absence of the multiband effect on the KW relation: In a 
usual heavy-fermion compound, there are one or two (relatively) large Fermi surfaces 
composed of heavy quasiparticles, and several (relatively) small Fermi surfaces composed 
of light-quasiparticles. In the presence impurities, AT 2 oc (7k) fs as explained in eq. 
(33). Therefore, heavy quasiparticles on the large Fermi surfaces give the dominant 
contribution to both the specific heat and the A-term. As a result, the grand-KW 
relation is universally realized even in multiband systems. Moreover, in many heavy- 
fermions, there is evidence that the heavy quasiparticles on the large Fermi surfaces give 
the dominant contribution to the conductivity. For example, |i?g| ~ 1/ne seems to be 
realized in many heavy-fermions away from AF-QCPs, although |i?g| ^> 1/ne should 
be realized when small Fermi surfaces are the most conductive. 

Figure 41 (ii) shows the mid-infrared peak energies E m \ r of the optical conductivity 
o~(u) in various heavy-fermion systems [305]. According to the Fermi liquid theory, 
E miT w \fzV is satisfied in the PAM, where V is the c-f mixing potential and 
1/z = m*/m is the mass enhancement factor. Okamura confirmed that the relation 
E m i r oc 7~ x / 2 is universally satisfied in various heavy-fermion compounds, as shown 
in Fig. 41 (ii). This scaling relation is well satisfied regardless of Nf. Their study 
established the validity of the PAM and the Fermi liquid theory in various heavy-fermion 
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compounds. 
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